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USER'S  MANUAL  FOR 
MIXING  ZONE  ANALYSIS  PROGRAMS 


SUMMARY 

This  User's  Manual  describes  the  calibration  and  application  of  a  personal 
computer  package  to  predict  the  mixing  zones  in  a  shallow  river  for  point  and 
diffuser  discharges.  The  package  is  set  up  in  an  interactive  (enquiry/response) 
mode.  The  required  site-specific  data  for  the  package  are  described.  The 
package  also  predicts  the  critical  points  in  any  river  transect  where  the 
provincial  water  quality  objectives  are  achieved.  The  package  outputs  include 
many  computer  graphic  options  to  assist  the  user.  In  the  appendices,  detailed 
technical  discussions  are  presented  on  the  various  package  components. 
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LOGICIELS  D' ANALYSE  DES  ZONES  DE  MELANGE 
MAITUEL  DE  L'UTILISATEUR 

SOMHAIRE 

Le  present  manuel  de  1 'utilisateur  montre  comment  etalonner  et 
utiliser  un  progiciel  servant  a  etablir  les  zones  de  melange 
dans  une  riviere  peu  profonde  en  cas  de  deversements  ponctuels 
et  de  deversements  diffus  de  polluants.  Le  progiciel 
fonctionne  en  mode  dialogue  (question-reponse) .  On  enonce  les 
donnees  particulieres  relatives  a  1 ' emplacement  etudie,  dont  on 
a  besoin  pour  utiliser  le  progiciel.  Le  progiciel  permet 
egalement  de  determiner  les  criteres  precis  selon  lesquels  une 
section  de  cours  d'eau  peut  repondre  auz  objectifs  de  qualite 
de  1 ' eau  de  la  province.  L' utilisateur  pourra  se  servir  de 
nombreuz  graphiques  figurant  parmi  les  etats  que  peut  creer  le 
progiciel.  Une  description  technique  des  elements  du  progiciel 
figure  en  ajineze. 


MIXING  ZONE  ANALYSIS  PROGRAMS 


INTRODUCTION 

Water  pollution  control  plant  effluents  discharged  into  receiving  streams  and 
rivers  often  contain  substances  such  as  chlorine  and  ammonia  which  are 
potentially  toxic  to  aquatic  biota.  Provincial  Water  Quality  Objectives  (PWQO) 
require  the  maintenance  of  a  portion  of  the  river  as  a  favourable  habitat  for 
the  biota  termed  the  "zone  of  passage"  (ZOP),  wherein  the  concentrations  of  the 
pollutants  comply  with  a  specified  water  quality  objective  (Cs).  The  remaining 
portion  of  the  river  where  the  pollutant  concentrations  do  not  con^ly  with  the 
specified  objective  is  called  the  limited  use  zone  (LUZ).  The  salient  features- 
of  the  mixing  zone  are  depicted  in  Figure  1.  Effluents  are  usually  discharged, 
through  outfalls  at  the  river  bank  (i.e.,  point  sources  located  at  the  shore), 
especially  in  shallow  rivers,  although  diffuser  outfalls  are  used  in  larger 
rivers. 


This  manual  discusses  the  integrated  use  of  a  number  of  programs  based  on  a 
steady-state  mathematical  model  to  predict  the  spatial  concentration 
distributions  in  the  mixing  zones  of  shallow  rivers.  The  model  is  based  on  the 
two-dimensional  convective  dispersion  equation  and  utilizes  the  stream  tube 
co-ordinate  transform  concept  developed  by  Yotsukura  &  Cobb  (1972).  In  this 
package,  four  programs,  developed  by  T.P.H.  Gowda  are  utilized.  These  programs 
have  been  combined  with  graphical  output  programs  and  interactive  inputs  for  use 
by  investigators  of  river  water  quality. 
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Figure  1   Schematic  Diagram  of  Mixing  Zone 
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FIELD  STUDIES  AND  DATA  COLLECTION 


The  collection  of  the  required  field  data  is  important  to  the  model  predictions. 
In  order  to  model  the  distributions  of  the  effluent  plume  in  a  receiving  river 
downstream  from  an  outfall,  measurements  of  salient  features  of  the  river  must 
be  made  at  several  river  cross-sections  or  transects  at  different  distances  from 
the  outfall  (Figure  2). 


Outfall 


Figure  2  Position  of  Transects 


For  this  package,  measurements  should  be  made  at  4  to  8  transects  downstream  of 
the  outfall.  The  position  of  the  first  transect  should  be  far  enough  downstream 
from  the  source  that  concentration  measurements  reflect  the  behavior  of  the  "far 
field"  mixing  regime  or  at  least  the  further  parts  of  the  "intermediate  field". 
The  terms  "near  field",  "intermediate  field"  and  "far  field"  are  used  to 
differentiate  mixing  regimes  in  the  mixing  zone  which  are  characterized  by  their 
time  rates  of  growth  of  variance  of  contaminant  clouds  or  plumes. 

The  "near  field"  mixing  processes  are  dominated  by  jetting  effects.  The 
"intermediate  field"  is  the  region  where  the  width  of  the  plume  is  smaller  than 
the  largest  turbulent  motions  present  in  a  receiving  water  (river,  lake  etc). 
Here  the  growth  rate  of  the  variance  of  the  plume  (dispersion)  is  larger  than 
that  of  the  "far  field"  (See  figure  3).  The  "far  field"  mixing  processes  are 
effected  by  all  sizes  of  turbulent  gyres  present  in  the  receiving  water  (ambient 
mixing  processes) . 
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Figure  3  Variance  versus  Time 

In  the  "far  field"  we  can  see  that  the  dispersion  D  =  1/2  do^/dt  is  a  constant. 
In  the  "intermediate  field"  the  dispersion  D  is  larger  than  in  the  far  field. 


Data  Required  at  Each  Transect 

The  data  used  to  drive  this  model  should  be  in  the  following  units. 

i)  distances  in  meters 

ii)  flow  rates  in  cubic  meters/second 

iii)  velocities  in  meters/second 

iv)  temperature  in  Celsius  degrees 

Concentration  may  be  expressed  in  any  unit  desired  since  the  equations  are 
linear  the  output  units  will  be  the  same  as  the  input  units. 

The  following  measurements  are  necessary. 


1.  Distance  downstream  from  outfall  in  meters. 

2.  Measurements  at  10  to  25  points  across  the  transect  of: 

i)  distance  y  from  the  reference  bank 

ii)  depth  at  y 

iii)  depth  averaged  velocity  at  y  (optional) 

iv)  concentration  of  pollutant  at  y 


Figure  4 


These  measurements  must  be  made  at  the  same  point  at  a  distance  y  from  the  bank. 
The  model  assumes  average  depths,  therefore  parameters  such  as  velocity  or 
concentration  must  be  taken  at  several  depths  at  each  point  and  averaged. 

Further  river  data  required  are: 

1.  The  contaminant  being  measured. 

2.  Flow  rate  of  the  river  upstream  of  the  outfall. 

3.  Flow  rate  of  the  effluent  from  the  outfall. 

4.  Background  concentration  of  the  contaminant. 
5.'  Effluent  concentration  of  the  contaminant. 

6.  The  decay  rate  of  the  contaminant  and  the  temperature  at  which  this  rate 

applies. 

7.  Temperature  of  the  river. 

Velocity  measurements  should  be  made  at  a  minimum  of  two  transects  by  following 
standard  streamflow  gauging  procedures.  Velocity  measurements  for  two  different 
flow  rates  is  desirable.  At  transects  where  velocities  are  not  measured,  the 
Manning's  equation  can  be  used  to  simulate  the  velocity  profiles  using  measured 
depth  profiles. 

Water  samples  must  be  collected  at  the  upstream  boundary,  at  the  effluent 
outfall  and  at  each  point  where  cross-sectional  depths  and  velocities  are 
measured.  However,  the  samples  can  be  collected  at  selected  points  at  each 
transect  (viz.,  less  points  outside  effluent  plume,  alternate  points,  etc.)  to 
reduce  sample  analyses  costs,  in  which  case,  the  concentrations  at  other  points 
are  obtained  by  interpolation. 
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In  order  to  account  for  fluctuations  in  effluent  water  quality  and  discharge  on 
the  instream  concentrations,  the  sampling  is  carried  out  either  by  following  the 
same  plug  of  water  beginning  at  the  outfall  and  proceeding  to  successive 
downstream  transects,  or  during  a  round-the-clock  intensive  survey  when  samples 
are  collected  at  each  point  at  3-4  hour  intervals.  Obviously,  the  selection  of 
a  sampling  methodology  depends  on  the  manpower,  time  and  other  resource 
constraints,  as  well  as  the  objectives  of  the  study. 

The  location  of  transects  can  be  based  on  preliminary  in-situ  measurements  of  a 
conservative  parameter  (eg.,  conductivity)  at  selected  access  points  to 
establish  the  approxiate  longitudinal  boundary  of  the  mixing  zone.  In-situ 
measurements  of  temperature,  pH  and  conductivity  must  be  taken  along  with 
collection  of  samples  which  are  to  be  analyzed  in  the  laboratory  for 
non-conservative  pollutants  of  concern.  In  some  cases,  it  may  be  desirable  to 
inject  a  solution  of  dye  continuously  to  gather  data  on  the  transverse 
distribution  characteristics  of  the  river.  This  is  particularly  useful  to 
simulate  effluent  discharge  from  proposed  outfall  locations  and  in  cases  where 
relocation  of  an  existing  outfall  is  being  considered.  The  dye  injection  must 
be  maintained  to  establish  steady-state  conditions  (about  2  or  3  hours).  The 
cross-sectional  distribution  of  dye  at  selected  transects  can  be  obtained 
directly  by  profiling  with  the  fluorometer. 

Generally,  two  surveys  should  be  carried  out  under  different  instream  hydraulic 
conditions  so  the  model  can  be  calibrated  with  one  survey  and  verified  with  the 
second  survey. 

Data  Analysis 

The  data  collected  during  one  of  the  field  surveys  is  utilized  to  determine  the 
parameters  required  for  modelling.  A  Fortran  computer  program  MIXANDAT  is 
utilized  to  perform  the  following  computations  using  the  survey  data  as  input: 

1.  Average  depth  and  velocity  at  each  transect. 

2.  Simulation  of  velocity  distributions  at  cross-sections  where  velocities  are 
not  measured. 

3.  Shape-velocity  factor  at  each  transect. 
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4.  Mass  flux  values  of  conservative  and  non-conservative  materials  at  each 
transect, 

5.  Variance  of  cross-sectional  distributions  of  conservative  materials  (used  to 
evaluate  the  dispersion  characteristics). 

The  Fortran  program  PREPARE  is  an  interactive  front-end  for  entering  the  survey 
data  to  be  analyzed  by  MIXANDAT. 

Mixing  zone  data  analysis  program  MIXANDAT  writes  the  analyzed  data  to  a  file 
called  OUTAN.DAT  which  can  be  printed.  An  example  of  OUTAN.DAT  is  in  Appendix 
D. 

The  plotting  program  PLTANDAT  gives  a  graphical  summary  of  the  field  data  as 
well  as  the  average  depth  and  average  velocity  at  each  transect  (Fig.  6  to  10). 

In  the  output  file  OUTAN.DAT,  the  variances  of  cross-section  distributions  are 
listed  for  each  transect. 

MIXANDAT  regresses  the  variances  against  different  parameters  to  determine  the 

dimensionless  dispersion  coefficient  6.  This  method  is  described  in  detail  in 

Reference  1.  The  program  MIXANDAT  chooses  the  regression  of  the  least  error  and 

equates  the  slope  of  the  linear  regression  with  an  expression  from  which  B  can 

be  determined.  This  value  of  B  is  returned  to  the  screen  for  use  in  the 
calibration  procedure. 

This  dimensionless  coefficient  B  is  varied  to  calibrate  the  predicted 
contaminant  spread  with  the  observed  concentration  distributions.  6  may  be 
considered  to  be  at  least  weakly  reach  dependent,  as  it  represents  the  rate  of 
spread  of  a  contaminant  plume  which  in  turn  depends  on  the  hydraulic  and 
geometric  properties  of  each  river  reach.  The  mathematical  basis  for  6  is 
discussed  in  more  detail  in  Appendix  A. 

Model  Calibration 

The  purpose  of  the  calibration  procedure  is  to  calculate  the  dispersion 
parameter  at  each  reach.  Program  MIXCALBN  calculates  a  cross-stream 
concentration  profile  at  each  transect  based  on  the  spread  Bi  at  transect  i. 
Program  COMPLOT  is  then  used  to  produce  a  plot  of  predicted  concentrations  and 
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observed  concentrations.  From  these  plots  the  user  can  adjust  the  value  of  6  at 
each  transect  to  match  the  predicted  with  observed  data.  The  user  then  runs 
MIXCALBN  again  with  new  estimates  of  B  at  each  reach  COMPLOT  is  run  after  each 
run  of  MIXCALBN.  When  the  distributions  of  contaminant  are  close  to  coinciding 
the  values  of  6i  are  recorded.  "!"hese  values  are  to  be  used  in  the  prediction 
models. 

Figures  11  to  13  show  an  example  of  the  sequence  required  to  obtain  a  good 
match. 

Applications 

Once  the  calibration  is  complete,  the  user  may  make  predictions  of  the  mixing 
zone  for  different  design  parameters.  Parameters  that  may  be  varied  include: 

1.  The  lateral  position  of  the  outfall. 

2.  Type  of  outfall  (pipe  or  diffuser). 

3.  Effluent  and  upstream  river  flow  rates. 

4.  Effluent  type  (different  decay  rates). 

5.  Effluent  concentration. 

6.  River  background  concentration. 

7.  River  pH  and  temperature. 

The  program  MIXPRED  makes  these  predictions  for  a  pipe  outfall  and  program 
MIXCADIF  does  so  for  a  diffuser  outfall.  The  program  CONMIX  plots  the  resulting 
concentration  distributions  in  the  receiving  water  in  three  different  formats 
(figures  14  to  16). 

These  programs  determine  the  effects  of  changing  certain  parameters  on  river 
distributions,  however  as  mentioned  in  the  Introduction,  the  PWQO  for  an 
effluent  discharge  are  based  on  the  concept  of  the  zone  of  passage.  Program 
MIXAPPLN  determines  the  "critical  points"  in  the  river  for  different  flow  rates, 
concentrations,  contaminants  and  management  options  such  as  outfall  location. 
"Critical  point"  are  locations  in  the  river  transect  where  PWQO  are  achieved. 
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MIXAPPLN  has  an  interactive  query/response  system  as  well  as  reading  files 
generated  by  MIXANDAT  and  PLTANDAT.  MIXAPPLN  handles  up  to  eight  upstream  flow 
rates,  six  effluent  flow  rates  and  six  temperature  values,  as  well  as  four  pH 
values  in  the  case  of  prediction  of  unionized  anmonia  levels.  The  output  is 
summarized  graphically  by  program  PLTCRIT  for  each  combination  of  design  options 
(see  Figure  17).  The  output  file  from  MIXAPPLN  (OUTAPP.DAT)  is  found  in 
Appendix  D. 

The  outputs  for  a  particular  case  in  MIXAPPLN  are: 

1.  Critical  values  of  concentrations  and  distances  from  the  outfall  at  .1,  .2, 
.3  and  .4  of  the  total  flow  of  the  river.  (That  is,  the  border  of  the  1st, 
2nd,  3rd  and  4th  stream  tube). 

2.  Mixing  Zone  Length  -  that  distance  from  the  outfall  to  the  point  downstream 
where  the  pollutant  is  fully  mixed  across  the  cross-section  and  the 
corresponding  average  zone. 

3.  Maximum  longitudinal   spread. 

A  more  detailed  description  of  the  critical  point  procedure  can  be  found  in 
Appendix  A. 
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DETAILED  USER'S  MANUAL 
Hardware  and  Software 

The  programs  in  this  model  package  are  written  in  FORTRAN-77  and  satisfy  all  the 
requirements  of  Microsoft  FTN-77  version  3.3.  The  mixing  zone  package  is 
currently  installed  on  a  COMPAQ  286  microcomputer  and  hence  will  run  on  IBM  PC 
and  any  lOOX  IBM  PC  compatibles  with  MS/DOS  operating  system.  256K  of  RAM  is 
required  to  run  the  largest  of  the  program's  executable  run  files.  The  plotting 
programs  use  the  PL0T88  graphics  package  which  is  licensed  copyrighted  software 
and  the  user  will  require  this  package  in  their  library.  The  user  will  receive 
the  mixing  zone  package  in  FORTRAN  source  code  and  must  compile  the  code  and 
load  the  library  modules  on  their  own  system.  The  compile/load  sequence  using 
Microsoft  FTN-77  is: 

FORI   program,,; 

PAS2 

LINK   /segments:256  program  ,,,  PL0T88  +  Altmath  +  Fortran; 

If  a  hard  disk  drive  is  present,  all  the  source  files  can  be  copied  from  the 
floppy  diskette  to  the  hard  disk  and  compiled  there.  If  there  is  no  hard  drive, 
all  the  executable  files  will  exceed  the  memory  of  a  single  diskette  and 
therefore  the  source  code  programs  should  be  copied  2  or  3  at  a  time  to  another 
floppy  and  compiled  there.  The  following  is  a  list  of  the  programs  in  the  order 
they  are  usually  used: 

PREPARE. FOR 

MIXANDAT.FOR 

PLTANDAT.FOR 

MIXCALBN.FOR 

COMPLOT.FOR 

MIXPRED.FOR 

MIXCADIF.FOR 

CONMIX.FOR 

MIXAPPLN.FOR 

PLTCRIT.FOR 
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MIXING  ZONE  PACKAGE  USER'S  MANUAL 

In  this  section  the  sequence  of  steps  involved  in  using  this  mixing  zone  package 
consists  of  a  text  description  with  pictures  of  the  corresponding  display  in  the 
computer  screen.  The  data  that  is  used  in  the  example  of  the  screens  is  just 
for  illustration  purposes  only  and  the  underlined  characters  are  those  that  are 
entered  by  the  user. 

1.  Key  PREPARE. 

This  program  will  respond  by  prompting  for  field  data,  specifically  velocity 
and  depth  data  for  at  least  one  transect.  It  is  most  favourable  to  have 
measurements  of  two  transects  and  at  least  flow  conditions  (i.e.  low  flow 
case  and  high  flow  case)  so  that  the  exponents  in  the  Leopold  -  Haddock 
equations,  as  well  as  the  total  river  flow  may  be  determined.  If  data  is 
available  at  only  one  flow  condition,  the  Leopold  -  Haddock  equations  will 
use  default  values  or  user-entered  values. 

The  following  diagrams  are  example  of  program  PREPARE  during  a  typical  run. 


DATA  UNITS 
DISTANCE  IN  METERS 
FLOWS  IN  CUBIC  METERS  PER  SECOND 
CONCENTRATIONS  IN  USER  CHOSEN.  CONSISTENT  UNITS 

ENTER  TITLE  OF  STUDY 

MISSISSIPPI  RIVER  MIXING  ZONE  STUDY 


HOW  MANY  FLOW  CONDITIONS  SURVEYED  (  1  OR  2  )  ?  :]_ 

Af  HOW  MANY  TRANSECTS  ARE  VELOCITY  MEASUREMENTS  AVAILABLE  ?  (1  OR  2) 
FOR  now  CONDITION  1  :  2 


11- 


/                     ENTER   THE 

rOLLOMlNG  DATA   FOR   TRANSECT      1 

U 

/               DEPTHS:     HOW 

MANY   MEASUREMENT   POINTS  ACROSS   THE   RIVER 

DEPTH      1 

0 

DEPTH      2 

1.67 

DEPTH      3 

8.69 

DEPTH      4 

7.78 

DEPTH      5 

6.11 

DEPIH     6 

5.56 

OfPTH      7 

4.44 

DEPTH     S 

3.0 

DEPIH     9 

2.?? 

DEPIH    10 

4.44 

/ 

\""-" 

1.11 

\ 

^ 

[^--      -      -    -                                                                                                                            ^, 

/                    DISTANCE    FROM 

SHORE  OF  EACH  OF   THE    11   POINTS                                            \ 

DISTANCE      1 

0. 

DISIANCE      2 

16.67 

DISTANCE      3 

100. 

DISTANCE      4 

168.89 

DISTANCE      5 

188.89 

DISTANCE     5 

205.56 

DISTANCE      7 

228.89 

DISTANCE      8 

300. 

DISTANCE      9 

341.11 

DISTANCE    10 

441.11 

DISTANCE    11 

L 

463.33 

\ 

/                                     ■              ^ 

ENTER  DEPTH  AVERAGED  VELOCITIES  AI  THE  11  POINIS 


VELOCITY 

1 

0. 

VELOCITY 

2 

.7 

VELOCITY 

3 

.9 

VELOCITY 

4 

VELOCITY 

5 

VELOCITY 

6 

VELOCITY 

7 

VELOCITY 

e 

VELOCITY 

9 

VELOCITY 

10 

.8 

VELOCITY 

11 

0. 
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ENTER   FIELD  DATA  COLLECTED  AT   RIVER   TRANSECTS 
ENTER   NUMBER  OF   TRANSECTS: _5_ 
NAME   OF   POLLUTANT:    AMMONIA 
FLOW  RATE  OF  EFFLUENT: 26 
CONCENTRATION  OF  EFFLUENT:    15. 
BACKGROUND  CONCENTRATION:  ^ 
BANK   OF   OUTFALL   (RIGHT  OR   LEFT):    RIGHT 
TEMPERATURE  OF  RIVER:   22.9 
DECAY  RATE  OF   CONTAMINANT    IN  RIVER:    .0000231 
TEMPERATURE   AT   WHICH  THIS   RATE    IS   KNOWN:    20^ 


TRANSECT  f   1 
DISTANCE   DOWNSTREAM   FROM  OUTFALL:    200 
HOW  MANY   MEASUREMENTS  ACROSS    THIS    TRANSECT:      U 


ENTER   DISTANCE   MEASURED  FROM  RIGHT  BANK 


0. 

16.67 

44.44 

63.33 

132.22 

150. 

177.78 

222.22 

333.33 

444.44 

551.11 
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Key  MIXANOAT. 

Program  MIXANDAT  performs  an  analysis  of  the  field  data, 
returns  a  value  of  BETA  to  be  used  in  subsequent  programs. 


This  program 


3.   Key  PLTANDAT. 

Program  PLTANDAT  gives  a  graphical  summary  of  the  field  data.  (See  Figures 
6  to  10  for  example).  PLTANDAT  is  responsible  for  writing  a  file  to  be  used 
by  Program  COMPLOT  and  hence  this  step  cannot  be  missed.  In  the  case  of 
plotting  to  the  screen,  strike  any  key  after  each  plot  to  continue  to  the 
next. 
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4.  Key  MIXCALBN. 

This  is  the  mixing  zone  calibration  program  and  is  used  to  adjust  the  model 
parameters  to  match  predicted  with  observed  concentration  data.  MIXCALBN  is 
typically  run  3  to  4  times  to  arrive  at  a  calibrated  situation.  In  order  to 
compare  the  observed  data  with  the  predicted  data  each  time,  run  the  program 
COMPLOT  after  each  run  of  MIXCALBN.  The  following  screen  diagrams 
illustrate  a  typical  first  run  through  MIXCALBN. 


HYDROOYNAHIC  DATA  ENTRY  AREA 


THE   EXPONENTS  FOR  THE   LEOPOLO-KADOOCK  EQNS  ARE 


HIOTH  EXP: 

.050 

DEPTH  EXP: 

.500 

VEL.     EXP: 

.450 

DO  YOU  WISH   TO  CHANCE   THEM?        YES-1      NO-0 
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OEC»r    OAff    DAT*   fHTRt    ADIA 

ENTER  A  DECAY   RATE   FOR   THE   RIVER   BACKGROUND: 

ENTER  A  DECA*   RATE   AT   EACH   TRANSECT 

TRANSECT  1  :  .0000231 

TRANSECT  2  :  .0000231 

TRANSECT  3  :  .0000231 

TRANSECT  4  :  .0000231 

TRANSECT  5  :  .0000231 


■0000231 


AT  WHAT   TEMPERATURE    IS  THIS  RATE  KNOWN?   IN  C   :   20 
WHAT    IS   THE   RIVER  TEMPERATURE:    IN  C   :   22.9 


DC   TOO  WISH   TO   CONSIDER  AMMONIA?      YES-l     NO-O        ] 
ENTER   PH  8.3 


•1 

/             ENTER  VALUES  Of  BETA 

] 

/             TRANSECT                    1    :    .00092 

TRANSECT                     2 

.00092 

IRANSECT                      3 

.00092 

IRAHSECT                     4 

.00092 

'RAHStCT                     b 

.00092 

IHE  OUTPUT  FILE  FROM  MIXCALBN. 

FOR    IS   CALLED 

•   CALOUT.DAT   • 

STOP   -   PRXRAM   TERMINATED. 

L 



^                                             ■                            \i 
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5.  Key  COMPLOT. 

Figures  11  to  13  show  an  example  of  the  graphics  output  by  COMPLOT.  At  each 
transect,  observed  concentrations  and  predicted  concentrations  are  graphed 
on  the  same  set  of  axes.  Values  for  mean  concentration  and  spread  for  both 
the  observed  and  predicted  concentrations  are  printed  as  well  as  the  value 
of  BETA  that  was  used  in  the  prediction.  It  is  suggested  that  calibration 
be  done  with  a  conservative  contaminant.  For  a  non-conservative 
contaminant,  its  decay  rate  must  be  obtained  from  literature.  The  value 
BETA  is  proportional  to  the  spread  of  the  predicted  distributions  and  hence 
by  increasing  BETA  at  a  transect  in  the  next  run  of  MIXCALBN,  will  increase 
the  spread  of  the  next  predicted  concentration  distribution.  The  mean 
concentrations  that  are  printed  correspond  to  the  concentration  that  would 
be  observed  if  the  same  amount  of  pollutant  were  fully  mixed  across  the 
river  transect.   It  should  be  equal  to: 

Ca  =  QEFL*CEFL 
QRS 

where  Ca  =  average  concentration 
QEFL  =  flow  rate  of  effluent 
CEFL  =  concentration  of  effluent 
QRS  =  total  flow  of  river  below  outfall 

The  predicted  curves  will  have  a  mean  concentration  Ca  at  all  transects 
(unless  the  spread  is  quite  small  and  then  it  may  overestimate  Ca) .  Hence 
deviations  from  this  mean  by  the  observed  data  reflect  the  quality  of  the 
data  and  it  is  not  necessarily  desirable  to  match  these  mean 
concentrations. 

A  second  estimate  of  BETA  is  made  for  each  transect  and  the  user  returns  to 
Step  4  (i-e.,  keys  MIXCALBN  and  enters  the  new  BETA  values).  For  a  true  far 
field  mixing  zone,  BETA  should  be  only  moderately  reach  dependent  and  lie  in 
a  range  of  .0025  to  .0001  approximately.  If,  however,  any  of  the  transects 
are  located  closer  to  the  outfall,  the  mixing  zone  area  may  fall  into  an 
intermediate  field  regime  where  the  BETA  values  are  dependent  on  distance 
from  the  source.  For  transects  near  the  outfall,  the  BETA  value  used  may 
have  to  be  larger  to  obtain  the  correct  spread.  MIXCALBN  is  run  until  the 
curves  and  spread  values  are  as  close  as  possible  then  the  calibration  is 
finished.  These  BETA  values  must  be  noted  for  use  in  subsequent  programs. 

Now  that  the  model  is  calibrated,  the  mixing  zone  package  may  be  used  for 
design  applications.  The  3  design  application  programs  are  MIXAPPLN, 
MIXPRED  and  MIXCADIF.  Any  of  these  three  programs  may  be  run  now. 

6.  Key  MIXAPPLN. 

MIXAPPLN  acts  interactively  to  input  different  combinations  of  river  flow, 
effluent  flow,  river  temperature  on  different  possible  pollutants.  An 
explanation  of  the  critical  point  method  is  given  in  Appendix  A.  The 
following  screen  diagrams  illustrate  how  to  run  MIXAPPLN. 
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REFERENCE  RIVER  PARAMETERS 


TOTAL  RIVER  FLOW  BELOW  OUTFALL  AT  TIME  OF  SURVEY    :       2374.020 

TRANSECT  DISTANCE  RIVER  WIDTH  AVERAGE  DEPTH  AVERAGE  VELOCITY 
3.97 
4.33 
2.19 
5.85 
3.40 

NOW   YOU  HAY   ENTER   DESIGN   PARAMETERS 
STRIKE   [ENTER]   TO  CONTINUE 


200.00 

551.11 

2000.00 

463.33 

4500.00 

677.78 

10150.00 

355.56 

17450.00 

750.00 

ENTER  STUDY   TITLE    :  MISSISSIPPI   RIVER  MUUG  ZONE  STUDY 


ENIER   POLLUTANT   NAME    :   AMMONIA 


r: y\ 

j                         ENTER   1  OF   UPSTREAM  FLOW  RATES 

9) 

.            \ 

ENTER   THESE   FLOW  RATES 

RATE     I:   3000 

RATE     2:  2700 

RATE     3:    2500 

RATE     4:    2000 

' 

RATE     5:    1500 

V 
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ENTER  #  OF  EFFLUENT  FLOW  RATES   (    7)    :    3 

ENTER   THESE   FLOW   RATES 
RATE     I:  _100 
RATE     2:   50_ 
RATE     3:    10 


ENTER  f  OF  RIVER  TEMPERATURES   (    7)    :J_ 

ENTER   THESE  TEMPERATURES 

TfMP.  1:   25. 

TEMP.  2:    22. 

TEMP.  3:    le 

TEMP.  i:    10 


ARE    you  CONSIDERING  AtflONIA?   UYES     0-NO   : 

ENTER  I  OF  PH  VALUES   (    5)    : _2_ 

ENTER   THESE   PH  VALUES 
PH   1:      7. 
PH  2:     8JI_ 

ENTER  EFFLUENT  CONCENTRATION   :    15 

ENTER  BACKGROUND  CONCENTRATION   :   0. 

ENTER  PROVINCIAL  WATER  QUALITY  OBJECTIVE   : 

ENfER   TEMPERATURE   COEFFICIENT    :    1.106 

ENTER  DECAY  RATE  OF  BACKGROUND   :    .0000231 


ENTER   THE   TEMPERATURE  THIS  RATE    IS  KNOWN  AT   :   20 
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r-" —             — \ 

/         ENTER   VALUES  OF   BETA                                                                                                           \ 

/                      TRANSECT         1:    .004                                                                                            \ 

/                       TRANSECT         2:    .0017 

TRANSECT         3:    .00078 

TRANSECT         4:    .001 

TRANSECT         5:    .0007 

ENTER  DECAY  RATES  AT  EACH  TRANSECT 

TRANSECT         1:    .0000231 

TRANSECT         2:    .0000231 

TRANSECT         3:    .0000231 

TRANSECT         4:    .0000231 

TRANSECT         5:    .0000231 

7.  Key  PLTCRIT. 

PLTCRIT  outputs  a  graphical  summary  of  the  output  from  MIXAPPLN  for  each  of 
the  design  option  combinations  (see  Figure  17). 

8.  Key  MIXPRED. 

MIXPRED  computes  the  two  dimensional  concentration  distribution  for  a  pipe 
outfall  located  some  distance  from  shore  for  user  defined  design  parameters. 
The  following  screen  diagrams  illustrate  how  to  use  MIXPRED. 


/                   THE   EXPONENTS  FOR  THE   UOPOLO-MADOOCIC 

WIDTH  EXP:                  .050 
DEPTH  EXP:                 .500 
VEL.     EXP:                  .450 

DO  rOU  WISH  TO  CHANGE   THEM?     rES-l     NO-O 

EQNS  ARE                                        \ 
:  0 

■20- 


RIVER   FLOW   RATE   ABOVE   OUTFALL 
OOTFALL   FLOW  RATE   -   76.00 


2298.02 


00   YOU  WISH   TO  CHANGE   EITHER?      YES-l     NO-0 
RIVER   FLOR   RATE    •   3000 
OUTFALL   FLOW  RATE  TloO 


EFFLUENT  CONCENTRATION   IS    :    16.00 
DO  YOU  WISH  TO  CHANGE    IT?     YES-l  NO-0 
THE   NEW  EFFLUENT  CONC.   •  20 


THE  BACKGROUND  CONCENTRATION   IS    :    .00 
00  YOU  WISH  TO  CHANGE    IT?     YES=1     NO  =0 
ENTER  NEW  BACKGROUND  CONC:      2. 
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ENICR  A  OECAY  RATE   FOR  THE   RIVER  SACKGROUNO:      . 0000231 
AT  yHAT  TEMPERATURE    IS  THIS  RATE  UNOHN?      IN  C   :   20. 
WHAT    IS   THE  RIVER   TEMPERATURE?    IN   C    :    22.9 


THE   OUTFALL    IS  AT  SHORE 
00  YOU  yiSH  TO  CHANGE    IT?     YES-1     NO-0  1 

ENTER  THE  DISTANCE  OF   THE  OUTFALL  FROM  THE  BANK   :    100. 


00  YOU  WISH  TO  COMSIOCR  AMMONIA?     VES-1     N0>0  l 

ENTER  PH  8.3 
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ENTER   5  VALUES  OF  DECAY 


.0000231 
.0000231 
.0000231 
.0000231 
.0000231 


ENTER   VALUES  OF   BETA 


TRANSECT  1 

TRANSECT  2 

TRANSECT  3 

TRANSECT  A 

TRANSECT  5 


.004 

.0017 

.00078 

.001 

.0007 


THE   OUTPUT   FILE   FROH  MIXPRED.    FOR    IS   CALLED 
"   PREDOUT.DAT   " 
STOP   -   PROGRAM   TERMINATED. 


9.  Key  CONMIX. 

CONMIX  gives  the  choice  of  3  different  graphical  representations,  of  which 
the  user  may  choose  1,  2  or  all  of  them  (Figure  14  to  16). 

10.  Key  MIXCADIF 

.Program  MIXCADIF  is  almost  identical  to  program  MIXPRED  except  that  MIXCADIF 
applies  to  predictions  of  distributions  when  diffuser  outfalls  are  used. 

11.  Key  CONMIX 

CONMIX  also  gives  graphical  summary  of  results  from  MIXCADIF. 


PlOf  THE  FOLLOWING  TYPE  OF  PLOT   UYES  0=N0 
CROSS  STREAM  CONCENTRATION  PROFILES:   1 

2  DIMENSIONAL  CONTOUR  PLOT:   1 

3  DIMENSIONAL  MESH  PLOT:   1 

ENTER   POLLUTANT   NAME    (30   LETTERS  MAX)    :   AMMONIA 
PLOT   lOCATION    1-SCREEN   2-PLOTTER   3-PRINTER    :    1 
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FIGURES 


DATA  COLLECTION 
'PREPARE' 


DATA  ANALYSIS 
'MIXANDAT' 


DATA  ANALYSIS  PLOTS 
'PLTANDA? 


i 

i 

CALIBRATION 

'MIXCALBN' 

EMTER  NEW  VALUES 

i 

OF  BETA 

CALIBRATION  PI  OTS 

'COMPLOT 

i 

I 

1 

' 

NO 

- 

PREDICTION  . 
MODELS 

1 

r 

1 

\ 

r 

DIFFUSER  OUTFALL 
'MIXCADIP 

PIPE  OUTFALL 
'MIXPRED' 

CRITICAL  POINTS 
'MIXAPPLN' 

\ 

I                                  1                                 1 

CONCENTRATION  PLOT 
'CONMIX' 

CONCENTRATION  PLOT 
'CONMIX' 

CRmCAL  POINT  PLOT 
'PLTCR[r 

Figure  5   Flow  Diagram  of  Mixing  Zone  Package 
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lo)      CONCENTRATION     PROFILES  (b)    LONGITUDINAL    PROFILE    ALONG    A    LATERAL    BOUNDARY 

AT  CROSS-SECTIONS  O  TO©  OF    LIMITED    USE    ZONE. 


Figure  18   Critical  Point  of  LUZ 
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APPENDIX  A 
BASIC  STREAM  TUBE  MODEL ^ 


The  fundamental  concept  of  the  stream  tube  model,  developed  by  Yotsukura  and 
Cobb,  (1972)  is  the  use  of  the  cumulative  partial  discharge,  q,  at  a  given 
cross-section  instead  of  the  lateral  distance,  y,  as  the  independent  variable. 
In  this  approach,  the  river  cross  section  is  divided  into  a  number  of  vertical 
strips  termed  "stream  tubes",  such  that  the  discharge  within  each  stream  tube  is 
the  same.  Thus,  the  cross  sectional  concentration  distributions,  c(x,q), 
predicted  by  the  stream  tube  model  will  be  functions  of  q.  These  distributions 
can  be  transformed  into  c(x,y)  as  a  function  of  distance  from  the  bank,  y,  by 
knowing  the  relation  between  q  and  y  at  each  transect. 

The  derivations  of  the  basic  equations  of  the  stream  tube  model  have  been 
presented  by  Yotsukura  and  Cobb  (ibid);  they  are  subject  to  the  following 
assumptions: 

1.  The  density  of  effluent  (or  solute)  is  the  same  as  that  of  the  receiving 
water.  This  assumption  is  reasonably  satisfactory  for  most  of  the  municipal 
effluent  discharges  to  rivers. 

2.  The  concentration  distributions  in  the  far  field  are  not  affected  by  the 
near  field  mixing  processes  (eg.,  dilution  due  to  initial  momentum  of  jet). 
Usually,  the  jet-induced  diffusion  approaches  the  ambient  diffusion  for  a 
short  distance  below  a  source  in  a  shallow  river. 

3.  The  depth  distribution  of  effluent  (or  solute)  in  the  river  channel  is 
uniform.  Generally,  the  longitudinal  distance  required  to  attain  depth 
uniformity  is  short  in  shallow  rivers,  being  of  the  order  of  50  to  100  times 
the  channel  depth;  thus,  the  assumption  is  justified. 

4.  The  transport  due  to  longitudinal  dispersion  is  negligible.  In  the  case  of 
continuous  effluent  discharge,  this  transport  is  very  small  in  relation  to 
that  due  to  convection  and  lateral  dispersion,  thus  justifying  the 
assumption. 

The  Convective-Dispersion  Equation 

The  two-dimensional  convective-dispersion  equation  for  a  non-conservative 
material  in  the  far  field  region  of  the  mixing  zone  can  be  written  in  the 
following  form: 


From  Reference  1. 


3c  =  Dy  9^C  -  ^x^^^ 


(2) 


ax 


aq^ 


■ 
■ 


where  x  =  distance  downstream  from  outfall 
q  =  partial  cumulative  discharge 

Dy  =  lateral  diffusion  factor 

Kd  =  decay  coefficient 
u  =  depth-averaged  local  velocity 
c  =  c(x,q)  the  2-D  concentration  field 


A  solution  of  Equation  2  is: 

c(*.p)  =  R'Ca  (4114.)  "^/^ 


exp  )-  (2n  -H  ps  -*■  p) 
I      4^ 


j^.[ 


exp  )-  (2n  +   PS  -  p)  ^ 
4* 


I   I  exp  )-  (2n  -  PS  -  p)'' 


+  exp  )-  (2n__ps  +  p) 


(3) 


where 


=  Dyx 


p  =  q/Q,   Ps  =  Qs/Q.   Cj 


CeQe 
Q 


where  n  =  number  of  images  required  to  account  for  the  effect  on 
concentration  of  reflection  from  channel  banks 
Ce  =  effluent  concentration 
Qe  =  effluent  flow  rate 
Q  =  discharge  of  river  below  outfall 
R'  =  exp  (-Kdx/U)  decay  factor 


Diffusion  Factor 


The  dimensionless  diffusion  factor,  4),  in  Equation  3  is  a  time  dependent 
parameter  but  since  we  are  considering  a  steady-state  situation,  it  is  directly 
dependent  on  the  distance  downstream  from  the  outfall  and  can  be  written  as 


4> 


B  X 
b 


where  b  =  surface  width  at  a  transect 

The  parameter  B  (BETA)  is  an  input  parameter  for  the  mixing  zone  model  and  is 
determined  at  first  by  the  program  MIXANDAT.  During  the  calibration  procedure  B 
may  be  changed  to  enable  predicted  concentration  distribution  to  agree  with 
observed  data.  A  more  thorough  description  of  this  parameter  as  well  as  the 
mathematical  basis  for  the  model  may  be  found  in  Reference  2.  The  Typical 
values  of  B  for  a  shallow  river's  far  field  mixing  zone  are  in  the  range  .0001  - 
.002,  however  for  transects  that  are  relatively  close  to  the  outfall,  6  may  be 
artificially  large  and  may  be  dependent  on  the  distance,  x,  from  the  outfall. 

Critical  Point  Analysis 

The  boundary  of  a  limited  use  zone  (LUZ)  is  generally  identified  by  lateral  and 
longitudinal  co-ordinates  with  respect  to  the  outfall.  Usually,  the  lateral 
boundary  of  a  LUZ  is  limited  to  the  range  0.2  to  0.4  times  the  river  discharge 
Q.  The  cumulative  partial  discharge  between  the  outfall  bank  and  the  accepted 
lateral  boundary  will  be  denoted  by  the  non-dimensional  parameter,  Pl  =  q|_/Q.  A 
specified  pollutant  concentration  criterion,  Cj  must  be  met  within  qL.  For  a 
given  set  of  values  of  Ce,  Qe,  qL  and  Q,  the  longitudinal  distribution,  c(xi, 
qL)  predicted  by  Equation  3  attains  a  maximum  value  at  some  x-j  and  then  follows 
a  decreasing  trend  as  shown  in  Figure  IS.  The  point  at  which  the  concentration 
attains  the  maximum  value  is  termed  the  "critical  point".  The  longitudinal 
distance  between  the  outfall  and  the  critical  point  is  termed  the  "critical 
distance",  Xi_,  and  the  maximum  concentration  is  termed  the  "critical 
concentration",  Cl-  Detailed  description  of  the  method  of  computation  of  X^  and 
Cl  are  presented  by  Gowda  (1980) (Reference  1). 

Knowing  Cl  and  Cj,  the  allowable  effluent  concentration,  CgAi  can  be  calculated 
from  the  following  expression: 

CeA  =  ^e_^ 


The  maximum  longitudinal  boundary  of  LUZ,  Xj,  occurs  along  the  discharge 
shoreline  in  the  case  of  a  bank  outfall  as  shown  in  Figure  1. 
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APPENDIX  B 
RELATIONSHIP  OF  DEPTH,  WIDTH 
AND  VELOCITY  TO  STREAMFLOW 


The  streamflow  Q,  at  any  stream  cross-section  is  directly  related  to  the 
cross-sectional  area  A  and  the  mean  velocity  by 

Q  =  AU  C-1 

The  area  of  a  cross-section  is  equal  to  the  product  of  mean  depth  H,  the  mean 
top  width  W,  and  hence 

Q  =  HUW  C-2 

When  the  streamflow  Q  is  changed,  the  width,  depth  and  velocity  will  be 
affected.  The  following  general  relationships,  as  developed  by  Leopold  and 
Haddock  may  be  used  to  derrive  the  new  H,  U  and  W 

H  =  aQ"^  C-3 

U  =  bQ^  C-4 

W  =  cQ9  C-5 

Where  a,  b  and  c  are  empirical  constants  and  d,  f  and  g  are  exponents  which  are 
functions  of  the  hydraulic  radius,  slope  and  roughness  of  the  channel.  It  can 
be  shown  that  by  substitution  of  equations  C-3,  C-4  and  C-5  into  C-2,  that 

abc  =  1  C-6 

d  +  f  +  g  =  1     .  C-7 

These  equations  provide  us  with  a  hydraulic  model  requiring  only  easily 
determined  empirical  exponents  (d,  f  and  g).  They  are  extremely  useful  to 
estimate  the  width,  depth  and  velocity  of  a  reach  under  flow  conditions  not 
surveyed. 

The  Manning  equation  provides  another  hydraulic  model  that  requires  a  knowledge 
of  3  of  4  parameters;  these  are  velocity  V,  channel  slope  S,  hydraulic  radius  R 
and  roughness  (n)  of  the  channel 


=  r2/3  sl/2 


The  Manning  equation  is  accepted  as  more  accurate  than  the  corresponding 
Leopold-Haddock  equation  (C-4);  in  water  quality  engineering,  however,  the 
latter  is  more  cotimonly  used  because  of  its  simplicity. 

The  mixing  zone  model  presented  in  this  manual  computes  the  Leopold-Haddock 
exponents  if  surveys  are  done  at  different  flow  conditions  by  program  FLOWCAL. 

In  the  absence  of  data  on  a  stream  under  study,  Leopold-Haddock  coefficients  can 
be  derrived  by  referring  to  the  Hanning  equation.  For  a  rectangular  channel, 
the  Hanning  equation  gives  f  =  0.4,  d  =  0.6  and  g  =  0.  In  the  mixing  zone 
model,  the  default  values  of  the  exponents  are  g  =  0.05,  f  =  0.45  and  d  =  0.50. 
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FORTRAN  SOURCE  CODE 


1. 

PREPARE. FOR 

2. 

MIXANDAT.FOR 

3. 

PLTANDAT.FOR 

4. 

MIXCALBN.FOR 

5. 

COMPLOT.FOR 

6. 

MIXAPPLN.FOR 

7. 

PLTCRIT.FOR 

8. 

MIXPRED.FOR 

9. 

MIXCADIF.FOR 

10. 

CONMIX.FOR 

APPENDIX  C 
Fortran  Source  Code 
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SDEBUG 

(^***t«/>««*/ttttAjm/nm****««  ******  ***«****t*/n>«*****/>*tt*/m/mtt**t*** 

C  ♦ 

C  * 

C  PROGRAM  PREPARE   READS  ALL  THE    INPUT  DATA  AND  MAKES  A  FILE  * 

C  TO  BE  READ  BY  THE  REST  OF  THE   PROGRAMS.  * 

C  PREPARE  ALSO  CALCULATES  THE  TOTAL  FLOW  FROM  VELOCITY  * 

C  MEASUREMENTS  AND  CALCULATES  THE   LEOPOLD-MADDOCK  EXPONENTS  * 

C  IF  TWO  FLOW  SITUATIONS  ARE  MEASURED.  * 

C  WRITTEN  BY   R.    JARVIS  * 

C  GORE  &  STORRIE      1986  * 

Q**-k**-t**1rk-kit-lrk1rk1rkirk1rk**1rk1rk1rk1rk1rk1rkirir1rk1rir1rk1rki^ 

DIMENSION  X(10),D(10.30),Y(10,30),C(10,30).ICR(10),V(10.30) 
*.MVEL(10) 
CHARACTER*80  TITLE 
CHARACTER*20  CPARM 
CHARACTER*5  OUTBNK 

0PEN(5. FILE= 'FIELD. DAT '.STATUS=' NEW') 
CALL  CLRSCN 
WRITE(*.5558) 
5558     FORMAT(/'  DATA  PREPARATION  PROGRAM'/ 

1'  7/ 

2'  DATA  UNITS'/ 

3'  DISTANCE    IN  METERS'/ 

4'  FLOWS   IN  CUBIC  METERS  PER  SECOND'/ 

5'  CONCENTRATIONS    IN  USER  CHOSEN,   CONSISTENT  UNITS') 

WRITE(*,1) 
1  FORMAT(//'  ENTER  TITLE   OF   STUDY'//'  '\) 


READ(*,2)TITLE 
2  FORMAT(A) 


C 


CALL   FLOWCAL(QRS) 


CALL   CLRSCN 
WRITE(*,3) 

3  FORMAT(//'  ENTER   FIELD  DATA  COLLECTED  AT  RIVER  TRANSECTS') 
WRITE(*,4) 

4  FORMAT(/'        ENTER  NUMBER  OF   TRANSECTS:    '\) 
READ(*,*)NTR 

WRITE(*,5) 

5  FORMAT(/'        NAME  OF  POLLUTANT:    '\) 
READ(*,2)CPARM 

WRITE(*,7) 

7  FORMAT(/'        FLOW  RATE  OF  EFFLUENT:    '\) 
READ(*,*)QEFL 

QRUP=QRS-QEFL 
WRITE(*,8) 

8  FORMAT(/'        CONCENTRATION  OF  EFFLUENT:    '\) 
READ(*,*)CEFL 

WRITE(*,9) 

9  FORMAT(/'        BACKGROUND  CONCENTRATION:    '\) 
READ(*,*)CBKG 

WRITE(*.10) 

10  FORMAT(/'        BANK  OF  OUTFALL    (RIGHT  OR   LEFT):    '\) 
READ(*,2)0UTBNK 

IF(OUTBNK.EQ. 'R '.OR. OUTBNK.  EQ.'Rr)OUTBNK=' RIGHT' 
IF(OUTBNK.EQ. 'RIG'. OR. OUTBNK. EQ.'r')OUTBNK='RIGHT' 
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IF(OUTBNK.EQ.'RIGH'.OR.OUTBNK.EQ.'RIGHT')OUTBNK='RIGHT' 
IF(OUTBNK.EQ. 'ri ' .OR.OUTBNK.EQ. 'rig ' )OUTBNK=' RIGHT' 
IF(OUTBNK.EQ. 'right ' .OR.OUTBNK.EQ. ' righ ' )OUTBNK=' RIGHT' 
IF(OUTBNK.EQ.'L'.OR.OUTBNK.EQ.'LE')OUTBNK='LEFT' 
IF(OUTBNK.EQ. '1'. OR.OUTBNK.EQ. 'le')OUTBNK=' LEFT' 
IF(OUTBNK.EQ. ' LEF ' .OR.OUTBNK.EQ. ' LEFT' )OUTBNK=' LEFT' 
IF(OUTBNK.EQ. ' lef .OR.OUTBNK.EQ. ' left ' )OUTBNK=' LEFT ' 
IF(OUTBNK.NE. 'RIGHT'. AND. OUTBNK.NE. 'LEFT' )THEN 
WRITE(*,11) 

11  FORMAT(/'        TEMPERATURE  OF  RIVER:    '\) 
READ(*,*)TEMP 

WRITE(*,12) 

12  FORMAT(/'    DECAY  RATE  OF  CONTAMINANT  IN  RIVER:  '\) 
READ(*,*)RKS 

WRITE(*,13) 

13  FORMAT(/'    TEMPERATURE  AT  WHICH  THIS  RATE  IS  KNOWN:  '\) 
READ(*,*)TMP 

C 

DO  100  1=1, NTR 
CALL  CLRSCN 
WRITE(*,14)I 

14  F0RMAT(//'    TRANSECT  #  ',12/) 
WRITE(*,15) 

15  F0RMAT(/'    DISTANCE  DOWNSTREAM  FROM  OUTFALL:  '\) 
READ(*,*)X(I) 

WRITE(*,300) 
300   FORMAT(/'    HOW  MANY  MEASUREMENTS  ACROSS  THIS  TRANSECT:  '\) 
READ(*,*)ICR(I) 
CALL  CLRSCN 
WRITE(*,17)0UTBNK 

17  FORMAT(/'    ENTER  DISTANCE  MEASURED  FROM  ' ,A, '  BANK  '/) 
DO  500  J=1,ICR(I) 

18  WRITE(*,19)J 

19  F0RMAT{2X,I2,':  '\) 
READ(*.*)Y(I,J) 

500  CONTINUE 
CALL  CLRSCN 
WRITE(*,21)0UTBNK 

21  FORMAT(/'   ENTER  DEPTH  FROM  ',A,'  BANK   '/) 
DO  501  J=1,ICR(I) 

22  WRITE(*,23)J 

23  F0RMAT(2X,I2,':  '\) 
READ(*.*)D(I,J) 

501  CONTINUE 
CALL  CLRSCN 

104   WRITE(*,25)0UTBNK 

25  FORMAT(/'   ENTER  CONCENTRATIONS  FROM  ',A,'  BANK   '/) 
DO  502  J=1,ICR(I) 

26  WRITE(*,27)J 

27  F0RMAT(2X,I2,':  '\) 
READ(*.*)C(I,J) 

502  CONTINUE 
MVEL(I)=0 

100   CONTINUE 

C 

C     NOW  WRITE  THE  DATA  INTO  THE  PROPER  FORMAT  INPUT  FILE 

C 
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WRITE(5,2)TITLE 
FEMT=1. 
Fl=l. 
F2=.67 

WRITE(5,200)NTR,QEFL,FEMT,F1,F2 
200   F0RMAT(2X,I2,2X,F8.3,2X,F6.2.2X,F6.2,2X,F6.2) 
WRITE(5,2)0UTBNK 
DO  201  1=1, NTR 
WRITE(5, 202)1 

202  FORMATC   TRANSECT  ',12) 
WRITE(5,203) 

203  FORMAT( 'BLANK') 
QRIVER=QRS 
WRITE(5,205)X(I),ICR(I),QRIVER,MVEL(I) 

205  FORMAT(2X,F10.2,2X,I3,2X,F10.3,2X,I2) 
WRITE(5,206)(Y(I,J),J=1,ICR(I)) 
WRITE(5,206)(D(I,J).J=1,ICR(I)) 
C1=0. 

WRITE(5,900)CPARM 
900        FORMAT(A) 

WRITE(5,*)CBKG,CEFL,C1,C1.C1,C1 
WRITE(5,206)(C(I,J),J=1,ICR(I)) 

206  F0RMAT(2X,25(1X,F9.2)) 
WRITE(5,207) 

207  FORMAT( 'NOCONC    ') 
WRITE(5,208) 

208  FORMATC  -999.0  0.0  0.0  0.0  0.0  0.0  0.0  ') 
201   CONTINUE 

CALL  CLRSCN 

WRITE(5,5559)CHAR(26) 
5559  FORMAT(A) 

CL0SE(5) 

STOP 

END 
C 
C 
C 
C 

SUBROUTINE  FLOWCAL(QRSl) 

COMMON/LM/BEX,HEX,UEX 

COMMON/IN/  D(30).Y(30).V(30),DAVE(2),VAVE(2) 

COMMON/QQ/  QRS(2),QAVE(2) 

COMMON/YM/  YMAX(2) 

CALL  CLRSCN 

WRITE(*,1) 

1  FORMAT(//'  FLOW  CALCULATION  PROGRAM') 
WRITE(*,2) 

2  FORMATC  ') 

WRITE(*,3) 

3  FORMAT(//'   HOW  MANY  FLOW  CONDITIONS  SURVEYED   (  1  OR  2  )  ?  :'\) 
READ(*,*)IFL 

WRITE(*,4) 

4  FORMAT(//'   AT  HOW  MANY  TRANSECTS  ARE  VELOCITY  MEASUREMENTS  AVAILA 
♦BLE  ?  (1  OR  2)  :  ') 

WRITE(*,5) 

5  FORMAT(/'   FOR  FLOW  CONDITION  1  :  '\) 
READ(*,*)NTR1 
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10 


12 


110 


16 


11 


13 


14 


IF(IFL.EQ.2)THEN 

WRITE(*,6) 

FORMAT(/'    FOR  FLOW  CONDITION  2  :  '\) 

READ(*,*)NTR2 
ENDIF 

CALL  FL0(1,NTR1,QRS1) 
IF(IFL.EQ.2)CALL  FL0(2,NTR2,QRS2) 
IF(IFL.EQ.2)CALL  LE0MAD(QRS1,QRS2) 
RETURN 
END 


SUBROUTINE  FLO( JFL.NTR.QRSA) 
COMMON/ IN/  D(30),Y(30).V(30),DAVE(2),VAVE(2) 
COMMON/QQ/  QRS(2),QAVE(2) 
COMMON/YM/  YMAX(2) 
DIMENSION  QQRS(2) 
DSUM=0. 
VSUM=0. 
11  =  1 

WRITE(*,10)JFL 
FORMAT(/' 
DO  99  11=1, NTR 
DSUM=0. 
VSUM=0. 
CALL  CLRSCN 
WRITE(*,1)II 
FORMAT(/' 
WRITE(*,2) 
FORMATC 
WRITE(*,3) 
FORMAT(//' 
■^R  ?'\) 

READ(*,*)IDEEP 
WRITE(*,12) 
FORMAT(/) 
DO  5  I=1,IDEEP 
WRITE(*,6)I 
FORMATC  DEPTH 
READ(*,*)D(I) 
CONTINUE 
CALL  CLRSCN 
WRITE(*,110)IDEEP 

FORMAT(//'   DISTANCE  FROM  SHORE  OF  EACH  OF  THE 
DO  11  I=1,IDEEP 
WRITE(*,16)I 
FORMATC  DISTANCE 
READ(*.*)Y(I) 
CONTINUE 
CALL  CLRSCN 
WRITE(*,13)IDEEP 
FORMAT(//'   ENTER  DEPTH  AVERAGED  VELOCITIES  AT  THE  '.12,'  POINTS' 


FLOW  CONDITION  ',12) 


ENTER  THE  FOLLOWING  DATA  FOR  TRANSECT  '  ,12) 

') 

DEPTHS:   HOW  MANY  MEASUREMENT  POINTS  ACROSS  THE  RIVE 


.12,' 


'\) 


.12.'  POINTS'//) 


12  ' 


'\) 


7/) 

DO  15  I=1,IDEEP 
WRITE(*,14)I 
FORMATC   VELOCITY 
READ(*,*)V(I) 


12 


'\) 
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15    CONTINUE 

JDEEP=IDEEP-1 

DO  60  I=1,J0EEP 

DSUM=DSUM+(Y(I+l)-Y(I))*(D(I+l)+0(I))/2. 

VSUM=VSUM+(Y(I+l)-Y(I))*(V(I+l)+V(I))/2. 
60    CONTINUE 

QQRS(II)=DSUM*VSUM/Y(IDEEP) 

IF(II.EQ.1)THEN 

DAVE(JFL)=DSUM/Y(IDEEP) 
VAVE(JFL)=VSUM/Y(IDEEP) 
YMAX(JFL)=Y(IDEEP) 

ENDIF 

CALL  CLRSCN 

WRITE(*,101)JFL 
101   FORMAT(/'    FLOW  CONDITION  :  ',12) 

WRITE(*,100)QQRS(II) 
100   F0RMAT(/'    TOTAL  FLOW  IN  RIVER  =',F10.2) 
99    CONTINUE 

IF(NTR.EQ.2)QRSA=(Q0RS(l)^QRS(2))/2. 

IF(NTR.EQ.1)QRSA=QQRS(1) 

WRITE(*,1010)JFL,QRSA 
1010  FORMAT(//'   AVERAGE  FLOW  RATE  FOR  RIVER 
CONDITION  ',I2,'=',F10. 
*3) 

RETURN 

END 
C 
C 

SUBROUTINE  LE0MAD(QRS1,QRS2) 

COMMON/LM/  BEX,HEX,UEX 

COMMON/ IN/  D(30),Y(30),V(30),DAVE(2),VAVE(2) 

COMMON/YM/  YMAX(2) 

COMMON/QQ/  QRS(2),QAVE(2) 

QRAT=L0G(QRS1/QRS2) 

DRAT=L0G(DAVE(1)/DAVE(2)) 

VRAT=L0G(VAVE(1)/VAVE(2)) 

WRAT=L0G(YMAX(1)/YMAX(2)) 

BEX=WRAT/QRAT 

HEX=DRAT/QRAT 

UEX=VRAT/QRAT 

WRITE(*,1) 

1  FORMAT(/'  LEOPOLD-MADDOCK  EXPONENTS') 
WRITE(*,2)BEX 

2  FORMAT(/'  WIDTH  EXPONENT 
',F5.3) 

WRITE(*,3)HEX 

3  FORMAT(/'    DEPTH  EXPONENT 
'.F5.3) 

WRITE(*.4)UEX 

4  FORMAT(/'  VELOCITY  EXPONENT 
'.F5.3) 

RETURN 

END 
C 
C 
C 

SUBROUTINE  CLRSCN 
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WRITE(*,101)CHAR(27),'[2J' 
101   FORMAT (IX, A, A\) 
RETURN 
END 
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f*  *-*■***-*•♦♦**  *■ 

♦* 

c  * 

C  PROGRAM  MIXANDAT. FOR  * 

C  DATA  ANALYSIS  PROGRAM  FOR  MIXING  ZONE   PACKAGE                                          * 

C  WRITTEN  BY  T.    P.   H.   GOWDA  ♦ 

C  PREPARED  FOR  USE  ON  PC  BY  R.   JARVIS                                                               * 

C  GORE  &  STORRIE      1986  ♦ 

C  * 

C  * 

COMMON/A/  X(8).Y(8,60).Z(8.60),VEL(8.60),C0NC(8.60) 
COHMON/B/  YL(8,6O),ZL(8,60).C0NL(8.6O),UVL(8.6O).VUF(5).VCN(5) 
♦.SUMA(60),DVEL(60).DELQ(60).SUMQ(60).DC0NC(60),SBF(5) 
*.FLUX(60),DELA(60).SUMF(60),SHF(5) 
*.Q(8),U(8).YW(60).UNIF(60).VCMX(5),SBS(5).SHS(5) 
*,CMX(5),VCQ(5),SBQ(5),SHQ(5).SQQ(5),VPQ(5),ZZAV(25) 
CHARACTER*80  TITLE 
CHARACTER*20  CPARM.TRNSCT 
CHARACTER*10  OUTBNK.REFLD 
C 

C  DATA   INPUT  *** 

C  DATA   IS  READ  FROM  FILE    'FIELD.DAT'    WHICH   IS  FROM  PROGRAM    'PREPARE' 

C 

C  'OUTAN.DAT'    IS   THE   FULL  OUTPUT  FILE 

C 

C  'PINCAL.DAT'    IS  THE    INPUT   FILE   TO  MIXCALBN.FOR 

C 

C  'BETA. DAT'    IS  THE   FILE   USED   IN  SUBROUTINE   REGRESS  TO  CALCULATE 

C  THE   OVERALL  BETA  VALUE   FOR  THE   RIVER 

c 

C  'PLOTAN.DAT'    IS   THE    INPUT  FILE   FOR   PROGRAM  PLTANDAT.FOR 

C 

0PEN(9, FI LE=' BETA. DAT ',STATUS=' NEW') 

0PEN(8,FILE='PL0TAN.DAT',STATUS='NEW') 

0PEN(7,FILE='PINCAL.DAT',STATUS='NEW') 

0PEN(5,FILE='FIELD.DAT',STATUS='0LD') 

0PEN(6,FILE='0UTAN.DAT',STATUS='NEW') 

IP-1 

READ(5,10)TITLE 

WRITE(6,10)TITLE 

WRITE(8,10)TITLE 
10         FORMAT(A) 

READ(5,*)NTR,QEFL,FMET,F1.F2 

WRITE(9.2299)NTR 
2299     F0RMAT(2X,I3) 

READ(5.10)0UTBNK 

WRITE(8.1331)NTR 
1331     F0RMAT(3X,I2) 

IF(FMET.LE.O.O)    FMET=1.0 

IF(Fl.LE.O.O)    Fl=1.0 

IF(F2.LE.0.0)    F2=0.67 

QEFL=QEFL*FMET**3 

DO   100   1=1, NTR 

READ(5,10)TRNSCT 

READ(5,10)REFLD 

READ(5,*)X(I),NYZ.Q(I),MVEL 
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READ(5,*)(YL(I.J),J=1.NYZ) 

READ(5,*)(ZL(I.J),J=1,NYZ) 

IF(REFLD.NE. 'BLANK')  GO  TO  38 

IF(OUTBNK.EQ.'LEFT')  REFLD='LEFT' 

IF(OUTBNK.EQ. 'RIGHT')  REFLD=' RIGHT' 
38    CONTINUE 

X(I)=X(I)*FMET 

Q(I)=Q(I)*FMET**3 

KYZ=NYZ-1 

JP=0 
99    READ(5,10)CPARM 

CBKGKBKG 

CEFL1=CEFL 

READ(5,*)  CBKG,CEFL,FLAGDY,XRNG,CIPT,TCORF 

WRITE(6,20)CPARM,CBKG,CEFL,FLAGDY,XFNG,CIPT,TC0RF 
20    F0RMAT(A,5X,6F10.5) 

IF(CBKG.EQ.-999.0.0R.JP.GE.5)G0  TO  150 

JP=JP+1 

READ(5,*)(C0NL(I,J),J=1,NYZ) 
C 

C     EXPRESS  DATA  W.R.T  ORIGIN  AT  OUTFALL  BANK  SIDE  OF  TRANSECT 
C 

CMX(JP)   =-1.0E+10 

DO  130  J=1,NYZ 

IF(0UTBNK.EQ. 'LEFT'. AND. REFLD.EQ. 'LEFT')  GO  TO  40 

IF(OUTBNK.EQ. 'RIGHT'. AND. REFLD.EQ. 'RIGHT')  GO  TO  40 

K=NYZ-J+1 

GO  TO  42 
40    K=J 
42    CONTINUE 

IF(JP.GE.2)  GO  TO  120 

Y(I,K)=YL(I,J)*FMET 

YW(K)=Y(I,K) 

Z(I.K)=ZL(I.J)*FMET 

IF(MVEL.EQ.99)  VEL( I,K)=UVL(I ,J)*FMET 
C 

C     CALCULATE  DYE  CONC'NS  WHEN  FLUORESCENCE  VALUES  ARE  INPUT 
C 

IF(FLAGDY.GE.1.)C0NL(I,J)=XRNG*(TC0RF*C0NL(I,J)-CIPT) 
120   CONTINUE 

CONC(I,J)=CONL(I,J) 
C 

C     SET  NEGATIVE  VALUES  TO  ZERO 
C 

IF(CONC(I.K).LT.O.O)  CONC(I .K)=0.0 
C 

C     FIND  PEAK  CONC.  &  ITS  POSITION  AT  TRANSECT. 
C 

IF(CMX(JP).GE.CONC(I,K))  GO  TO  130 

KP=K 

CMX(JP)=CONC(I.K) 
130   CONTINUE 
C 

BKFX=CBKG*(Q(I)-QEFL) 

EFLX=CEFL*QEFL 

TFLX=BKFX+EFLX 

IF(I.GT.1)WRITE(6.198) 
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198   F0RMAT(1H1/) 

WRITE(6,200)TRNSCT,X(I).JP.CPARM,Q(I),CBKG,QEFL,CEFL,BKFX,EFLX, 

*TFLX 
200   F0RMAT(/5X,A,2X,F8.1,'  METERS  FROM  OUTFALL '/lOX, ' PARAMETER' , 

llX.Il,':  ', A/5X,'QRIVER=',F9. 3, 5X,' BACKGROUND  CONC.=' .F8.3/5X, 

2'QEFL  =',F9. 3, 5X, 'EFFLUENT  CONCN.= ' ,F9.3/5X, 'UPSTREAM  FLUX=  ', 

3F10. 2, 3X. 'EFFLUENT  FLUX=    ' ,F8.2,3X, 'TOTAL  FLUX  =    '.F9.2) 
IF(MVEL.NE.99)WRITE(6,210) 
210       F0RMAT(10X, 'VELOCITIES  SIMULATED  FROM  RESISTANCE   EQN.'/) 

IF(MVEL.EQ.99)WRITE(6,212) 
212   F0RMAT(10X, 'MEASURED  VELOCITIES  CORRECTED  TO  GET  Q=SUMQ(NYZ) ' /) 

WRITE(6,202) 
202        F0RMAT(5X. 'Y' ,6X, 'Z' ,5X. 'VEL' ,5X. 'CONC ,6X. 'SUMA' ,5X, 'SUMQ' . 
*5X. 'SUMF' ,4X. 'Y/B' ,4X. 'QY/QT' .3X, 'C/CAVG' ,3X. 'C/CATRN'/) 

IF(JP.GE.2)  GO  TO  204 
C 

C  COMPUTE  AREA  &  DISCHARGE. 
C 

SUMA(1)=0. 

SUMQ(1)=0. 

DO  22  J=1,KYZ 

JJ=J+1 

DELA(J)=0.5*(Y(I.JJ)-Y(I,J))*(Z(I,JJ)+Z(I.J)) 

SUMA(JJ)=SUMA(J)+DELA(J) 

IF(MVEL.NE.99)  GO  TO  22 

DVEL(J)=0.5*(VEL(I,JJ)+VEL(I,J)) 

DELQ(J)=DELA(J)*DVEL(J) 

SUMQ(JJ)=SUMQ(J)+DELQ(J) 
21  CONTINUE 

ZAV=SUMA(NYZ)/Y(I,NYZ) 

ZZAV(I)=ZAV 

IF(MVEL.NE.99)  U( I)=Q(I )/SUMA(NYZ) 

IF(MVEL.EQ.99)U(I)=SUMQ(NYZ)/SUMA(NYZ) 
C   VELOCITY  SIMUUTION  USING  RESISTANCE(EG.  MANNING'S)  EQN. 

IF(MVEL.EQ.99)  GO  TO  106 

00  104  J-1,NYZ 
104   VEL(I,J)=F1*U(I)*(Z(I.J)/ZAV)**F2 
C 

C  ESTIMATE  DISCHARGE  FROM  SIMULATED  VEL.  DISTR'N. 
C 

SUMQ(1)=0.0 

DO  108  J=1,KYZ 

JJ=J+1 

DVEL(J)=0.5*(VEL(I.JJ)+VEL(I,J)) 

DELQ(J)=DELA(J)*DVEL(J) 

108  SUMQ(JJ)=SUMQ(J)+DELQ(J) 
106   CONTINUE 

C 

C   VELOCITY  CORRECTION  TO  CONFORM  WITH  SUMQ(NYZ)=Q(I ) . 

DO  109  J=1,NYZ 

109  VEL(I.J)=VEL(I.J)*Q(I)/SUMQ(NYZ) 
SUMQ(1)=0.0 

DO  110  J=1,KYZ 

JJ=J+1 

DVEL(J)=0.5*(VEL(I.JJ)+VEL(I,J)) 

DELQ(J)=DELA(J)*DVEL(J) 
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110   SUMQ(JJ)=SUMQ(J)+DELQ(J) 

U(I)=SUMQ(NYZ)/SUMA(NYZ) 
C 

C  SHAPE-VELOCITY  FACTOR. 
C 

CALL  CADIS(Z.VEL,ZAV,U,Q,DELQ,I,NYZ.SHAPE) 
C 

C   COMPUTE  FLUX  OF  TRACER  OR  POLLUTANT. 
C 
204   SUMF(1)=0. 

ARCY=0. 

DO  132  J=1,KYZ 

JJ=J+1 

IF(J.GT.KP.AND.CONC(I,J).LE.O.O)CONC(I,JJ)=0.0 

DC0NC(J)=0.5*(C0NC(I,JJ)+C0NC(I,J)) 

DCY=OCONC(J)*(Y(I.JJ)-Y(I,J)) 

ARCY=ARCY+DCY 

FLUX(J)=DCONC(J)*DELQ{J) 

SUMF(JJ)=SUMF(J)+FLUX(J) 

UNIF(J)=FLUX(J)/(Y(I,JJ)-Y(I.J)) 
132   CONTINUE 
C 

C   COMPUTE  AVG.  CONC'NS  IN  RIVER  AT  OUTFALL  &  TRANSECT. 
C 

CAVG=CEFL*QEFL/Q(I) 

CATRN=SUMF(NYZ)/SUMQ(NYZ) 
C 
C   PRINT  OUTPUT  MATRIX  OF  Y.Z.VEL.SUMQ  &  SUMF 

c 

DO  134  J=1,NYZ 
RYB=Y(I,J)/Y(I,NYZ) 
RQ=SUMQ(J)/SUMQ(NYZ) 
RC=CONC(I,J)/CAVG 
RCTRN=CONC(I,J)/CATRN 

WRITE(8.220)Y(I,J).Z(I,J),VEL(I,J),C0NC(I,J),SUMA(J), 
*SUMQ(J) ,SUMF(J) ,RYB,RQ,RC,RCTRN 
134   WRITE(6,220)Y(I.J),Z(I,J),VEL(I,J).C0NC(I.J).SUMA(J), 

*SUMQ(J),SUMF(J),RYB,RQ.RC,RCTRN 
220   F0RMAT(1X,F7.2,1X,F6.2,2X,F5.2,2X,3(F7.2,2X),F8.2,4(2X,F6.3)) 
FXTR=BKFX+SUMF(NYZ) 
WRITE(6,224)CAV6,CATRN,FXTR 
224   FORMAT(  /5X,'AVG.  CONC.  JUST  BELOW  OUTFALL,  CAVG= ' ,F8.3/5X, 'AVG.  C 
lONC.  AT  THE  TRANSECT,  CATRN  =' ,F8.3,5X, 'TOTAL  FLUX  AT  TRANSECT=  ' 
2,F9.2) 
IF(JP.LE.1)WRITE(6.312)  ZAV,U(I ) .SHAPE 
312   FORMAT (5X, 'MEAN  DEPTH=' ,F6.3,5X, 'MEAN  VELOCITY= ' ,F6.3,5X, 'SHAPE-VE 
*LOCITY  FACT0R=',F6.3) 
IF(JP.LE.1)WRITE(8,1122)ZAV,U(I) 
1122  FORMAT(2X,F10.5,3X,F10.5) 
C 

C      VARIANCE  COMPUTATION  *** 
C   VARIANCE  FROM  PEAK  CONC'N. 

c 

RLCN=CMX(JP)/ARCY 

VCMX(JP)  =1.0/(6.2836*RLCN*RLCN) 
C 
C  VARIANCE  FROM  2ND  MOMENT  OF  C-Y  DISTR'N. 


MIXANDAT  NOV.    14,    1986       15:54  PAGE     5 


CALL  VARANC(DCONC.NYZ.YW,VCN,ARCY,I,JP,KP,ZAV,SBS,SHS) 
C 

C  VARIANCE  FROM  2ND  MOMENT  OF  UNIT  FLUX  DISTR'N. 
C 

FLXPK=-1.0E+10 

DO  142  J=1.KYZ 

IF(FLXPK.GE.FLUX(J))  GO  TO  142 

KF=J 
142   CONTINUE 

ARUF=SUMF(NYZ) 

CALL  VARANC(UNIF,NYZ,YW,VUF,ARUF,I,JP,KF,ZAV,SBF,SHF) 
C 

C  VARIANCE  OF  C-Q  DISTRIBUTIONS 
C 

CALL  VARANC(DCONC,KYZ,SUMQ,VCQ,ARUF,I,JP,KF,ZAV,SBQ,SHQ) 

SQQ(JP)=VCQ(JP)/SUMQ(NYZ)**2 

RLCQ=CMX(JP)/ARUF 

VPQ(JP)=1.0/(6.2836*RLCQ*RLCQ) 

GO  TO  99 
150   CONTINUE 

WRITE(6.230)TRNSCT 
230   F0RMAT(/5X,A,':  VARIANCE  FROM  DIFFERENT  METHODS: 74X, 'PARAMETER' 
*, 4X , ' VCMAX ' ,7X, ' VCN ' ,8X , ' VUF ' ,8X, ' VCQ ' ,8X , ' VPQ ' /) 

DO  102  N=1,JP 

WRITE (6.234)  N,VCMX(N) ,VCN(N) ,VUF(N) ,VCQ(N) ,VPQ(N) 
234   F0RMAT(7X,I2,3X,5(F10.2.1X)) 
102   CONTINUE 

RXH=X(I)/ZAV 

RXB=X(I)/Y(I,NYZ) 

WRITE(6,240)RXB,RXH 
240   F0RMAT(///5X,'N0NDIMENSI0NAL  VARIANCE ' .5X, ' X/B= ' ,F8. 2, 5X, 'X/H=' , 
1F9. 1//2X . ' PARAMETER ' ,4X , ' VCN/BB ' , 5X , ' VCN/HH ' . 5X , ' VUF/BB ' , 5X , 
2'VUF/HH',5X,'VCQ/QQ'/) 

WRITE(9,3032)SBS(IP),SHS(IP),SBF(IP).SHF(IP),SQQ(IP),RXB.RXH 
3032  F0RMAT(2X,2(F8.4,2X,F9.3,3X),F9.4,2X,F8.4,2X,F10.4) 

DO  100  N=1,JP 

WRITE(6,242)N,SBS(N).SHS(N),SBF(N),SHF(N),SQQ(N) 
242   F0RMAT(5X,I2,6X,2(F8.4,2X,F9.3,3X),F9.4) 
100   CONTINUE 
9999  CLOSE (5) 

CLOSE (6) 

CL0SE(9) 

CL0SE(8) 
C 

C  HERE    IS  WHERE  THE   INPUT  FILE  FOR  MIXCALBN   IS  WRIHEN 

C  TO    INCAL.DAT 

C 

WRITE(7,10)TITLE 

BPWR=0.05 

HPWR=0.5 

UPWR=0.45 

THETA1=1.03 

WRITE(7.9000)Q(1),BPWR.HPWR,UPWR.THETA1 
9000     F0RMAT(2X,F10.2.2X,F4.2,2X,F4.2.2X,F4.2,2X,F6.3,2X,F7.2.2X,F10.8) 

QRUP=Q(1)-QEFL 

WRITE(7,9001)QRUP,QEFL,CEFL1,CBKG1.TMP1 


MIXANDAT  NOV.  14,  1986   15:54  PAGE   6 


9001  F0RMAT(2X.F10.2,2X,F10.2,2X,F8.3.2X,F8.3,2X,F8.3) 
IF(OUTBNK.EQ. 'RIGHT' )QCP=0. 
IF(0UTBNK.EQ.'LEFT')QCP=Q(1) 

NYZ=10 
WRITE(7,9002)NTR,NYZ,QCP 

9002  F0RMAT(3X,I3.3X,I3,3X.F10.3) 
DO  9004  III=1,NTR 
BS1=Q(1)/(ZZAV(III)*U(III)) 
WRITE(7,9003)X(III),BS1.ZZAV(III),U(III) 

9003  F0RMAT(2X.F10.2,2X,F10.2.2X.F6.2,2X,F6.2) 

9004  CONTINUE 
C 

C 

CL0SE(7) 

CALL  REGVAR 

STOP 

END 
C 

C  COMPUTATION  OF  SHAPE-VELOCITY  FACTOR. 
C 

SUBROUTINE  CADIS(H,V,HA,VA,QR,DLQ.L,N,SHP) 

DIMENSION  H(8,60) .V(8,60) ,QR(8) ,DLQ(60) ,VA(8) 

SUMD=0. 

NS=N-1 

DO  12  J=1,NS 

JJ=J+1 

HR=((H(L,JJ)+H(L,J))/HA)**2 

VR=((V(L,JJ)+V(L,J))/VA(L)) 

BHUY=DLO(J)*HR*VR 


12 

SUMD=SUMD+BHUY 

SHP=SUMD/(8.0*QR(L)) 

RETURN 

END 

C 

C 

SUBROUTINE  TO  COMPUTE  VARIANCE 

c 

VALUES  FOR  BANK  OUTFALL  CASE 

SUBROUTINE  VARANC(P.NY,R,V,SUMD,IC, J,KI .Z.SB.SH) 

DIMENSION  P(60).R(60),V(5).SB(5),SH(5) 

SUMN=0. 

K=NY-1 

M=l 
10    L=M+1 

IF(M.GT.KI. AND. P(M).LE. 0.00011)00  TO  12 

SUMN=SUMN-K).25*P(M)*(R(L)-R(M))*(R(M)+R(L))**2 

IF(M.GE.K)GO  TO  12 

M=M+1 

GO  TO  10 
12    V(J)=SUMN/SUMD 

SB(J)=V(J)/R(NY)**2 

SH(J)=V(J)/Z**2 

RETURN 

END 
C 
C 
C 

C     PROGRAM  REGVAR. FOR 
C     TAKES  VALUES  COMPUTED  BY  MIXANDAT. FOR  TO  DETERMINE  IF  THE 
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11 
2 

3 


102 
4 
C 
C 


c 
c 

78 
12 
13 
14 
15 


16 
200 
C 
C 


NONDIMENSIONAL  PARAMETER  BETA  IS  DEPTH  OR  WIDTH  DEPENDANT. 
IT  REGRESSES  THE  VARIOUS  VARIANCE  ESTIMATES  AGAINST  WIDTH, 
DEPTH  AND  FLOW  VARIABLES  AND  CHOOSES  BETA  BY  CONSIDERING  THE 
EQUATION  WITH  THE  SMALLEST  RESIDUAL  ERROR. 

SUBROUTINE  REGVAR 

DIMENSION  VCB(20) ,VCH(20) ,VUB(20) ,VUH(20) ,VCQ(20) ,8(20) .H(20) 
*.Q(20) 

CHARACTER*30  SNAME 

0PEN(5,FILE='0UTVAR.DAT',STATUS='NEW') 

0PEN(6,FILE='BETA.DAT',STATUS='0LD') 

WRITE(5,1) 

FORMATC   ENTER  THE  NUMBER  OR  TRANSECTS  USED:   '\) 

READ(6,*)NTR 

WRITE(5,11)NTR 

F0RMAT(/,I2) 

WRITE(5,2) 

FORMAT(/'   ENTER  THE  FOLLOWING  VARIANCES  AND  PARAMETERS'//) 

WRITE(5,3) 

FORMAT('TRAN     VCN/BB    VCN/HH    VUF/BB    VUF/HH    VCQ/QQ 
•^  BB    HH  ') 

DO  4  J=1,NTR 

WRITE(5,5)J 

FORMAT (2X, I 2\) 

READ(6,*)VCB(J).VCH(J),VUB(J),VUH(J),VCQ(J).B(J),H(J) 

WRITE(5,102)VCB(J),VCH(J),VUB(J),VUH(J),VCQ(J),B(J),H(J) 

F0RMAT(2X,F7.4,2X,F7.2,2X,F7.4,2X,F7.2,2X.F7.4,2X.F7.3.2X,F7.2) 

CONTINUE 


DO  200  KK=1,4 
NNY=NTR-KK+1 
IF(NNY.LE.2)G0  TO  200 


WRITE(5,78)NNY 

F0RMAT(/'   REGRESSION  COEFFICIENTS  FOR  ',12,'  TRANSECTS') 

CALL  REGO(VCB,B,RCB,ACB,BCB,NTR,NNY) 

WRITE(5,12)ACB,BCB,RCB 

FORMATC    ACB=' ,F10.5,'   BCB= ' ,F10.5, '    RCB=',F6.4) 

CALL  REGO(VCH,H,RCH,ACH,BCH,NTR,NNY) 

WRITE(5,13)ACH,BCH,RCH 

FORMATC    ACH=',F10.5,'   BCH= ' ,F10.5, '    RCH=',F6.4) 

CALL  REGO(VUB,B,RUB,AUB,BUB,NTR,NNY) 

WRITE(5,14)AUB,BUB,RUB 

FORMATC    AUB=',F10.5,'   BUB= ' .F10.5, '    RUB=',F6.4) 

CALL  REGO(VUH,H,RUH,AUH,BUH,NTR,NNY) 

WRITE(5.15)AUH,BUH,RUH 

FORMATC    AUH=' ,F10.5,'    BUH=' ,F10.5, '    RUH=',F6.4) 

CALL  REGO(VCQ,B.RCQ,ACQ,BCQ,NTR,NNY) 

WRITE(5.16)ACQ,BCQ.RCQ 

FORMATC    ACQ=' ,F10.5,'   BCQ=' .F10.5, '    RCQ=',F6.4) 

CONTINUE 


BETA=BC0/2. 
WRITE(5,2206)BETA 
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CALL  CLRSCN 

WRITE(*,5558) 
5558     FORMAT(///) 

WRITE(*,2206)BETA 
2206     FORMATC  THE  CALCULATED  VALUE  OF  BETA   IS:    ',F10.7) 

WRITE(*,9912) 
9912     FORMAT(//'  REMEMBER  THIS  VALUE  FOR  USE   IN  THE  CALIBRATION  PROG 

*RAM  "MIXCALBN"') 

C 
C 

WRITE(5,9911)CHAR(26) 

WRITE(6,9911)CHAR(26) 
9911     FORMAT(A) 

CL0SE(5) 

CLOSE (6) 

RETURN 

END 
C 
C 
C 
C 

SUBROUTINE   RE60(V,X,R,A,B,N,NNY) 

DIMENSION  V(20),X(20) 

SMX=0. 

SMV=0. 

SMXX=0. 

SMVV=0. 

SMXV=0. 

A=0. 

B=0. 

N0N=N-NNY+1 

DO   50    I=NON,N 

SMX=SMX+X(I) 

SMV=SMV+V(I) 

SMXX=SMXX+X(I)*X(I) 

SMVV=SMVV+V(I)*V(I) 

SMXV=SMXV+X(I)*V(I) 
50  CONTINUE 

B1=SMXV-SMX*SMV/N 

B2=SMXX-SMX*SMX/N 

B=B1/B2 

A=(SMV-B*SMX)/N 

R1=B1 

R2=SQRT(B2) 

R3=SMVV-SMV*SMV/N 

R3=SQRT(R3) 

R=R1/(R2*R3) 

RETURN 

END 
C 
C 
C 

SUBROUTINE  CLRSCN 

WRITE(*,101)CHAR(27),'[2J' 
101   F0RMAT(1X,A,A\) 

RETURN 

END 
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PLTANDAT 

NOV. 

14, 

1986 

15:55 

PAGE 

1 

SOEBUG 

**** 

******< 

r***********  »»< 

r*  *♦■*-*•*-**-**-*  * 

** 

c 

C  PROGRAH  PLTANOAT.FOR  * 

C  * 

C  PLOTS  THE  OUTPUT  FROM  PROGRAM  MIXANDAT  AND  * 

C  PRODUCES  INPUT  FILES  FOR  MIXCALBN.FOR  AND  MIXCADIF.FOR  * 

C  WRIHEN  BY  R.  JARVIS  * 

c 

C     GORE  &  STORRIE  1986  * 

COMMON/A/  Y(35.10).Z(35.10).V(35.10).C(35.10),QY(35.10) 

COMMON/B/  A1(35).A2(35).A3(35).A4(35).A5(35).DAVA.VAVA 

COMMON/D/  A1R(35) ,A2R(35) .A3R(35) .A4R(35) .A5R(35) 

COMMON/C/  DAV(10),VAV(10) 

COMMON/WW/  YQ(35).ZZ(35).VV(35),CC(35) 

COMMON/MAX/  YMAX.ZMAX.CMAX.VMAX 

COMMON/NAME/  TITLE 

INTEGER  N(35) 

CHARACTER*20  FILIN 

CHARACTER*80  TITLE 

CHARACTER*30  SNAME 

YMAX=0. 

ZMAX=0. 

CMAX=0. 

VMAX=0. 

0PEN(3,FILE='PL0T0UT.DAT',STATUS='NEW') 

0PEN( 5 , FI LE= ' PLOTAN .DAT ' , STATUS= ' OLD ' ) 

0P€N(6,  FI  LE=' SCALE. DAT ',STATUS=' NEW) 

READ(5.333)TITLE 

WRITE(3,333)TITLE 
333   FORMAT(A) 

READ(5,*)NTR 

WRITE(3.525)NTR 
525   F0RMAT(2X,I2) 

NYZ=10 

WRITE(3,252)NYZ 
252       F0RMAT(2X,I2) 

DO  2    1=1, NTR 

J=0 

N(I)=0 
99         J=J+1 

N(I)=N(I)+1 

?r/5!^?!^ii'^^'^(^'^)'^^^'^)'C('J'^)'Dl.D2.D3.D4.QY(J.I) 

IF(Y(J,I).GT.YMAX)YMAX=Y(J,I)  "^   v    .    ^ 

IF(Z(J,I).GT.ZMAX)ZMAX=Z(J,I) 

IF(C(J,I).GT.CMAX)CMAX=C(J.I) 

IF(V(J.I).GT.\/MAX)VMAX=V(J,I) 

IF(QY(J.I).NE.1.0)G0  TO  99 

READ(5,*)DAV(I),VAV(I) 

2  CONTINUE 
CL0SE(5) 

C 

WRITE(*.3) 

3  F0RMAT(/'        PLOT  LOCATION       SCREEN=1     PL0nER=2     PRINTER=3      '\) 

READ(*,*)IPL 
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IF(IPL.EQ.3)CALL  PL0TS(0,0,11) 
IF(IPL.EQ.2)CALL  PL0TS(0,9600,80) 
IF(IPL.EQ.1)CALL  PL0TS(0,0,99) 

DO  8  1=1,35 

A1(I)=0. 

A2(I)=0. 

A3(I)=0. 

A4(I)=0. 

A5(I)=0. 
8     CONTINUE 

DO  6  1=1, NTR 

NN=N(I) 

DO  7  J=1,NN 

A1(J)=Y(J.I) 

A2(J)=Z(J,I) 

A3(J)=V(J,I) 

A4(J)=C(J,I) 

A5(J)=QY(J,I) 
7     CONTINUE 

DAVA=DAV(I) 

VAVA=VAV(I) 

CALL  PLOOT(NN,I,IPL) 

WRITE(6,*)(YQ(J).J=1,11) 
6     CONTINUE 

CALL  PL0T(0.,0.,999) 
900   CONTINUE 

CL0SE(3) 

CL0SE(6) 

STOP 

END 
G 
C 

SUBROUTINE  PLOOT(NN,ITR,IPL) 

COMMON/B/  Al(35) ,A2(35) ,A3(35) ,A4(35) ,A5(35) .DAVA.VAVA 

COMMON/D/  A1R(35) ,A2R(35) ,A3R(35) .A4R(35) ,A5R(35) 

COMMON/MAX/  YMAX,ZMAX,CMAX,VMAX 

COMMON/NAME/  TITLE 

COMMON/WW/  YQ(35),ZZ(35),VV(35).CC(35) 
•  CHARACTER*80  TITLE 
■  SCY=YMAX/16. 

SCZ=ZMAX/4. 

SCC=CMAX/4. 

SCV=VMAX/4. 

DO  99  J=1,NN 

A2(J)=-A2(J) 
99    CONTINUE 

IF(IPL.EQ.3)CALL  FACT0R(.25) 

IF(IPL.EQ.2)CALL  FACT0R(.25) 

IF(IPL.EQ.1)CALL  FACT0R(.2) 

CALL  SIMPLX 

CALL  STAXIS(.2,.30,.05,.07.1) 
C 

CALL  PL0T(15.,20..-3) 

A1(NN+1)=0. 

A1(NN+2)=SCY 

A2(NN+1)=-ZMAX 

A2(NN+2)=SCZ 
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A3(NN+1)=0. 

A3(NN+2)=SCV 

A4(NN+1)=0. 

A4(NN+2)=SCC 

A5(NN+1)=0. 

A5(NN+2)=.2 

c 

CALL  SYMB0L(-14.5,5.0..4,TITLE,0.,80) 

CALL  SYMB0L(-14. 5, 4.0, .3. 'CROSS-SECTIONAL  VALUES  FOR  TRANSECT' ,0. , 
♦35) 

SITR=FLOAT(ITR) 

CALL  NUMBER(-4.0,4.,.3,SITR,0.,-1) 

CALL  SYMB0L(-14.5,3.,.3,'(  WITH  10  STREAM  TUBES  )',0.,24) 
C 

CALL  AXIS(0.,0.,'FLOW  FRACTION' ,13,5. ,90. ,0. , .2) 

CALL  STAXIS(0.,.3,.05,.07,1) 

CALL  AXIS(0.,0.,'   '.-2,16.,0.,Al(NN+l),Al(NN+2)) 

CALL  STAXIS(.2,.3,.05,.07,1) 

CALL  C0L0R(4,IERR) 

CALL  LINE(A1,A5,NN, 1,1,1) 

CALL  C0L0R(0,IERR) 
C 

111=1 

DO  939  11=2,11 

F=FL0AT(II-1)/10. 
494   IF(F.GT.A5(III))III=III+1 

IF(F.GT.A5(III))G0  TO  494 

L=III-1 

Y1=A1(L) 

Y2=A1(III) 

DELY=Y2-Y1 

AAQ=A5(III)-F 

AAAQ=F-A5(L) 

DELQ=AAO+AAAQ 

YQ(II)=AAAQ*DELY/DELQ+Y1 

SS=YQ(II)/SCY 
C 

Z1=A2(L) 

Z2=A2(III) 

DELZ=Z2-Z1 

ZZ(II)=AAAQ*DELZ/DELQ+Z1 

V1=A3(L) 

V2=A3(III) 

DELV=V2-V1 

VV(II)=AAAQ*DELV/DELQ+V1 

C1=A4(L) 

C2=A4(III) 

DELC=C2-C1 

CC(II)=AAAQ*DELC/DELQ+C1 

c 

11=1 

IF(II.EQ.2)WRITE(3,444)I1,A1(1),A2(1),A3(1).A4(1) 

ZZZ=-ZZ(II) 

WRITE(3.444)II.YQ(II).ZZZ,VV(II),CC(II) 
444       F0RMAT(2X.I2.F7.2,2X,F7.3.2X,F7.3,2X,F7.4) 
C 

FF=5.*F 
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CALL  STDASH(.2,.2) 

CALL  PL0TD(0.,FF.3) 

CALL  PL0TD(SS,FF,2) 

CALL  PL0TD(SS,0..2) 

CALL  PL0TD(SS.-1.5,3) 

CALL  PLOTD(SS.-5.5.2) 

CALL  PLOTO(SS,-7.0,3) 

CALL  PL0TD(SS.-11.0,2) 

CALL  PL0TD(SS,-12.5,3) 

CALL  PL0TD(SS.-16.5,2) 
939   CONTINUE 
C 

CALL  PL0T(0.,2.,-3) 

CALL  PL0T(-10.,-2.,3) 

CALL  PL0T(-10.,-18.,2) 

XST=-10.25 

XEN=-9.75 

DO  9938  IP=1.11 

QQ=FL0AT(IP-1)/10. 

YPL=-18.+1.6*FL0AT(IP-1) 

CALL  PL0T(XST,YPL.3) 

CALL  PL0T(XEN,YPL,2) 

XNUM=XST-1.0 

XMUN=XEN+.5 

YYQ=YQ(IP) 

CALL  NUMBER(XNUM,YPL-.15,.2,QQ,0.,1) 

CALL  NUMBER(XMUN,YPL-.15,.2,YYQ,0.,1) 
9938  CONTINUE 

CALL  SYMBOL (-13., -14.,. 3, 'PARTIAL-  DISCHARGE ' ,90. , 17) 

CALL  SYMB0L(-7. 0,-14., .3, 'LATERAL  DISTANCE  (M)',90.,20) 

CALL  SYMB0L(-14., -1.0, .3, 'STREAM  TUBE/DISTANCE  SCALE ' ,0. ,26) 

CALL  PL0T(0.,-2.,-3) 
C 

CALL  COLOR(0,IER) 

CALL  PLOT(0. ,-5.5,-3) 

CALL  AXIS(0.,0., 'DEPTH' ,5,4. ,90. ,A2(NN+1) ,A2(NN+2) ) 

CALL  STAXIS(0.,.3,.05,.07,1) 

CALL  AXIS(0.,0.,'   ',-2,16.,0.,Al(NN+l),Al(NN+2)) 

CALL  STAXIS(.2,.3,.05,.07,1) 

NNN=NN+2 

CALL  C0L0R(4,IERR) 

CALL  LINE(A1,A2,NN, 1,1,1) 

CALL  C0L0R(0,IERR) 

SCZ=A2(NN+2) 

DA=4.-DAVA/SCZ 

CALL  PL0T(0.,DA.3) 

CALL  PL0T(17.,DA,2) 

CALL  SYMB0L(17.5,DA-.5,.35.'AVE.  DEPTH' ,0. ,10) 

CALL  NUMBER(17.75,DA,.35,DAVA,0.,2) 
C 

CALL  PLOT(0. ,-5.5.-3) 

CALL  AXIS(0., 0.,' VELOCITY' ,8.4. ,90., A3(NN+l),A3(NN+2)) 

CALL  STAXIS(0...3,.05..07,1) 

CALL  AXIS(0.,0.,'   ',-2,16.,0..Al(NN+l),Al(NN+2)) 

CALL  STAXIS(0.2,.3,.05..07.1) 

CALL  C0L0R(4.IERR) 

CALL  LINE(A1.A3,NN. 1.1.1) 
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CALL  COLOR(0,IERR) 

SCV=A3(NN+2) 

VA=VAVA/SCV 

CALL   PL0T(0.,VA,3) 

CALL   PL0T(17..VA,2) 

CALL  SYMB0L(17.5.VA-.5,.35.'AVE.    VELOCITY ' ,0. ,13) 

CALL  NUMBER(17.75,VA,.35,VAVA,0.,2) 

CALL  PLOT(0. ,-5.5,-3) 

CALL  AXIS(0.,0., 'CONCENTRATION' ,13,4. ,90. .A4(NN^1) ,A4(NN+2) ) 

CALL  AXIS(0.,0.,' LATERAL  DISTANCE   Y ' ,-18, 16. ,0. ,A1(NN+1) ,Al(NN+2) ) 

CALL  C0L0R(4,IERR) 

CALL  LINE(A1,A4,NN, 1,1,1) 

CALL  C0L0R(0,IERR) 

CALL  PL0T(0.,0.,-999) 

RETURN 

END 
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SDEBUG 

C******^****  ***************************************************  ****** 

C     PROGRAM  NAME  MIXCALBN. FOR  DEVELOPED  FROM  MIXCALBN  TO  RUN  ON     * 

C     A  MICROCOMPUTER  * 

C****  ************************************************************  ****** 

C  PROGRAM  NAME:  MIXCALBN  *  *  STREAMTUBE  MODEL  FOR  PIPE  OUTFALL  * 

C  DEVELOPED  BY  T.  P.  H.  GOWDA,  WATER  RESOURCES  BRANCH,  MOE.  * 

C  THIS  PROGRAM  PREDICTS  LAT ' L  &  LONG'L  DISTRN.  OF  CONSERVATIVE  * 

C  AND  NONCONSERVATIVE  MATERIALS  DISCHARGED  INTO  A  RIVER  FROM  * 

C  A  PIPE  OUTFALL  LOCATED  AT  BANK  OR  IN  RIVER(VERT.  LINE  SOURCE).  * 

C  PROGRAM  MODIFIED:  JUNE  1983  FOR  DILUTION  FACTOR  AND  TERMINATE  * 

C  CALC'NS  IF  CONCN  <  1.0E-04.  * 

C  * 

c  * 

C     GORE  &  STORRIE  1986  ♦ 

C**  ***♦**♦**♦****************♦♦♦*«♦*********♦*****♦****♦  **********^^** 

DIMENSION  C(50,102),CUI(50,102) 

REAL*8  X(50),XX(50),P1,P2,P3,P4.T1,T2,T3,T4,RKS(50).QY(502). 
*THETA,BPWR,HPWR,UPWR,QRTO,QT,RBT,QRS,QRUP,QEFL,RBK,CTDP,PHDR, 
*RF(50),BS(50),HS(50),US(50),BETA(50),BW,HW,UW,RKT,R,PHI,TMP, 
*B(50),H(50),U(50),BSUM(50),T0T(50),V0L(50), TEMPS, PAX1,PAX2, 
*A3,QCP.DELQ 

CHARACTER*80  TITLE 

CHARACTER*20  FILIN,FIL0UT 
C 

C     INPUT  DATA 

C     MIXCALBN  READS  FILE  "PINCAL.DAT"  FROM  SUBROUTINE  SETUP 
C     THE  FIRST  TIME  THROUGH  AND  MAKES  FILE  "PPINCAL.DAT"  TO 
C     BE  USED  ON  SUBSEQUENT  RUNS. 
C 

CALL  CLRSCN 

WRITE(*,1290) 
1290  F0RMAT(/////'     FIRST  TIME  THROUGH  MIXCALBN    1=YES   0=N0') 

READ(*.*)ISET 

IF(ISET.EQ.I)CALL  SETUP 

0PEN(5,FILE='PPINCAL.DAT',STATUS='0LD') 

0PEN(4,FILE='CAL0UT.DAT',STATUS='NEW') 

0PEN(6,FILE='PLCALPC' ,STATUS='NEW') 
C 

C     ICAL  IS  THE  FLAG  IN  PLCALPC  THAT  INDICATES  A  PIPE  OUTFALL 
C     TO  THE  PLOTTING  PROGRAM  CONMIX.FOR 

c 

ICAL=1 

WRITE(6,*)ICAL 

WRITE(4,2) 

2  FORMAT(/'      ENTER  TITLE  OF   STUDY') 
READ(5,3)   TITLE 
WRITE(6,3)TITLE 

3  FORMAT(A) 
WRITE(4,3)TITLE 

35         WRITE(4.4) 

4  FORMATC        ENTER  QRS,BPWR,HPWR, UPWR,THETA, TEMPS, RBK') 
READ(5,*)    QRS,BPWR,HPWR,UPWR,THETA, TEMPS, RBK 

WRITE (4,400)   QRS,BPWR,HPWR,UPWR,THETA, TEMPS, RBK 
400       F0RMAT(F7.2,2X.F5.2,2X,F5.2,2X,F4.2,2X,F6.3,2X,F6.2,2X.F4.2) 
43         WRITE(4,19) 
19  FORMATC        ENTER  DESIGN  CASE:   QRUP.QEFL,CEFL,CBKG,TMP' ) 
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READ(5,*)QRUP,QEFL,CEFL,CBK6,TMP 
WRITE(4.401)QRUP,QEFL.CEFL.CBKG.TMP 

401  F0RMAT(2X,F7.2.2X,F8.3,2X,F8.2.2X.F5.2,2X,F6.2) 
33    WRITE(4  5) 

5  FORMAT('   ENTER  NTR.NYZ.QCP' ) 
READ(5,*)  NTR.NYZ.qCP 
WRITE(6.*)NTR 
WRITE(6.*)NYZ 
WRITE(4,402)NTR,NYZ,QCP 

402  F0RMAT(2X,I3.3X,I3,F5.1) 
WRITE(4,6)NTR 

6  FORMATC   ENTER  '.12,'  VALUES  OF  X,BS,HS,US') 
READ(5,*)  (X(I),BS(I),HS(I),US(I),I=1,NTR) 
READ(5,*)(RKS(I),I=1,NTR) 

DO  405  1=1, NTR 

WRITE(4,403)I,X(I).BS(I),HS(I),US(I) 
WRITE(6, 8234)1, X(I) 
8234  F0RMAT(2X,I3,2X,F8.2) 

403  F0RMAT(2X,I3,2X,F8.2,3X,F8.2,3X,F6.2,3X,F7.4) 

405  CONTINUE 

47    WRITE(4,55)NTR 

55    FORMATC    ENTER  ',12,'  VALUES  OF  DECAY') 

DO  408  1  =  1, NTR 

WRITE(4, 406)1, RKS(I) 

WRITE(6, 406)1, RKS(I) 

406  F0RMAT(2X,I3,F9.7) 

408  CONTINUE 
45    CONTINUE 

WRITE(*,7777) 
WRITE(4,7777) 

7777  F0RKAT(/'   ENTER  VALUES  OF  BETA'/) 
DO  409  1=1, NTR 

WRITE(*, 7778)1 
WRITE(4, 7778)1 

7778  FORMATC   TRANSECT  ',13,'  :  '\) 
READ(*,*)BETA(I) 
WRITE(4,411)BETA(I) 

411  F0RMAT(5X,F9.7) 
WRITE(6, 410)1, BETA(I) 

410   F0RMAT(2X,I3,F9.7) 

409  CONTINUE 
WRITE(4,52) 

52    FORMATC   ARE  YOU  CONSIDERING  UN-IONIZED  AhWONIA:  YES=1   N0=0') 

READ(5,*)  AMONIA,PH 

IF(AMONIA.EQ.O)GO  TO  101 

WRITE(4,102) 
102   FORMATC   ENTER  THE  PH  OF  UN-IONIZED  AhWONIA') 

WRITE(4,412)PH 

412  FORMATC   THE  PH  OF  UN- IONIZED  AMMONIA  IS   ',F3.1) 
101   CONTINUE 

C 

C  CALCULATE  FLOW  &  TEMP'R  PARAMETERS 

C 

DELTA=O.OOOI 

QT=QRUP+qEFL 

QRTO=QT/QRS 

CTDP=THETA**(TMP-TEMPS) 
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RBT=RBK*CTDP 

DELQ=QT/NYZ 

NQ=NYZ+1 

CA=(CEFL*QEFL  /QT) 

KCP=QCP/DELQ+1.5 

IF(KCP.LE.l)KCP=KCP+3 

DO  12  1=1, NTR 

BSUM(I)=0. 
12  TOT(I)=0. 

DO  14  1=1. NTR 
C 

C       CALCULATE   B.H.U   FOR  FLOW=QT.    FROM  LEOPOLD-MADDOCK  EQNS. 
C 

B(I)=BS(I)*QRTO**BPWR 

H(I)=HS(I)*QRTO**HPWR 

U(I)=US(I)*QRTO**UPWR 
C 

C        CALCULATE  WEIGHTED  MEAN  VALUES  BW.HW.UW  FROM  OUTFALL  TO  TRANSECT(I) 
C 

IF(I.GE.2)   GO  TO  60 

XX(1)=X(1) 

BW=B(1) 

HW=H(1) 

UW=U(1) 

BSUM(1)=B(1)*XX(1) 

V0L(1)=B(1)*H(1)*XX(1) 

T0T(1)=XX(1)/U(1) 

GO  TO  62 
60   11=1-1 

XX(I)=X(I)-X(I1) 

BSUM(I)=BSUM(I1)+0.5*XX(I)*(B(I1)+6(I)) 

V0L(I)=V0L(Il)+0.25*XX(I)*(B(Il)+B(I))*(H(Il)++i(I)) 

T0T(I)=T0T(I1)+XX(I)/U(I) 

BW=BSUM(I)/X(I) 

HW=VOL(I)/(X(I)*BW) 

UW=QT/(BW*HW) 
C 

C   CALCULATE  PRODUCT  FUNCTION  FOR  DECAY.  RF(I). 
C 

62  RKT=CTDP*RKS(I) 

CBKX=CBKG*DEXP( -RBT*TOT( I ) ) 

A3=(RKT*XX(I))/U(I) 

R=0EXP(-A3) 

IF(I.GE.2)   GO  TO  64 

RF(1)=R 

GO  TO  66 
64  RF(I)=RF(I1)*R 
66  CONTINUE 

PHI=BETA(I)*X(I)/BW 

PHDR=4.0*PHI 

CRPX=0.5*CA*RF(I)/DSQRT(3.1416*PHI) 

BGX=PHI*AL0G(1. /DELTA) 

SBG=SQRT(BGX) 

WRITE(4,40)    I.BETA(I).RKS(I) 
40         F0RMAT(/5X. 'TRANSECT:    ' ,I2.2X. ' BETA=' ,F9.6,2X. ' RKS= ' .F9.6.2X. 

♦/7X,'X'  .9X,'BW'  ,9X.'HW'  .BX.'UW) 

WRITE(4.23)    X(I). BW.HW.UW 
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23  F0RMAT(2X,4(F9.3,1X)) 

WRITE(4,42) 
42         F0RMAT(4X,'QY'.5X.'C(X.QY)',5X.'CUr,9X.'C/CA'.7X,'QY/QT'.6X, 
*'DIL  FAC7) 

DO   16  K=1,NQ 

QY(K)=FL0AT(K-1)*DELQ 

IF    (QY(K).GT.QT)   QY(K)=QT 

PAX1=(QY(K)-QCP)/QT 

PAX2=(QY(K)-K)CP)/QT 
C 

C     DETERMINE  NO.   OF   IMAGES  REQUIRED 
C 

ANl=(0.5*PAXl-SBG)-0.5 

AN2=(0.5*PAXl+SB6)+0.5 

AN3=-AN2 

AN4=-AN1 

NM1=IFIX(AN1) 

NM2=IFIX(AN2) 

NM3=IFIX(AN3) 

NM4=IFIX(AN4) 

NNl  =  l-fNM2+IABS(NMl) 

NN2=1+IABS(NM3)+IABS(NM4) 

IF(NN1.GE.NN2)NN=NN1+1 

IF(NN1.LT.NN2)NN=NN2+1 
C 

C  COMPUTATION  OF  CONC'N  DISTR'NS. 
C 

SUM=0. 

DO  32  J=1,NN 

N=J-1 

P1=(PAX1-2.*N)**2/PHDR 

P2=(PAX2+2.*N)**2/PHDR 

CALL  PDET(P1,T1) 

CALL  PDET(P2,T2) 

IF(N.LE.O)  GO  TO  30 

P3=(PAX1+2.*N)**2/PHDR 

P4=(PAX2-2.*N)**2/PHDR 

CALL  PDET(P3,T3) 

CALL  PDET(P4,T4) 

GO  TO  32 
30  T3=0. 

T4=0. 
32  SUM=SUM+T1+T2+T3+T4 

C(I,K)=CRPX*SUM+CBKX 
C 

C  CALCULATE  UN-IONIZED  AMMONIA  CONCENTRATIONS*  OPTIONAL  * 
C 

IF(AMONIA.LE.O)  GO  TO  15 

PK;^=0. 090 18+2729. 92/(TMP+273. 2) 

PF=PKA-PH 

PCTU=1./(1.+10.**PF) 

CUI(I,K)=C(I,K)*PCTU 

15  CONTINUE 

16  CONTINUE 
NQO=NQ 

18  CONTINUE 
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C  PRINT  OUTPUT 
C 

DO  20  K=1,NQ 

RC=C(I,K)/CA 

RQ=QY(K)/QT 

CNET=C(I,K)-CBKX 

I F(CNET.LE. 0.000009)  CNET=-CEFL 

DLF=CEFL/CNET 

IF(AMONIA.LE.O)  CUI(I,K)=0.0 

WRITE(6,25)K,QY(K),C(I,K),CUI(I,K),RC,RQ,DLF 
20    WRITE(4,25)K,QY(K),C(I,K),CUI(I.K),RC,RQ,DLF 
25    F0RMAT(I2,2X,F10.2,2X,4(F9.4,2X),F10.2) 
14    CONTINUE 
999   CLOSE (4) 

CL0SE(5) 

CL0SE(6) 

CALL  CLRSCN 

WRITE(*,1234) 
1234  F0RMAT(//'    THE  OUTPUT  FILE  FROM  MIXCALBN. FOR  IS  CALLED'//' 
*         "  CALOUT.DAT  "') 

STOP 

END 
C 

SUBROUTINE  PDET(P,T) 

REAL*8  P.T 

IF(P.GE.40.0)G0  TO  10 

T=DEXP(-P) 

GO  TO  12 
10  T=0.0 
12  CONTINUE 

RETURN 

END 
C 
C 

SUBROUTINE  SETUP 

DIMENSION  RKS(20),BS(20),ZAV(20),VAV(20),X(20),DY(15) 

CHARACTER*80  TITLE 

0PEN(7,FILE='PINCAL.DAT',STATUS='0LD') 

READ(7,1)TITLE 

1  FORMAT(A) 
READ(7,*)QRS,BP,HP,UP,THETA 
READ(7,*)QRUP.QEFL,CEFL,CBKG.TMP 
READ(7,*)NTR,NYZ.QCP 

DO  2  1=1, NTR 
READ(7,*)X(I).BS(I),ZAV(I),VAV(I) 

2  CONTINUE 
CALL  CLRSCN 
WRITE(*.4005) 

4005  FORMAT(/'     HYDRODYNAMIC  PARAMETER  ENTRY  AREA'/'    *****#*****i 

WRITE(*,3)BP,HP,UP 

3  FORMAT(/////'    THE  EXPONENTS  FOR  THE  LEOPOLD-MADDOCK  EQNS  ARE'// 
♦'   WIDTH  EXP:  ',F5.3/'  DEPTH  EXP:  ',F5.3/'   VEL.   EXP:  ',F5.3) 

WRITE(*,4) 

4  FORMAT(///'    DO  YOU  WISH  TO  CHANGE  THEM?   YES=1   N0=0   '\) 
READ(*,*)ICH 
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11  IF(ICH.EQ.1)THEN 

WRITE(*.5) 

5  FORMAT(/'      WIDTH  EXP=      '\) 
READ(*,*)BP 

WRITE(*,6) 

6  FORMAT(/'      DEPTH  EXP=      '\) 
READ(*,*)HP 

WRITE(*,7) 

7  FORMAT(/'      VEL.    EXP=        '\) 
READ(*.*)UP 

ENDIF 

CALL   CLRSCN 
TOT=BP+HP-KJP 
IF(T0T.NE.1.)THEN 
WRITE(*,10) 
10  FORMAT(/////'  THE   EXPONENTS  MUST  SUM  TO   1.0'/'      RE-ENTER  TH 

*EM') 

GO  TO   11 
ENDIF 

CALL  CLRSCN 
WRITE(*.4004) 
4004     F0RMAT(//'  DECAY   RATE   DATA  ENTRY   AREA'. 

C 

WRITE(*,20) 

20  FORMAT(///'  ENTER  A  DECAY   RATE   FOR  THE   RIVER  BACKGROUND:    '\) 
READ(*,*)RBK 

C 
C 

WRITE(*,4000) 

4000  FORMAT{//'        ENTER  A  DECAY   RATE   AT   EACH  TRANSECT    '/) 
DO  4001    1=1, NTR 

WRITE(*, 4002)1 
4002     FORMATC  TRANSECT    ',12,'    :    '\) 

READ(*,*)RKS(I) 

4001  CONTINUE 
C 

c 
c 
c 

WRITE(*,21) 

21  FORMAT(/'  AT  WHAT   TEMPERATURE    IS  THIS   RATE   KNOWN?    IN  C   :    '\) 
READ(*,*)TEMP 

WRITE(*,22) 

22  FORMAT(/'    WHAT  IS  THE  RIVER  TEMPERATURE?  IN  C  :   '\) 
READ(*,*)TMP 

CALL  CLRSCN 
C 

WRITE(*,50) 

50  FORMATC/////'    DO  YOU  WISH  TO  CONSIDER  AMMONIA?  YES=1  N0=0   '\ 

*) 

READ(*,*)IAM 
IF(IAM.EQ.1)THEN 
WRITE(*,51) 

51  FORMAT(/'   ENTER  PH   '\) 
READ(*,*)PH 

ELSE 
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PH=7.0 
ENDIF 

C 

0PEN(7,FILE='PPINCAL.DAT',STATUS='NEW') 

WRITE(7,1)TITLE 

WRITE(7,111)QRS,BP,HP,UP.THETA,TEMP,RBK 

111  FORMAT(2X,F9.2.4(2X.F5.3),2X.F6.2,2X,F10.7) 
WRITE(7,112)QRUP,QEFL.CEFL,CBKG,TMP 

112  F0RMAT(2X,F9.2,2X,F8.2,2X,F8.3,2X,F6.3,2X,F6.2] 
WRITE(7,*)NTR,NYZ,QCP,Y0UT 

DO  40  1=1, NTR 
WRITE(7,*)X(I),BS(I),ZAV(I).VAV(I) 

40  CONTINUE 

DO  41  1=1, NTR 
WRITE(7,*)RKS(I) 

41  CONTINUE 
WRITE(7,*)IAM,PH 
WRITE(7,42)CHAR(26) 

42  FORMAT(A) 
CL0SE(7) 
CALL  CLRSCN 
RETURN 

END 
C 
C 

c 

SUBROUTINE  CLRSCN 
WRITE(*,101)CHAR(27),'[2J' 
101   F0RMAT(1X,A,A\) 
RETURN 
END 
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c  * 

C  PROGRAM  COMPLOT. FOR  COMPARES  THE  OBSERVED  CONCENTRATIONS  * 

C  WITH  THOSE  CALCULATED  BY  MIXCALBN.FOR  TO  ADJUST  THE   VALUE  * 

C  OF  BETA  * 

C  WRIHEN  BY   R.    JARVIS  * 

C  * 

C  GORE  &  STORRIE     1986 

C********ir***************************************************^r********** 

COMMON/A/  C0(10.30),YST(10.30).X(20).CP(10,30),BETA(10).SKD(10) 

COMMON/B/  QQY(10,30).CMAX(10),YY(32).CC(32).CD(32).CFACT(10) 

DIMENSION  P(4),PP(4),CSAV(10).CSAV2(10),C2SAV(10).C2SAV2(10) 

CHARACTER*70  TITLEl ,TITLE3 

CHARACTER*80  TITLE 
C 

TITLEU 'COMPARISON  OF  OBSERVED  CONCENTRATIONS  TO  CALIBRATION  CONCE 
♦NTRAT IONS' 

DO  922   1=1,10 

C2SAV(I)=0. 

C2SAV2(I)=0. 

CSAV2(I)=0. 

CSAV(I)=0. 

CMAX(I)=0.0 
922       CONTINUE 

CALL   CLRSCN 

WRITE(*,801) 
801        F0RMAT(/////'      HOW  DO  YOU  WANT   LATERAL  UNITS:      1=  STREAM  TUBE   COOR 
*D.S' 

*/'  2=   REAL  DISTANCE    IN  METERS') 

READ(*,*)IUN 

IF(IUN.EQ.1)TITLE3='LATERAL  COORDINATES    IN  STREAM  TUBE   UNITS' 
IF(IUN.EQ.2)TITLE3='LATERAL   COORDINATES    IN   REAL  UNITS' 
C 

0PEN(5,FILE='PLCALPC' ,STATUS= 'OLD ' ) 

READ(5,*)ICAL 

READ(5,610)TITLE 
610        FORMAT(A) 

READ(5,*)NTR 

READ(5.*)NYZ 

NYZP1=NY2+1 

DO   79    1=1, NTR 

READ(5.*)ID,XX 
79  CONTINUE 

DO   78   1  =  1, NTR 

READ(5,*)ID,SKD(I) 
78  CONTINUE 

DO   77    1=1, NTR 

READ(5.*)ID,BETA(I) 
77  CONTINUE 

DO  66   1  =  1, NTR 

DO  88  J=1,NYZP1 

READ(5.*)ID,QQY(I,J).C0(I,J) 

IF(CO(I,J).GT.CMAX(I))CMAX(I)=CO(I.J) 
88         CONTINUE 

DO   1920  J=1,NYZ 

CSAV(I)=CSAV(I)  +  .l*(C0(I.J-Hl)+C0(I.J))/2. 
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C2SAV(I)=C2SAV(I)-H(.l*FL0AT(J))**2*(C0(I,J+l)^0(I,J))/2. 
1920  CONTINUE 

C2SAV(I)=C2SAV(I)/CSAV(I) 
66    CONTINUE 

CL0SE(5) 

0PEN(5,FILE='PL0T0UT.DAT'.STATUS='0LD') 

READ(5,610)TITD 

READ(5,*)NTRR 

READ(5,*)NYZZ 

DO  33  1=1, NTR 

DO  44  J=1,NYZP1 

READ(5.*)ID,YST(I,J),D,D,CP(I,J) 

IF(CP(I,J).GT.CMAX(I))CMAX(I)=CP(I,J) 
44    CONTINUE 

DO  1902  J=1,NYZ 

CSAV2(I)=CSAV2(I)-t-.l*(CP(I,J+l)+CP(I,J))/2. 

C2SAV2(I)=C2SAV2(I)t-(.l*FL0AT(J))**2*(CP(I,J+l)+CP(I,J))/2. 
1902  CONTINUE 

C2SAV2( I ) =C2SAV2( I ) /CSAV2 ( I ) 
33    CONTINUE 

CL0SE(5) 
C 

WRITE(*,444) 
444   FORMATC    PLOT  LOCATION:  SCREEN=1  PL0TTER=2  PRINTER=3   '\) 

READ(*,*)IPL 

IF(IPL.EQ.3)CALL  PL0TS(0,0,11) 

IF(IPL.EQ.3)CALL  WIND0W(0. ,0. ,20. ,20. ) 

IF(IPL.EQ.1)CALL  PL0TS(0,0,99) 

IF( IPL.EQ.2)CALL  PL0TS(0,9600,80) 

IF(IPL.EQ.2)CALL  FACT0R(.16) 

IF(IPL.EQ.1)CALL  FACT0R(.13) 

IF(IPL.EQ.3)CALL  FACT0R(.15) 

CALL  DUPLX 


CALL  SYMB0L(5.,42.,1.,TITLE,0.,80) 

CALL  SYMB0L(5.,40.,.75,TITLE1,0.,70) 

CALL  SYMBOL(5.,38.,.75,TITLE3,0.,50) 

CALL  SYMB0L(5., 36. ,.5, 'OBSERVED  CONCENTRATIONS' ,0. ,23) 

CALL  SYMB0L(5. ,35. ,.5. 'PREDICTED  CONCENTRATIONS ' ,0. ,24) 

P(1)=0. 

P(2)=0. 

PP(1)=0. 

PP(2)=10. 

CALL  SCALE(P,2.,2,1) 

CALL  SCALE(PP,5..2,1) 

CALL  STDASH(.5,.2) 

CALL  PL0T(18.,36.5,-3) 

CALL  STLINE(-1,.15,0.) 

CALL  C0L0R(2,IERR) 

CALL  LINE(PP,P.2, 1,1.1) 

CALL  PL0T(0.,-l.,-3) 

CALL  STLINE(+1..15,0.) 

CALL  C0L0R(1,IERR) 

CALL  LINE(PP,P,2, 1,1,1) 

CALL  STLINE(+1,.15,0.) 

CALL  PL0T(-18. .-35.5,-3) 
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CALL  COLOR(O.IERR) 
C 

DO  300  1=1, NTR 

CSCALE=CMAX(I)/8. 

DO  400  J=1,NYZP1 

CC(J)=C0(I,J) 

CD(J)=CP(I,J) 

IF(IUN.EQ.2)YY(J)=YST(I.J) 

IF(IUN.EQ.1)YY(J)=FL0AT(J-1)/10. 
400   CONTINUE 

CC(NYZPl-t-l)=0. 

CC(NYZP1+2)=CSCALE 

CD(NYZP1+1)=0. 

CD(NYZP1+2)=CSCALE 

CALL  SCALE(YY,10.,NYZP1,1) 

IF(I.EQ.1.0R.I.EQ.5)X0R=2.5 

IF(I.EQ.2.0R.I.EQ.6)X0R=15.0 

IF(I.EQ.3.0R.I.EQ.7)X0R=27.5 

IF(I.EQ.4.0R.I.EQ.8)X0R=40. 

IF(I.G£.1.AND.I.LE.4)Y0R=19.5 

IF(I.GE.5.AND.I.LE.8)Y0R=2.5 

CALL  PL0T(X0R,Y0R.-3) 

CALL  STAXIS(.3,.5..10,.3,1) 

IF(IUN.EQ.l)  CALL  AXIS(0. ,0. , 'QY/QT ' ,-5,10. ,0. , YY(NYZPl-fl) ,YY(NYZP 
*U2)) 

IF(IUN.EQ.2)CALL  AX  IS (0.  ,0. ,' LATERAL  DISTANCE ', -16, 10. ,0. .YY(NYZP1 
*+l),YY(NYZPl+2)) 

CALL  AXIS(0.,0. .'CONCENTRATION' ,13,8. ,90. ,0.,CSCALE) 

CALL  STDASH(.5,.2) 

CALL  STLINE(+1,.2,0.) 

CALL  C0L0R(1,IERR) 

CALL  LINE(YY,CC,NYZP1, 1,1,1) 

CALL  STLINE(-1,.2,0.) 

CALL  C0L0R(2,IERR) 

CALL  LINE(YY,CD,NYZP1,1,1,2) 

CALL  COLOR(0,IERR) 

CALL  SYMB0L(2. 0,13., .5, 'TRANSECT',  0.  ,8) 

SNUM=FL0AT(I) 

CALL  NUMBER(6.0,13.,.5,SNUM,0.,-1) 

CALL  SYMB0L(0.0.12.,.35,'BETA=',0.,5) 

SNUM=BETA(I) 

CALL  NUMBER(2.5.12.,.35.SNUM,0.,5) 

CALL  SYMB0L(5.0,12.0,.35,'Kd=',0.,3) 

SNUM=SKD(I) 

CALL  NUMBER(6.25,12.0,.35,SNUM,0.,6) 

CALL  SYMB0L(4. 5, 11., .35, 'OBSERVED   PREDICTED ' ,0. ,20) 

CALL  SYMB0L(0.,10.5,.35,'MEAN  CONC. ' ,0. ,10) 

CALL  SYMBOL(0., 10. ,.35, 'SPREAD    ',0.,10) 

SNUM=CSAV(I) 

CALL  NUMBER(9.5,10.5,.35,SNUM,0.,3) 

SNUM=CSAV2(I) 

CALL  NUMBER(4.5,10.5..35,SNUM,0.,3) 

SNUM=C2SAV(I) 

CALL  NUMBER(9.5,10.O,.35,SNUM,0.,3) 

SNUM=C2SAV2(I) 

CALL  NUMBER(4.5,10.,.35,SNUM,0.,3) 

CALL  PL0T(-X0R,-Y0R,-3) 
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300   CONTINUE 

CALL  PL0T(0.,0.,999) 

STOP 

END 

C 
C 

c 

SUBROUTINE  CLRSCN 
WRITE(*,101)CHAR(27),'[2J' 
101   F0RMAT(1X,A,A\) 
RETURN 
END 
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C  PROGRAM  MIXAPPLN. FOR   IS  THE   PC  VERSION  OF  MIXAPPLN  * 

C  BY   ROB  JARVIS  * 
C 

C  PROGRAM  NAME:      MIXAPPLN     *  *  STREAMTUBE  MODEL   FOR   PIPE  OUTFALL  * 

C  THIS   PROGRAM   IS  SET  UP  FOR  CONSERVATIVE,    NONCONSERVATIVE  WITH  ♦ 

C  FIRST  ORDER  DECAY(VI2. .RESIDUAL  CHLORINE,    PHENOL,    RADIONUCLIDES,    * 

C  INDICATOR  BACTERIA),   AND  UN-IONIZED  AMMONIA  CONSTITUENTS.  * 

C  THIS   PROGRAM   INCLUDES  OPTIONS   FOR  DESIGN    :   QRIVER,QEFFL,TEMP  &  PH  * 

C  PROGRAM  DEVELOPED  BY  T.    P.    H.   GOWDA,   WATER  RESOURCES  BRANCH  * 

C  DATE:   JUNE   1980.  * 

C  * 

C  GORE  &  STORRIE     1986  * 

IMPLICIT  REAL*8   (A-H,0-Z) 

COMMON/A/  C(10,50),CUI(10,50),PH(4) 

COMMON/B/  ARAY1(20,15),ARAY2(20,4),ARAY3(20,4),ARAY4(20,4) 

COMMON/D/  X(10),XX(10),PBK,TBK,PXWC,TXWC,A3,RKS(10),QY(50    ,QRS. 

1THETA,BPWR,HPWR,UPWR,ORTO,QT,QRUP(6),QEFL(6),CTDP,PHDG,VOL(10). 

2B(10),H(10),U(10),BSUM(10),T0T(10),TMP(6),RQ(11),TEMPS,PAX1,PAX2. 

3RF(10),BS(10),HS(10),US(10),BETA(10),BW(10),HW(10).UW(10),RKT(10). 

4PHI(10).R,QCP,PQX.EY(10),XCRIT,PCRIT,XEK,XMZ,PMZ,TMZ,RBK,RCRT   REK, 

5RXS,RKAV,AWCP,RFWC.XWCP,RBT,RBKG,XSCE,XSCEA,PW1,PW2.QCR,PHWC.DELQ. 
6XL,XEST  .        .y      .  .        w, 

CHARACTER*80  TITLE,PARAM1 ,PARAM2 

CHARACTER*20  FILIN,FILOUT,POLLU 

CALL  CLRSCN 
C 

CALL  SETUP(JRUN) 
C 

0PEN(1,FILE='APINCAL.DAT',STATUS='0LD') 
0PEN( 5 , FI LE= ' OUTAPP .DAT ' , STATUS= ' NEW ' ) 
0PEN(6,FILE='PLAPPLN',STATUS='NEW') 
37  WRITE(5,2) 

2  FORMAT(/'      ENTER  TITLE  OF  STUDY') 
READ(1,3)    TITLE 

3  FORMAT(A) 
WRITE(6,*)JRUN 
WRITE(6,3)TITLE 
WRITE(5,3)TITLE 
WRITE(5,1992) 

1992     FORMAT(/'      ENTER  POLLUTANT  NAME') 
READ(1,3)P0LLU 
WRITE(6,3)P0LLU 
WRITE(5,3)POLLU 
WRITE(5,4) 

4  FORMATC      ENTER  QRS,BPWR,HPWR. UPWR, TEMPS. NTR') 

READ(1,*)   QRS,BPWR,HPWR, UPWR, TEMPS, NTR 

WRITE (5, 901 )QRS,BPWR,HPWR, UPWR, TEMPS, NTR 

F0RMAT(3X,F7.1,4X,F5.1,4X,F5.2,4X,F5.2.F5.1.2X.I2) 
WRITE(6,*)NTR  '      '      > 

WRITE(5.6) 

FORMATC      ENTER   NTR  VALUES  OF   X,BS,HS.US') 

DO    15    1=1, NTR 

READ(1,*)    X(I),BS(I).HS(I),US(I) 

WRITE(5,902)X(I),BS(I),HS(I),US(I) 


901 
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15    WRITE(6.902)X(I).BS(I).HS(I),US(I) 

902  F0RMAT(3X,F7.1,3X.F7.1,3X.F6.2.3X,F6.3) 
WRITE(5,8) 

8     FORMATC   ENTER  NTR  VALUES  OF  BETA') 
READ(1.*)  (BETA(I),I=1,NTR) 
DO  904  1=1, NTR 
WRITE(5,903)BETA(I) 

903  F0RMAT(3X,F8.5) 

904  CONTINUE 
WRITE(5,80) 

80    FORMATC   ENTER  MQ  &  QRUP  VALUES' ) 
READ(1,*)  MQ.(QRUP(J),J=1,MQ) 
DO  906  1=1, MQ 
WRITE(5, 907)1, QRUP(I) 

907  FORMATC    QRUP( ' ,  14. '  )=  '  ,F9.2) 
906   CONTINUE 

WRITE(6,*)MQ,(QRUP(I),I=1,MQ) 
WRITE(5.82) 
82    FORMATC   ENTER  MT  &  TEMP  VALUES' ) 
READ(1,*)MT,(TMP(L),L=1,MT) 
DO  908  L=1,MT 
WRITE(5,909)L,TMP(L) 

909  FORMATC    TMP( '  ,13.  '  )  =  ' .F6. 1) 

908  CONTINUE 
WRITE(6.*)MT.(TMP(I).I=1.MT) 
WRITE(5,86) 

86    FORMATC   ENTER  MF  &  QELF  VALUES' ) 
READ(1,*)  MF,(QEFL(L),L=1,MF) 
DO  910  L=1,MF 
WRITE(5,911)L,QEFL(L) 

911  FORMATC  QEFLC.I2.')=',F6.2) 

910  CONTINUE 
WRITE(6,*)MF,(QEFL(I),I=1,MF) 

35    WRITE(5.52) 

52    FORMATC   UN-IONIZED  AMMONIA:  ENTER  1  FOR  YES 
0  FOR  NO') 

READ(1,*)  AMONIA 
WRITE(5,912)AMONIA 

912  FORMATC   AMONIA=' ,F2.0) 
IF(AMONIA.LE. -888.8)  GO  TO  888 
IF(AMONIA.LE.O)GO  TO  70 
WRITE(5,57) 

57    FORMATC   ENTER  MPH  AND  PH  VALUES') 
READ(1.*)MPH,(PH(JPH),JPH=1,MPH) 
DO  914  JJ=1.MPH 
WRITE(5,913)JJ,PH(JJ) 

913  FORMATC   PH( ' ,  12. '  )=  '  ,F7.3) 

914  CONTINUE 
WRITE(6.*)MPH.(PH(I).I=1,MPH) 

70  CONTINUE 

WRITE(5,56) 
56         FORMATC      ENTER  QCP.CEFL.CBKG.CS.THETA.RBK.XWCP ' ) 

READ(1,*)QCP.CEFL.CBKG.CS,THETA.RBK,XWCP 

WRITE(5.915)QCP,CEFL,CBKG,CS,THETA,RBK.XWCP 

WRITE(6.915)QCP.CEFL.CBK6,CS.THETA.RBK.XWCP 

915  F0RMAT(2X,F5.2.3X,F9.2.3X,F6.2,3X.F7.2.3X.F5.2.3X.F5.2.F8.1) 
WRITE(5,55) 
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55 


917 
916 


92 

C 
C 


34 
88 


54 


FORMATC   ENTER  NTR  VALUES  OF  RKS 

READ(1,*)  (RKS(I),I=1,NTR) 

DO  916  1  =  1, NTR 

WRITE(5, 917)1, RKS(I) 

F0RMAT(2X,'RKS(',I2,')=',F8.6) 

CONTINUE 

WRITE(6.*)(RKS(I),I=1,NTR) 


10 

24 

C 
C 

c 


CALCULATE  FLOW  &  TEMP'R  SCALE-UP  PARAHETERS 

DELTA=0.0001 
AL0G(1/DELTA)=(4. 0*2. 3026)  FOR  CALCN.  OF  BGX 

NQR=11 

NW=5 

DO  10  L=1,NQR 

RQ(L)=(L-1)/10.0 

WRITE(5,24)  TITLE 

F0RMAT(1H1/////6X, 'PREDICTIONS  OF  RUNS  FOR  MANAGEMENT  OPTIONS'/A) 


BEGIN  COMPUTATIONS  FOR  THE  INPUT  OPTIONS. 

DO  20  JF=1,MF 

DO  20  J0=1,MQ 

QT=0RUP(JQ)+OEFL(JF) 

DELQ=QT/(NQR-1) 

QRTO=QT/QRS 

CA=(CEFL*QEFL(JF)/QT) 

DO  20  JT=1,MT 

CTDP=THETA**(TMP(JT) -TEMPS) 

RBT=RBK*CTDP 

IF(AM0NIA.LE.0)MPH=1 

DO  20  JPH=1,MPH 

IRUN=IRUN+1 

IR=IRUN 

IF(AM0NIA.LE.0)PH(JPH)=7.0 

ARAY1(IR,1)=IRUN 

ARAY1(IR,2)=QEFL(JF) 

ARAY1(IR,3)-QRUP(JQ) 

ARAY1(IR,4)=TMP(JT) 

ARAY1(IR,5)=PH(JPH) 

ARAY1(IR,12)=CS 

ARAY1(IR,15)=CA 

WRITEf5,92)IRUN 


FORMAT(/' 


RUN  NO. 


.14) 


CALCULATE  AMMONIA  IONIZATION  PARAMETER 

IF(AMONIA.LE.O)  GO  TO  34 

PICA=0. 09018+2729. 92/(TMP(JT)+273. 2) 
PF=PKA-PH(JPH) 
PCTU=1./(1.+10.**PF) 
CSTTL=CS/PCTU 

WRITE(5.88)QEFL(JF),QRUP(JQ),TMP(JT).PH(JPH),CEFL.CS 
FORMAT(2X.'QEFL=',F8.3,2X,'ORUP=',F10.3,2X,'TEMPR=',F4.1  2X  'PH= 
1F4.1./2X,'CEFL='.F8.2.2X.'CS=',F5.2) 
IF(AMONIA.GE.l)WRITE(5.54)CSTTL 
F0RMAT(5X.'   CRITERION  FOR  TOTAL  AhWONIA,  CS=',F6.3) 
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WRITE(5,90)  (RQ(L),L=1.NW) 
90    F0RMAT(/6X,'X',6X,'EY',7X,5(F3.1,6X),//) 

DO  12  1=1, NTR 

BSUM(I)=0. 

VOL(I)=0.0 
12    TOT(I)=0. 
C 
C     BEGIN  COMPUTATIONS  AT  TRANSECT, I.  ' 

c 

DO  14  1  =  1, NTR 
C 

C        CALCULATE   B,H,U  FOR  DESIGN   FLOW=QT,    FROM  LEOPOLD-HADDOCK  EQNS. 
C 

B(I)=BS(I)*QRTO**BPWR 

H(I)=HS(I)*QRTO**HPWR 

U(I)=US(I)*QRTO**UPWR 
C 

C        CALCULATE  WEIGHTED  MEAN   VALUES  BW,HW,UW  FROM  OUTFALL  TO  TRANSECT(I) 
C 

IF(I.GE.2)   GO  TO  60 

XX(1)=X(1) 

BW(1)=B(1) 

HW(1)=H(1) 

UW(1)=U(1) 

BSUM(1)=B(1)*XX(1) 

V0L(1)=XX(1)*B(1)*H(1) 

T0T(1)=XX(1)/U(1) 

GO  TO  62 
60         11=1-1 

XX(I)=X(I)-X(I1) 

BSUM(I)=BSUM(I1)+0.5*XX(I)*(B(I1)+6(I)) 

V0L(I)=V0L(I1)+0.25*XX(I)*(B(I1)+B(I))*(H(I1)+H(I)) 

T0T(I)=T0T(I1)+XX(I)/U(I) 

BW(I)=BSUM(I)/X(I) 

HW(I)=V0L(I)/(X(I)*BW(I)) 

UW(I)=QT/(BW(I)*HW(I)) 
C 

C   CALCULATE  PRODUCT  FUNCTION  FOR  DECAY,  RF(I),  &  DISPERSION  FACTOR 
C 
62    RKT(I)=CTDP*RKS(I) 

PBK=RBT*TOT(I) 

CALL  PDET(PBK,TBK) 

CBKX=CBKG*TBK 

A3=(RKT{I)*XX(I))/U(I) 

R=DEXP(-A3) 

IF(I.GE.2)   GO  TO  64 

RF(1)=R 

GO  TO  66 
64         RF(I)=RF(I1)*R 
66         EY(I)=BETA(I)*B(I)*U(I) 

PHI(I)=BETA(I)*X(I)/BW(I) 

PHDG=4.0*PHI(I) 

CMAX=0.5*CA*RF(I)/DSQRT(3.1416*PHI(I)) 
C 

C     LATERAL  CONC.   DISTR'N  AT  TRANSECT,!. 
C 

DO   16  K=1,NW 
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QY(K)=(K-1)*0ELQ 

IF  (QY(K).GT.QT)  QY(K)=QT 

PAX1=(QY(K)-QCP)/QT 

PAX2=(QY(K)-K3CP)/QT 

CALL  SUMSRS(PAX1,PAX2,PHDG,SUHT) 

C(I,K)=CMAX*SUMT+CBKX 

16  CONTINUE 

WRITE (5,25)  X(I),EY(I).(C(I,K),K=1,NW) 

25  F0RMAT(2X,F8.1,1X,F8.3.11F9.3) 
14    CONTINUE 

C 

C     UN-IONIZED  AMMONIA  CONCENTRATION  DISTRN. 

C 

IF(AMONIA.LE.O)GO  TO   191 

WRITE(5,26) 

26  FORMAT (2X, 'TOXIC  AMMONIA') 
DO  19  1=1, NTR 

DO  17  K=1,NW 
CUI(I,K)=C(I,K)*PCTU 

17  CONTINUE 

WRITE (5,25)  X(I),EY(I),(CUI(I.K),K=1,NW) 
19    CONTINUE 
191   CONTINUE 
C 

C     END  OF  COMPUTATIONS  AT  TRANSECT, I. 
C 

IF(AMONIA.LE.O)CSL=CS 

IF(AMONIA.GE.l)CSL=CSnL 

c 

C     COMPUTE  BACKGROUND  AVG.  CONC.  AT  D/S  WPCP. 
C 

RKAV=-DLOG(RF(NTR))/TOT(NTR) 

AWCP=RKAV*XWCP/UW(NTR) 

CALL  PDET(AWCP,RFWC) 

CBA=CA*RFWC 
76    PXWC=RBT*XWCP/UW(NTR) 

CALL  PDET(PXWC,TXWC) 

CBB=CBKG*TXWC 

CAWP=CBA+CBB 
C     COMPUTE  BANK  CONC.  AT  D/S  WPCP. 

PW1=0. 

PW2=0. 

PHWC=4.0*BETA(NTR)*XWCP/BW(NTR) 

CALL  SUMSRS(PW1,PW2,PHWC.SUMWC) 

CWCP=CBA*SUMWC/DSQRT(3.1416*PHWC)+CBB 
C 

C  COMPUTATIONS  FOR  MIXING  ZONE  PARAMETERS 
C 

XMZ=B(NTR)/BETA(NTR) 

PMZ=RKAV*XMZ/UW(NTR) 

CALL   PDET(PMZ.TMZ) 

CMZ=CA*TMZ 

c 

C  CALCULATE  XSCE 

C 

IF(CWCP.GT.CSL)  GOTO  220 

CALL  PARSPR(C,CSL,X,NTR,KXS.XEST,CXS) 
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RXS=-DLOG(RF(KXS))/TOT(KXS) 

CALL  SPREAD(CSL,CA.CBKG.RBT,CXS,XEST,XSCE.KXS,BW,RXS.UW.BETA) 

GO  TO  222 
220   XSCE=-999.0 

ITRN=0 
222   WRITE(5,50)XSCE.XMZ,CMZ,XWCP,CWCP,CAWP 

WRITE(6,*)XMZ 
50    F0RMAT(/5X.'XS  (WITH  CE)= ' ,F10. 1,/5X, 

♦'MIXING  ZONE  LENGTH='.F9.1,5X,'C0NC=',F8.2/5X,'DIST.  TO  D/S  WPCP=' 
*,F8.1,5X, 'SHORE  CONC.  AT  D/S  WPCP=' .F6.2/8X, ' AVG.  CONC.  AT  D/S  WPC 
*P=  ',F6.2) 

ARAY1(IR,6)=CAWP 

ARAY1(IR,7)=XSCE 
C 

C  CRITICAL  POINT  RESULTS  FOR  QRATI0=0.1  ,  0.2  &  0.3 
C 

WRITE(5,94) 
94    F0RMAT(/9X.'     CRITICAL  POINT  METHOD  RESULTS  '/9X, 'QY/QT' ,6X, 
*'XL',8X,'CL',9X,'CEA',6X,'XSCEA') 

DO  20  K=2,5 

KK=K-1 

CCRIT=-1.0D+10 
C 

C   SEARCH  FOR  TRANSECT  NEAR  WHICH  CRIT.  CONC.  OCCURS  TO  FIND  MOVING 
C     AVG.  VALUES  FOR  CRITICAL  POINT  COMPUTATIONS 
C 

DO  18  1=1, NTR 

IF(CCRIT.GE.C(I,K))  GO  TO  18 

ICR  =  I 

CCRIT=C(I,K) 
18    CONTINUE 

XCRIT=X(ICR) 

PCRIT=PHI(ICR) 

RCRT=-DLOG(RF(ICR))/TOT(ICR) 
C 

C   CALCULATE  XL,CL  &  CEA  BY  CRIT.  POINT  METHOD  USING  MOVING  AVG.  VALUES 
C 
40    QCR=4.0*RCRT*XCRIT/(PCRIT*UW(ICR))*(QY(K)/QT)**2 

IF(RF(ICR).EQ.1.0)G0  TO  9007 

XL=0.25*UW(ICR)*(-1.0+OSQRT(1.+QCR))/RCRT 

GO  TO  9008 

9007  XL=.5*XCRIT*QY(K)**2/(PCRIT*QT**2) 

9008  XLDIF=DABS(XL-XCRIT) 
XPCT=100.*XLDIF/XL 
IF(XPCT.LE.5.)  GO  TO  42 
XCRIT=XL 

PCRIT=BETA(ICR)*XL/BW(ICR) 
GO  TO  40 

42    PQX=(RCRT*XL/UW(ICR))+(0.25/PCRIT*(QY(K)/QT)**2) 

IF(RF(ICR).EQ.1.0)GO  TO  9010 
CL=(CA/DSQRT(3.1416*PCRIT))*DEXP(-PQX)+CBKG*DEXP(-RBT*XL/UW(ICR)) 

GO  TO  9011 

9010  CL=CA*QT/(QY(K)*2.07) 

9011  CEA=CEFL*CSL/CL 
C 

C     COMPUTE  XSCEA 
C 
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IF(CWCP.GT.CL)GO  TO  114 

CALL  PARSPR(C,CL,X,NTR.IEK,XEK,CXK) 

REK=-DLOG(RF(IEK))/TOT(IEK) 

CALL  SPREAD(CL,CA.CBKG,RBT.CXK,XEK,XSCEA,IEK.BW,REK,UW,BETA) 

GO  TO  116 
114   XSCEA=-999.0 

ITRA=0 
116   WRITE(5,96)RQ(K),XL,CL,CEA,XSCEA 

WRITE(6,996)XL,CL 
996   F0RMAT(2X,F8.1.2X,F8.3) 

96    F0RMAT(9X,F4.2,2X,F8.1,2X,F8.3,4X,F8.2,3X,F8.1) 
C 

C     STORE  OUTPUT  IN  ARRAYS 
C 

GO  TO  (45, 46, 47, 477), KK 

45  ARAY1(IR,8)=RQ(2) 
ARAY1(IR,9)=CEA 
ARAY1(IR,10)=XSCEA 
GO  TO  20 

46  ARAY2(IR,1)=RQ(3) 
ARAY2(IR,2)=CEA 
ARAY2(IR,3)=XSCEA 
GO  TO  20 

47  ARAY3(IR,1)=RQ(4) 
ARAY3(IR,2)=CEA 
ARAY3(IR,3)=XSCEA 
GO  TO  20 

477   ARAY4(IR,1)=RQ(5) 

ARAY4(IR,2)=CEA 

ARAY4(IR,3)=XSCEA 
20    CONTINUE 
C      WRITE(6,*)MZL 
888   CONTINUE 

WRITE(5,138) 
138   F0RKAT(1H,//) 

WRITE(5,78)TITLE 
78    FORMAT (1H1//6X, 'SUMMARY  OF  RUNS  FOR  MANAGEMENT  OPTIONS'/A) 

WRITE(5,100) 
100       F0RMAT(//1X,'RUN  #' ,2X, 'QEFL' ,5X, 'QRUP' ,3X, 'TEMP' ,3X, 'PH' ,4X, 
* ' CAWP ' , 3X , ' XSCE ' , 4X , ' QY /QT ' , 4X , ' CEA ' , 4X , ' XSCE A ' 
*,5X,'CBKG',4X,'CSIJC',5X,'CBI0T',4X,'CDRNK',4X,'CDILN',//) 

DO  48  N=1,IR 

WRITE(5,102)(ARAY1(N,J),J=1,15) 
102       F0RMAT(1X,F4.0,1X,F7.3,1X,F7.1,2X,2F5.1.1X,F7.3, 

*1X,F7.1,2X,F4.2,2X,F8.2,1X,F8.1,2X,F6.2,3(2X,F7.1),3X,F7.3) 

WRITE(5,104)(ARAY2(N,J),J=1.3) 
104        F0RMAT(51X,F4.2,2X,F8.2,1X,F8.1) 

WRITE(5,104)(ARAY3(N,J),J=1,3) 

48  WRITE(5,104)(ARAY4(N,J),J=1,3) 
WRITE(5,9)    PARAMl 

9  F0RMAT(//1X,A) 

IF(NPARM.GE.7)    WRITE(5,3)    PARAM2 

WRITE(5.138) 

CLOSE(l) 
999       CLOSE (1) 

WRITE(5.5558)CHAR(26) 
5558     FORMAT(A) 
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WRITE(5,5558)CHAR(26) 

CLOSE (5) 

CLOSE (6) 

WRITE(*.1092) 
1092     FORMAT(//'     THE  OUTPUT  FILE  FROM  PROGRAM  MIXAPPLN   IS'/'      "OUTAPP.D 
*AT" ' ) 

STOP 

END 
C 

C  SUBROUTINE  FOR  SUMMATION  OF  EXPONENTIAL  SERIES  TERMS. 
C 

SUBROUTINE  SUMSRS(QAX1,QAX2,PFDR,SUM) 

IMPLICIT  REAL*8  (A-H,0-Z) 

BGX=2.3026*PFDR 

SBG=DSQRT(BGX) 
C 

C  DETERMINE  NO.  OF  IMAGES  REQUIRED 
C 

ANl=(0.5*QAXl-SBG)-0.5 

AN2=(0.5*QAXl+SBG)+0.5 

AN3=-AN2 

AN4=-AN1 

NM1=IDINT(AN1) 

NM2=IDINT(AN2) 

NM3=IDINT(AN3) 

NM4=IDINT(AN4) 

NN1=1+NM2+IABS(NM1) 

NN2=1+IABS(NM3)+IABS(NM4) 

IF(NN1.GE.NN2)NN=NN1+1 

IF(NN1.LT.NN2)NN=NN2+1 
C 

C   COMPUTE  SUM  OF  EXPONENTIAL  SERIES  TERMS 
C 

SUM=0. 

DO  32  J=1,NN 

N=J-1 

P1=(QAX1-2.*N)**2/PFDR 

P2=(QAX2+2.*N)**2/PFDR 

CALL  PDET(P1,T1) 

CALL  PDET(P2,T2) 

IF(N.LE.O)  GO  TO  30 

P3=(QAX1+2.*N)**2/PFDR 

P4=(QAX2-2.*N)**2/PFDR 

CALL  PDET(P3,T3) 

CALL  PDET(P4,T4) 

GO  TO  32 
30  T3=0. 

T4=0. 
32  SUM=SUM+T1+T2+T3+T4 

RETURN 

END 
C 

C  COMPUTE  EXPONENTIAL  TERMS:  SET  (P.LE.40.0)  TO  AVOID  ERROR  208,  SO 
C    THAT  EXP(-P)=4.3E-18. 
C 

SUBROUTINE  PDET(P,T) 
IMPLICIT  REAL*8  (A-H,0-Z) 
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IF(P.GE.40.)G0  TO  10 

T=DEXP(-P) 

60  TO  12 
10  T=0.0 
12  CONTINUE 

RETURN 

END 
C 

C  COMPUTE  PARAMETERS  IN  THE  SUBROUTINE  'SPREAD'. 
C 

SUBROUTINE  PARSPR(C,CL,X,NTR, lEK.XEK.CXK) 

IMPLICIT  REAL*8  (A-H,0-Z) 

DIMENSION  C(1.1),X(1) 

XEK=0. 

IF(C(NTR,1).GE.CL)  GO  TO  75 

DO  72  I=2,NTR 

11=1-1 

IF(C(1,1).LE.CL)G0  TO  73 

IF(C(II.1).GT.CL.AND.C(I,1).LE.CL)  GO  TO  73 

72  CONTINUE 

73  IEK=II 
GO  TO  79 

75   IEK=NTR 
79  XEK=X(IEK) 

CXK=C(IEK,1) 

RETURN 

END 
C 

C   COMPUTATION  OF  MAX.  LONGL.  SPREAD  XS  ALONG  OUTFALL  BANK  WHERE 
C   (C(XS,0)-CS)=5  PERCENT(ABSOLUTE). 
C 

SUBROUTINE  SPREAD(CST,CAV,CB6,RB,CXY,XEY,XST,M,BI ,RS,US,BTA) 

IMPLICIT  REAL*8  (A-H.O-Z) 

DIMENSION  BI(1),BTA(1),US(1) 

CXX=CXY 

XS=XEY 

IT=0 
18   IT=ITfl 

DIFF=(CXX-CST) 

RCAB=OABS(DIFF/CST) 

PRCNT=100.*RCAB 

IF(PRCNT.LE.5.)G0  TO  26 

IF(IT.GT.30)  GO  TO  27 

FRX=RCAB/(1.0+PS) 

IF(DIFF.LE.O.O)  XS-XS*( 1-FRX) 

IF(DIFF.GT.O.O)  XS=XS*(  1-t-FRX) 

PHB=4.0*BTA(M)*XS/BI(M) 

PS=RS*XS/US(M) 

CALL  PDET(PS.TS) 

CMM=CAV*TS/DSQRT(3. U16*PHB) 

PX1=0. 

PX2=0. 

PG=RB*XS/US(M) 

CALL  PDET(PG.TG) 

CALL  SUMSRS(PX1,PX2.PHB.SUMP) 

CXX=CMH*SUMP+CBG*TG 

60  TO  18 
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27  XS=-888.0 
26  XST=XS 

RETURN 

END 
C 

C     SUBROUTINE  SETUP  READS  A  RAW  INPUT  FILE  AND  ASKS  THE  USER  FOR 
C     ADDITIONAL  INFORMATION  AND  WRITES  ALL  THE  DATA  IN  A  FILE 
C     APINCAL.DAT  TO  BE  READ  BY  THE  MAIN  OF  MIXAPPLN. FOR 
C 

SUBROUTINE  SETUP(IRUN) 

DIMENSION  X(10),BS(10).HS(10),US(10),QRUP(10),QEFL(10),TMP(10) 
*,PH(10),BETA(10),RKS(10),CFACT(10) 

CHARACTER*80  TITLE, CONNAME 

CHARACTER*20  POLLU 

CHARACTER*2  PPP 

0PEN(5,FILE='PINCAL.DAT',STATUS='0LD') 

IRUN=0 

READ(5,1)TITLE 

I  FORMAT(A) 
READ(5,*)QRS,BPWR,HPWR,UPWR,THETA 
READ(5,*)0RUP1,QEFL1,CEFL,CBKG,TEMP 
READ(5,*)NTR,NYZ,QCP 

DO  3  1=1, NTR 
READ(5,*)X(I),BS(I),HS(I),US(I) 

3  CONTINUE 
CL0SE(5) 
CALL  CLRSCN 
WRITE(*.4) 

4  F0RMAT(/'         SUMMARY  OF  INPUT  DATA') 
WRITE(*,5) 

5  FORMATf'         **-***■*-**-*•*-**-**-*•*-***•*-*  1  \ 

WRITE(*,12) 

12  FORMAT(//'         REFERENCE  RIVER  PARAMETERS'//) 
WRITE(*,11)QRS 

II  FORMAT(/'   TOTAL  RIVER  FLOW  BELOW  OUTFALL  AT  TIME  OF  SURVEY  :  ',F 
*10.3) 

WRITE(*,13) 

13  FORMAT(/'   TRANSECT  DISTANCE ',5X,' RIVER  WIDTH' ,5X, 'AVERAGE  DEPTH' 
*,5X, 'AVERAGE  VELOCITY') 

DO  100  1=1, NTR 
WRITE(*,14)X(I),BS(I),HS(I),US(I) 

14  F0RMAT(6X,F10.2,8X,F8.2,12X,F6.2,11X,F6.2) 
100   CONTINUE 

WRITEf*   15) 

15  FORMAT(/'        NOW  YOU  MAY  ENTER  DESIGN  PARAMETERS'/'    STRIKE     [ENTER] 
*  TO  CONTINUE'X) 

READ(*,1)PPP 

CALL  CLRSCN 

WRITE(*,10) 
10    FORMAT(/'   ENTER  STUDY  TITLE  :  '\) 

READ(*,1)C0NNAME 

WRITE(*.1293) 
1293  FORMAT(/'   ENTER  POLLUTANT  NAME  :  '\) 

READ(*,1)P0LLU 
C 

CALL  CLRSCN 

WRITE(*.6) 
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6  FORMAT(/'   ENTER  #  OF  UPSTREAM  FLOW  RATES  (<9)  :  '\) 
READ(*.*)IFR 

WRITE(*,16) 

16  FORMAT(/'   ENTER  THESE  FLOW  RATES') 
DO  101  1  =  1,  IFR 

WRITE(*,17)I 

17  F0RMAT(10X/'    RATE    ',12,':      '\) 
READ(*,*)QRUP(I) 

101  CONTINUE 
C 

CALL  CLRSCN 
WRITE(*,7) 

7  F0RMAT(//'   ENTER  #  OF  EFFLUENT  FLOW  RATES  (<7)  :  '\) 
READ(*,*)IEFL 

WRITE(*,18) 

18  FORMAT(/'   ENTER  THESE  FLOW  RATES') 
DO  102  I=1,IEFL 

WRITE(*.19)I 

19  F0RMAT(10X/'    RATE    ',12,':      '\) 
READ(*.*)QEFL(I) 

102  CONTINUE 
C 

CALL  CLRSCN 
WRITE(*,20) 

20  FORMAT(//'   ENTER  #  OF  RIVER  TEMPERATURES  (<7)  :  '\) 
READ(*,*)ITEM 

WRITE(*,21) 
2-1    FORMAT(/'   ENTER  THESE  TEMPERATURES') 
DO  103  1  =  1, ITEM 
WRITE(*,22)I 

22  F0RMAT(10X,'  TEMP.  ',12,':   '\) 
READ(*,*)TMP(I) 

103  CONTINUE 
C 

CALL   CLRSCN 
WRITE(*,23) 

23  FORMAT(//'      ARE   YOU  CONSIDERING  AMMONIA?   1=YES     0=N0   :    '\) 
READ(*,*)IAM 

IF(IAM.EQ.1)THEN 
WRITE(*,24) 

24  FORMAT(/'   ENTER  #  OF  PH  VALUES  (<5)  :  '\) 
READ(*,*)IPH 

WRITE(*,25) 

25  FORMAT(/'   ENTER  THESE  PH  VALUES') 
DO  104  1=1, IPH 

WRITE(*,26)I 

26  F0RMAT(10X,'  PH  ',12,':   '\) 
READ(*,*)PH(I) 

104  CONTINUE 
ENDIF 

C 

WRITE(*,27) 

27  FORMAT(//'   ENTER  EFFLUENT  CONCENTRATION  :  '\) 
READ(*,*)CEFL 

WRITE(*,28) 

28  FORMAT(//'   ENTER  BACKGROUND  CONCENTRATION  :  '\) 
READ(*,*)CBKG 


I 
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WRITE(*,29) 

29  FORMAT(//'   ENTER  PROVINCIAL  WATER  QUALITY  OBJECTIVE  :  '\) 
READ(*,*)CS 
WRITE(*  30) 

30  FORMAT(//'  ENTER  TEMPERATURE  COEFFICIENT  :  '\) 
READ(*,*)THETA 
WRITE(*,31) 

31  FORMAT(//'   ENTER  DECAY  RATE  OF  BACKGROUND  :  '\) 
READ(*,*)RBK 
WRITE(*,32) 

32  FORMAT(//'   ENTER  THE  TEMPERATURE  THIS  RATE  IS  KNOWN  AT  :  '\) 
READ(*,*)TEMP 
WRITEf*  33) 

33  FORMAT(//'   ENTER  THE  DOWNSTREAM  BOUNDARY  DISTANCE  :  '\) 
READ(*,*)XWCP 
CALL  CLRSCN 
WRITE(*,34) 

34  FORMAT(//'   ENTER  VALUES  OF  BETA  '/) 
DO  105  1=1, NTR 
WRITE(*,35)I 

35  F0RMAT(10X,'   TRANSECT  ',12,':   '\) 
READ(*,*)BETA(I) 

105  CONTINUE 
WRITE(*,36) 

36  FORMAT(//'   ENTER  DECAY  RATES  AT  EACH  TRANSECT  '/) 
DO  106  1=1, NTR 
WRITE (*   37)1 

37  F0RMAT(10X,'   TRANSECT  ',12,':   '\) 
READ(*,*)RKS(I) 

106  CONTINUE 
C 

C     NOW  WRITE  ALL  THIS  DATA  TO  THE  INPUT  FILE  FOR  MIXAPPLN. FOR 
C     MAIN  PRXRAM. 
C 

0PEN(5,FILE='APINCAL.DAT',STATUS='NEW') 

WRITE(5,1)C0NNAME 

WRITE(5,1)P0LLU 

WRITE (5, *)QRS,BPWR,HPWR,UPWR, TEMP, NTR 

DO  50   1=1, NTR 

WRITE(5,*)X(I),BS(I),HS(I),US(I) 
50         CONTINUE 

WRITE(5,*)(BETA(I),I=1,NTR) 

WRITE(5,*)IFR,(QRUP(I),I=1,IFR) 

WRITE(5,*)ITEM,(TMP(I), 1=1. ITEM)  n 

WRITE(5,*)IEFL,(QEFL(I),I=1,IEFL)  I 

WRITE(5  *)IAM 

IF(IAM.EQ.1)WRITE(5,*)IPH,(PH(I),I=1,IPH) 

WRITE(5,*)QCP,CEFL,CBKG,CS,THETA,RBK,XWCP  | 

WRITE(5,*)(RKS(I),I=1,NTR)  I 

WRITE(5,55) 
55         FORMATC    3')  » 

WRITE(5,5558)CHAR(26)  | 

5558     FORMAT(A)  * 

CL0SE(5) 

IRUN=IFR*ITEM*IEFL*IPH 

RETURN 

END 
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C 
C 

SUBROUTINE  CLRSCN 
WRITE(*,101)CHAR(27),'[2J' 
101   F0RMAT(1X,A,A\) 
RETURN 
END 


I 
I 
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SDEBUG 

r**»**t **********************************  ******************************* 

c 

C  PROGRAM  PLTCRIT. FOR  USES  DATA  GENERATED  BY  THE  MIXING  ZONE      * 

C  APPLICATION  PROGRAM  "  MIXAPPLN.FOR  "  AND  PLOTS  CRITICAL         * 

C  POINTS  FOR  VARIOUS  MANAGEMENT  OPTIONS.  * 

C  WRITTEN  BY  R.  JARVIS  * 

C  GORE  &  STORRIE  1986  * 

C  * 
C************************ ***♦♦*********♦*♦*♦*♦*♦*******♦***♦*♦**♦****♦♦* 

COMMON/A/  X(10),XL(4,20),CL(4,20).XLS(4),XS(10) 
COMMON/B/  QRUP(10).QEFL(10),PH(10),RKS(10),TMP(10),ZML(10) 
CHARACTER*80  TITLE 
CHARACTER*20  POLLU 
C 

IAM=0 
MPH=1 

c 

C     OPEN  THE  FILE  "PLAPPLN"  PRODUCED  BY  "MIXAPPLN.FOR" 
C 

OPEN( 5 , FI LE= ' PLAPPLN ' ,STATUS= ' OLD ' ) 

READ(5,*)JRUN 

READ(5,2)TITLE 

READ(5,2)P0LLU 
2     FORMAT(A) 

IF(POLLU.EQ. 'AMMONIA     ')IAM=1 

READ(5,*)NTR 

DO  101  1=1, NTR 

READ(5,*)X(I) 

101  CONTINUE 
READ(5.*)MQ,(QRUP(I),I=1,MQ) 
READ(5,*)MT,(TMP(I),I=1,MT) 
READ(5.*)MF,(QEFL(I),I=1,MF) 
IF(IAH.EQ.1)READ(5,*)MPH.(PH(I),I=1,MPH) 
READ(5,*)QCP,CEFL,CBKG,CS,THETA,RBK,XWCP 
READ(5,*)(RKS(I).I=1.NTR) 

DO  1002  I=1,JRUN 
READ(5,*)ZML(I) 
DO  102  J=l,4 
READ(5,*)XL(J.I),CL(J,I) 

102  CONTINUE 
1002  CONTINUE 

CL0SE(5) 
C 

WRITE(*,103) 

103  FORMATC  PLOT  LOCATION   1=SCREEN  2=PL0TTER  3=PRINTER   :    '\) 
READ(*,*)IPL 

IF(IPL.EQ.1)CALL  PL0TS(0,0,99) 
IF(IPL.EQ.2)CALL  PL0TS(0,9600.80) 
IF(IPL.EQ.3)CALL  PL0TS(0,0. 11) 
C 

MM=0 

MRUN=0 

DO  20  JF=1,MF 

DO  20  JQ=1,MQ 

DO  20  JT=1,MT 

DO  20  JPH=1,MPH 
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CALL  SIMPLX 

IF(IPL.EQ.1)CALL  FACT0R(.75) 

IF(IPL.EQ.2)CALL  FACT0R(.9) 

MRUN=MRUN+1 

MM=MM+1 

IF(MM.EQ.2)h*1=0 

IF(MM.EQ.1)CALL  PLOT( .5,4.0,-3) 

IF(MM.EQ.O)CALL  PLOT( .0,-3.5,-3) 

CALL   PL0T(0.,0.,3) 

CALL  PL0T(0.,2.,2) 

DO   104  J=l,ll 

SJ=FL0AT(J-1)*.2 

DJ=-.15 

IF(J.EQ.2)CALL  COLOR( 1, lERR) 

IF(J.EQ.11)CALL  COLOR(0,IERR) 

CALL  PL0T(DJ,SJ,3) 

CALL  PL0T(10.,SJ,2) 

104  CONTINUE 

DO  1004  J=l,ll 

SN=FL0AT(J-1)/FL0AT(10) 

SJ=FL0AT(J-1)*.2 

SSJ=SJ-.l 

CALL  NUMBER(-.2,SSJ,.06,SN,0.,1) 
1004  CONTINUE 

CALL  SYMB0L(-.25,.15,.1,'FL0W  FRACTION  QY/QT' ,90. .19) 

IF(X(NTR).GT.XL(4,MRUN))SCALE=8.0/X(NTR) 

IF(X(NTR).LE.XL(4,MRUN))SCALE=8.0/XL(4,MRUN) 

XLS(1)=XL(1,MRUN)*SCALE 

XLS(2)=XL(2,MRUN)*SCALE 

XLS(3)=XL(3,MRUN)*SCALE 

XLS(4)=XL(4,MRUN)*SCALE 
44    CONTINUE 

DO  105  J=1,NTR 

XS(J)=X(J)*SCALE 

XSX=XS(J) 

CALL  PL0T(XSX,0.,3) 

CALL  PL0T(XSX,2.05,2) 

STR=FLOAT(J) 

CALL  NUMBER(XSX-.05,2.12,.075,STR,0.,-1) 

105  CONTINUE 

DO  106  J=1.4 

YS=FL0AT(J)*.2 

XP=XLS(J) 

CALL  C0L0R(3,IERR) 

CALL  PL0T(XP,0.,3) 

CALL  PL0T(XP,YS,2) 

CALL  SYMB0L(XP,YS,.1,1,0.,-1) 

CC=CL(J,MRUN) 

XX=XL(J,MRUN) 

CALL  C0L0R(0,IERR) 

CALL  SYMBOL(XP-.25,YS+.06,.07,'CONC=',0.,5) 

CALL  NUMBER(XP+.10,YS+.06,.07,CC,0.,2) 

CALL  SYMB0L(XP-.25,YS+.21,.07,'DIST='.0..5) 

CALL  NUMBER(XP+.10,YS+.21,.07,XX,0..2) 

106  CONTINUE 
IF(XL(4,MRUN).GT.X(NTR))SS=10.*XL(4,MRUN)/8.0 
IF(XL(4,MRUN).LE.X(NTR))SS=10.0*X(NTR)/8.0 
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SS=SS/10. 

CALL  STAXIS(.07,.07,.04,.05.2) 

CALL  AXIS(0.,0.,' DISTANCE  DOWNSTREAM  FROM  OUTFALL '.-32, 10  0  0  S 
*S) 

CALL  PL0T(10.,0.,3) 
CALL  PL0T(10.,2..2) 

CALL  SYMB0L(3..2.2..07.'TRANSECT  NUMBER' .0. .15) 
IF(MM.EQ.1)THEN 

CALL  SYMB0L(1.,3.6,.15,TITLE.0.,80) 

CALL  SYMB0L(1.,3.4..1.'PL0T  OF  CRITICAL  POINTS  FROM  PROGRAM  MIXAPP 
♦LN' .0..45) 
CALL  SYMBOL (1.. 3. 2.. 1, 'POLLUTANT  :  '.0..12) 
CALL  SYMB0L(2.5.3.2..1.P0LLU,0..20) 
ELSE 
ENDIF 

CALL  SYMB0L(0.,2.8..1.'RUN  NUMBER  '.0.,11) 
RUN=FL0AT(MRUN) 

CALL  NUMBER(1.2.2.8..1,RUN,0..0) 
CALL  SYMB0L(3. ,2.8. .075, 'RIVER  FLOW  RATE  =',0.,17) 
QQ=QRUP(JQ) 

CALL  NUMBER (5. 0.2. 8.. 075. OQ.O.. 3) 

CALL  SYMB0L(6.. 2. 8. .075, 'EFFLUENT  CONCENTRATION  ='.0..24) 
CALL  NUMBER(8.0.2.8,.075.CEFL.0..2) 

CALL  SYMB0L(6. .2.6. .075, 'MIXING  ZONE  LENGTH  ='.0..22) 
ZMM=ZML(MRUN) 

CALL  NUMBER(8.0.2.6,.075,ZMM,0.,1) 

CALL  SYMB0L(3. 0.2. 7. .075. 'EFFLUENT  FLOW  RATE  ='.0..22) 
QQ=QEFL(JF) 

CALL  NUMBER (5. 0.2. 7.. 075, QQ,0., 3) 

CALL  SYMB0L(6.0,2.7..075,'RIVER  TEMPERATURE  =',0..21) 

TTMP=TMP(JT) 

CALL  NUMBER(8.0,2.7..075.nMP.0..1) 

CALL  SYMB0L(3.0.2.6..075.'RIVER  ACIDITY 
PH  ='.0..21) 

PPH=PH(JPH) 

CALL  NUMBER(5.,2.6,.075.PPH.0..1) 

IF(^W.EQ.O)CALL  PL0T(0. .0. .-999) 
20    CONTINUE 
999   CALL  PL0T(0.,0..999) 

STOP 

END 
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SDEBUG 
(;♦♦**♦***********♦*****«♦***♦♦**♦*****♦***********♦*♦***********«**« 

C     PROGRAM  NAME  MIXPRED. FOR  DEVELOPED  FROM  MIXCALBN  TO  PREDICT     * 
C     MIXING  ZONES  IN  A  RIVER  WITH  A  PIPE  OUTFALL  * 

C**********  ******♦**♦***♦****♦*  ***^  ****♦*«****♦*****♦♦♦*******  ******** 

C  PROGRAM  NAME:  MIXPRED   *  *  STREAMTUBE  MODEL  FOR  PIPE  OUTFALL  * 

C  DEVELOPED  BY  T.  P.  H.  GOWDA,  WATER  RESOURCES  BRANCH,  MOE.  * 

C  THIS  PROGRAM  PREDICTS  LAT'L  &  LONG'L  DISTRN.  OF  CONSERVATIVE  * 

C  AND  NONCONSERVATIVE  MATERIALS  DISCHARGED  INTO  A  RIVER  FROM  * 

C  A  PIPE  OUTFALL  LOCATED  AT  BANK  OR  IN  RIVER(VERT.  LINE  SOURCE).  * 

C  PROGRAM  MODIFIED:  JUNE  1983  FOR  DILUTION  FACTOR  AND  TO  * 

C  TO  TERMINATE  CALC'NS  IF  CONCN  <  l.OE-04.  * 

c  * 

C  GORE  &     STORRIE      1986  ♦ 

C  ♦ 

C********  ********************************************************  ****** 

DIMENSION  C(50,102),CUI(50,102) 

REAL*8  X(50).XX(50).P1,P2.P3.P4,T1.T2,T3,T4,RKS(50),QY(502), 
*THETA,BPWR,HPWR,UPWR,QRTO,QT,RBT,QRS,QRUP,QEFL,RBK,CTDP,PHDR, 
*RF(50),BS(50),HS(50),US(50),BETA(50),BW,HW,UW,RKT,R,PHI,TMP, 
*B(50),H(50),U(50),BSUM(50),T0T(50).V0L(50), TEMPS, PAX1,PAX2, 
*A3,QCP,DELQ 
CHARACTER*80  TITLE 
CHARACTER*20  FILIN.FILOUT 
C 

C     INPUT  DATA 

C  MIXPRED  READS  FILE  "PINCAL.DAT"  FROM  SUBROUTINE  SETUP 
C  THE  FIRST  TIME  THROUGH  AND  MAKES  FILE  "PPINCAL.DAT"  TO 
C     BE  USED  ON  SUBSEQUENT  RUNS. 

c 

CALL  SETUP 

0PEN(5,FILE='PPINPRE.DAT',STATUS='0LD') 

0PEN(4,FILE=' PREDOUT.DAT' ,STATUS= 'NEW' ) 

0PEN(6,FILE='PLTPRED',STATUS='NEW') 
C 

C     ICAL  IS  THE  FL^G  IN  PLCALPC  THAT  INDICATES  A  PIPE  OUTFALL 
C     TO  THE  PLOTTING  PROGRAM  CONMIX.FOR 
C 

ICAL=1 

WRITE(6,*)ICAL 

WRITE(4,2) 

2  FORMAT(/'      ENTER  TITLE  OF   STUDY') 
READ(5,3)   TITLE 
WRITE(6,3)TITLE 

3  FORMAT(A) 
WRITE(4.3)TITLE 

35         WRITE(4,4) 

4  FORMATC        ENTER  QRS. BPWR.HPWR,UPWR,THETA. TEMPS, RBK  ' ) 
READ(5.*)   QRS, BPWR,HPWR,UPWR,THETA, TEMPS. RBK 

WRITE (4,400)    QRS, BPWR,HPWR,UPWR,THETA, TEMPS, RBK 

400  F0RMAT(F7.2,2X,F5.2.2X,F5.2.2X,F4.2,2X,F6.3.2X,F6.2.2X,F4.2) 
43         WRITE(4,19) 

19         FORMATC        ENTER  DESIGN  CASE:   ORUP,QEFL,CEFL,CBKG,TMP') 
READ(5,*)QRUP.QEFL,CEFL,CBKG.TMP 
WRITE(4.401)QRUP.QEFL,CEFL.CBKG,TMP 

401  F0RMAT(2X,F7.2,2X,F8.3.2X.F8.2,2X,F5.2,2X,F6.2) 
33         WRITE(4,5) 
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5  FORMATC        ENTER  NTR.NYZ.QCP' ) 
READ(5,*)    NTR.NYZ.QCP 
WRITE(6,*)NTR 
WRITE(6.*)NYZ 
WRITE(4,402)NTR,NYZ,QCP 

402  FORMAT(2X,I3,3X,I3,F5.1) 
WRITE(4,6)NTR 

6  FORMATC   ENTER  ',12,'  VALUES  OF  X,BS,HS,US') 
READ(5,*)  (X(I),BS(I),HS(I),US(I).I=1,NTR) 

DO  405  1=1, NTR 

WRITE(4,403)I,X(I),BS(I),HS(I),US(I) 
WRITE(6, 8234)1, X(I) 
8234  F0RMAT(2X,I3,2X,F8.2) 

403  F0RMAT(2X,I3,2X,F8.2.3X,F8.2.3X,F6.2,3X,F7.4) 

405  CONTINUE 

47    WRITE(4,55)NTR 

WRITE(*,55)NTR 
55    FORMAT(/'   ENTER  ',12,'  VALUES  OF  DECAY'/) 

DO  9911  1=1. NTR 

WRITE(*,5656)I 
5656  FORMATC   TRANSECT  ',12,'  :  '\) 

READ(*,*)  RKS(I) 
9911  CONTINUE 

DO  408  1=1. NTR 

WRITE(4, 406)1, RKS(I) 

WRITE(6, 406)1, RKS(I) 

406  F0RMAT(2X,I3,F9.7) 

408  CONTINUE 
45    CONTINUE 

WRITE(*,7777) 
WRITE(4,7777) 

7777  FORMAT(/'   ENTER  VALUES  OF  BETA'/) 
DO  409  1=1. NTR 

WRITE(*, 7778)1 
WRITE(4, 7778)1 

7778  FORMATC   TRANSECT  ',13,'  :  '\) 
READ(*,*)BETA(I) 
WRITE(4,411)BETA(I) 

411  F0RMAT(5X,F9.7) 
WRITE(6, 410)1, BETA(I) 

410   F0RMAT(2X.I3.F9.7) 

409  CONTINUE 
WRITE(4,52) 

52    FORMATC   ARE  YOU  CONSIDERING  UN-IONIZED  AKWONIA:  YES=1   N0=0 ' ) 

READ(5.*)  AMONIA.PH 

IF(AMONIA.EQ.O)GO  TO  101 

WRITE(4,102) 
102   FORMATC   ENTER  THE  PH  OF  UN-IONIZED  AMMONIA') 

WRITE(4,412)PH 

412  FORMATC  THE  PH  OF  UN-IONIZED  AMMONIA  IS   ',F3.1) 
101   CONTINUE 

C 

C     CALCULATE  FLOW  &  TEMP'R  PARAMETERS 

C 

DELTA=0.0001 

QT=QRUP-K)EFL 

QRTO=QT/QRS 
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CTDP=THETA**(TMP-TEMPS) 

RBT=RBK*CTDP 

DELQ=QT/NYZ 

NQ=NYZ+1 

CA=(CEFL*QEFL  /QT) 

KCP=QCP/DELCK1.5 

IF(KCP.LE.l)KCP=KCP+3 

DO   12   1=1, NTR 

BSUM(I)=0. 
12     T0T(I)=0. 

DO   14    1=1, NTR 
C 

C       CALCULATE  B,H,U  FOR  FLOW=QT,   FROM  LEOPOLD-HADDOCK  EQNS. 
C 

B(I)=BS(I)*QRTO**BPWR 

H(I)=HS(I)*QRTO**HPWR 

U(I)=US(I)*QRTO**UPWR 
C 

C       CALCULATE  WEIGHTED  MEAN  VALUES  BW,HW,UW  FROM  OUTFALL  TO  TRANSECT(I) 
C 

IF(I.GE.2)   GO  TO  60 

XX(1)=X(1) 

BW=B(1) 

HW=H(1) 

UW=U(1) 

BSUM(1)=B(1)*XX(1) 

V0L(1)=B(1)*H(1)*XX(1) 

T0T(1)=XX(1)/U(1) 

GO  TO  62 
60   11=1-1 

XX(I)=X(I)-X(I1) 

BSUM(I)=BSUM(I1)+0.5*XX(I)*(B(I1)-^6(I)) 

VOL(I)=VOL(Il)+0.25*XX(I)*(B(Il)-t-B(I))*(H(Il)+H(I)) 

T0T(I)=T0T(I1)+XX(I)/U(I) 

BW=BSUM(I)/X(I) 

HW=VOL(I)/(X(I)*BW) 

UW=QT/(BW*HW) 
C 

C   CALCULATE  PRODUCT  FUNCTION  FOR  DECAY,  RF(I). 
C 

62  RKT=CTDP*RKS(I) 

CBKX=CBKG*DEXP( -RBT*TOT ( I ) ) 

A3=(RKT*XX(I))/U(I) 

R=DEXP(-A3) 

IF(I.GE.2)  GO  TO  64 

RF(1)=R 

GO  TO  66 
64  RF(I)=RF(I1)*R 
66  CONTINUE 

PHI=BETA(I)*X(I)/BW 

PHDR=4.0*PHI 

CRPX=0.5*CA*RF(I)/DSQRT(3.1416*PHI) 

BGX=PHI*AL0G(1. /DELTA) 

SBG=SQRT(BGX) 

WRITE(4,40)    I,BETA(I),RKS(I) 
40         F0RMAT(/5X, 'TRANSECT:    ' , I2.2X , ' BETA=' ,F9.6,2X. ' RKS= ' ,F9.6,2X, 
*/7X,  'X'  ,9X,  'BW  ,9X,  'HW  ,8X,  'UW  ) 
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WRITE(4,23)  X(I),BW,HW,UW 
23    F0RMAT(2X.4(F9.3,1X)) 

WRITE(4.42) 
42    F0RMAT(4X.'QY',5X.'C(X,QY)'.5X,'CUr.9X.'C/CA'.7X,'QY/QT'.6X. 

*'DIL  FAC'/) 

DO  16  K=1,NQ 

QY(K)=FL0AT(K-1)*DELQ 

IF  (QY(K).GT.QT)  QY(K)=QT 

PAX1=(QY(K)-QCP)/QT 

PAX2=(QY(K)-K)CP)/QT 
C 

C  DETERMINE  NO.  OF  IMAGES  REQUIRED 
C 

ANl=(0.5*PAXl-SBG)-0.5 

AN2=(0.5*PAXl+SBG)+0.5 

AN3=-AN2 

AN4=-AN1 

NM1=IFIX(AN1) 

NM2=IFIX(AN2) 

NM3=IFIX(AN3) 

NM4=IFIX(AN4) 

NN1=1+NM2+IABS(NM1) 

NN2=1+IABS(NM3)+IABS(NM4) 

IF(NN1.GE.NN2)NN=NN1+1 

IF(NN1.LT.NN2)NN=NN2+1 
C 

C  COMPUTATION  OF  CONC'N  DISTR'NS. 
C 

SUM=0. 

DO  32  J=1,NN 

N=J-1 

P1=(PAX1-2.*N)**2/PHDR 

P2=(PAX2+2.*N)**2/PHDR 

CALL  PDET(P1,T1) 

CALL  PDET(P2,T2) 

IF(N.LE.O)  GO  TO  30 

P3=(PAX1+2,*N)**2/PHDR 

P4=(PAX2-2.*N)**2/PHDR 

CALL  PDET(P3,T3) 

CALL  PDET(P4.T4) 

GO  TO  32 
30  T3=0. 

T4=0. 
32  SUM=SUM+T1+T2+T3+T4 

C(I,K)=CRPX*SUM+CBKX 
C 

C  CALCULATE  UN-IONIZED  AMMONIA  CONCENTRATIONS*  OPTIONAL  ♦ 
C 

IF(AMONIA.LE.O)  60  TO  15 

PKA=0. 09018+2729. 92/(TMP+273. 2) 

PF=PKA-PH 

PCTU=1./(1.+10.**PF) 

CUI(I,K)=C(I.K)*PCTU 

15  CONTINUE 

16  CONTINUE 
NQQ=NQ 

18  CONTINUE 
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C 

C  PRINT  OUTPUT 

C 

DO  20  K=1,NQ 

RC=C(I,K)/CA 

RQ=QY(K)/QT 

CNET=C(I,K)-CBKX 

I F(CNET.LE. 0.000009)  CNET=-CEFL 

DLF=CEFL/CNET 

IF(AMONIA.LE.O)  CUI(I,K)=0.0 

WRITE(6.25)K.QY(K).C(I.K).CUI(I,K).RC.RQ.DLF 
20    WRITE(4.25)K.QY(K).C(I.I().CUI(I.K).RC.RQ.DLF 
25    F0RMAT(I2.2X,F10.2.2X.4(F9.4.2X).F10.2) 
14    CONTINUE 
999   CONTINUE 

WRITE(4,585)CHAR(26) 
585   FORMAT(A) 

WRITE(6,585)CHAR(26) 
CL0SE(4) 
CL0SE(5) 
CL0SE(6) 
WRITE(*,1234) 
1234  FORMAT(//'    THE  OUTPUT  FILE  FROM  MIXPRED. FOR  IS  CALLED'//' 
*        "  PREDOUT.DAT  "') 

STOP 

END 
C 

SUBROUTINE  PDET(P,T) 

REAL*8  P,T 

IF(P.GE.40.0)G0  TO  10 

T=DEXP(-P) 

GO  TO  12 
10  T=0.0 
12  CONTINUE 

RETURN 

END 


SUBROUTINE   SETUP 

DIMENSION  RKS(20),BS(20),ZAV(20).VAV(20).X(20).DY(15) 
CHARACTER*80  TITLE  ' 

0PEN(7,FILE='PINCAL.DAT',STATUS='0LD') 

READ(7.1)TITLE 

FORMAT(A) 

READ(7.*)QRS,BP,HP,UP.THETA 

READ(7,*)QRUP,QEFL,CEFL.CBKG,TMP 

READ(7,*)NTR,NYZ.QCP 

DO  2  1=1, NTR 

READ(7.*)X(I).BS(I),ZAV(I),VAV(I) 

CONTINUE 

CALL  CLRSCN 

WRITE(*.3)BP,HP.UP 

FORMAT(//////'        THE   EXPONENTS   FOR  THE   LEOPOLD-MADDOCK  EONS  ARE'/ 

WIDTH  EXP:  '.F5.3/'   DEPTH  EXP:  '.F5.3/'   VEL   EXP-  '  F5  3^ 
WRITE(*.4)  '   ^^' 

FORMAT(///'   DO  YOU  WISH  TO  CHANGE  THEM?   YES=1   N0=0   :  '\) 
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READ(*,*)ICH 

11  IF(ICH.EQ.1)THEN 

WRITE(*,5) 

5  FORMAT(/'  WIDTH  EXP=  '\) 
READ(*,*)BP 

WRITE(*,6) 

6  FORMAT(/'   DEPTH  EXP=   '\) 
READ(*,*)HP 

WRITE(*.7) 

7  FORMAT(/'   VEL.  EXP=   '\) 
READ(*.*)UP 

ELSE 

CONTINUE 
ENDIF 

TOT=BP+HP+UP 
CALL  CLRSCN 
IF(T0T.NE.1.)THEN 

WRITE(*,10) 
10  F0RMAT(/////'  THE   EXPONENTS  MUST  SUM  TO   1.0'/'      RE-ENTER  TH 

♦EM') 

GO  TO   11 
ELSE 

CONTINUE 
ENDIF 

CALL  CLRSCN 
WRITE(*,8)QRUP,QEFL 

8  FORMAT(/////'  RIVER  FLOW  RATE  ABOVE  OUTFALL= ' ,F10.2/ '  OUTFAL 
*L  FLOW  RATE=  ',F10.2//'  DO  YOU  WISH  TO  CHANGE  EITHER?  YES= 
*1  N0=0   '\) 

READ(*,*)IC 
IF(IC.EQ.1)THEN 
WRITE(*,9) 

9  FORMATC   RIVER  FLOW  RATE=   '\) 
READ(*,*)QRUP 

WRITE(*,12) 

12  FORMATC   OUTFALL  FLOW  RATE=   '\) 
READ(*.*)QEFL 

QRS=QRUP+QEFL 
ENDIF 

CALL  CLRSCN 
C 

WRITE(*,28)CEFL 
28    FORMAT(/////'    THE  EFFLUENT  CONCENTRATION  IS  :',F10.2.//' 
*D0  WISH  TO  CHANGE  IT?  YES=1  N0=0  '\) 
READ(*,*)IC 
IF(IC.EQ.1)THEN 
WRITE(*,30) 

30  FORMAT(/'   THE  NEW  EFFLUENT  CONC.  =  '\) 
READ(*,*)CEFL 

ENDIF 

CALL  CLRSCN 

WRITE(*,31)CBKG 

31  FORMAT(/////'     THE  BACKGROUND  CONCENTRATION  IS  :'.F10.2.//' 
*D0  YOU  WISH  TO  CHANGE  IT?  YES=1  NO  =0  '\) 

READ(*,*)IC 
IF(IC.EQ.1)THEN 
WRITE(*,32) 


MIXPRED  NOV.   21,   1986       12:06  PAGE     7 


32  FORMAT(/'        ENTER  NEW  BACKGROUND  CONC:      '\) 

READ(*,*)CBKG 

ENDIF 

CALL  CLRSCN 
C 

WRITE(*.20) 

20  FORMAT(/////'  ENTER  A  DECAY  RATE   FOR  THE   RIVER  BACKGROUND:    '\) 
READ(*,*)RBK 

WRITE(*,21) 

21  FORMAT(/'    AT  WHAT  TEMPERATURE  IS  THIS  RATE  KNOWN?  IN  C  :  '\) 
READ(*,*)TEMP 

WRITE(*.22) 

22  FORMAT(/'    WHAT  IS  THE  RIVER  TEMPERATURE?  IN  C  :   '\) 
READ(*,*)TMP 

CALL  CLRSCN 
WRITE(*,23) 

23  FORMAT(/////'     THE  OUTFALL  IS  AT  SHORE  '// '    DO  YOU  WISH  TO  C 
*HANGE  IT?  YES=1   N0=0   '\) 

READ(*,*)IC 
IF(IC.EQ.1)THEN 
WRITE(*,24) 

24  FORMAT(/'   ENTER  THE  DISTANCE  OF  THE  OUTFALL  FROM  THE  BANK  :  '\ 

*) 

READ{*,*)YOUT 

0PEN(9.FILE=' SCALE. DAT' ,STATUS= 'OLD ' ) 
READ(9,*)(DY(I).I=1,11) 
CL0SE(9) 
IG=0 
9079    IG=IG+1 

IF(YOUT.LE.DY(IG))THEN 

A=DY(IG)-Y0UT 

DELY=DY(IG)-DY(IG-1) 

A1=FL0AT(IG-1)/10. 

QCP=(A1-.1*(A/DELY))*QRS 

ELSE 

GO  TO  9079 
ENDIF 
ENDIF 
C 

CALL  CLRSCN 
WRITE(*,50) 

50  FORMAT(/////'  DO  YOU  WISH  TO  CONSIDER  AMMONIA?     YES=1     N0=0      '\ 
*) 

READ(*.*)IAM 
IF(IAM.EQ.1)THEN 
WRITE(*.51) 

51  FORMAT(/'        ENTER  PH     '\) 
READ(*,*)PH 

ELSE 
PH=7.0 

ENDIF 

CALL  CLRSCN 
C 

OPEN( 7, FILE=' PPINPRE.DAT', STATUS- 'NEW') 

WRITE(7.1)TITLE 

WRITE(7.111)QRS,BP.HP,UP.THETA,TEMP,RBK 
111       F0RMAT(2X,F9.2,4(2X.F5.3),2X,F6.2,2X.F10.7) 
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WRITE(7,112)QRUP,QEFL,CEFL,CBKG.TMP 
112   F0RMAT(2X.F9.2.2X.F8.2,2X.F8.3,2X,F6.3,2X,F6.2) 
WRITE(7,*)NTR,NYZ,QCP,Y0UT 
DO  40  1=1, NTR 

WRITE(7.*)X(I),BS(I),ZAV(I),VAV(I) 
40    CONTINUE 

WRITE(7,*)IAM,PH 
WRITE(7,585)CHAR(26) 
585   FORMAT(A) 
CL0SE(7) 
RETURN 
END 
C 
C 

SUBROUTINE  CLRSCN 
WRITE(*.100)CHAR(27),'[2J' 
100     F0RMAT(1X,A,A\) 
RETURN 
END 
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SDEBUG 


C 

C  PROGRAM  MIXCADIF. FOR  IS  THE  MICRO-EDITION  OF  MIXCADIF  BY  H  GOWDA 

C 

C  PROGRAM  NAME:  MIXCADIF  *  * 

C  DEVELOPED  BY  T.  P.  H.  GOWDA,  WATER  RESOURCES  BRANCH,  MOE. 

C  THIS  PROGRAM  PREDICTS  LAT ' L  &  LONG'L  DISTRN.  OF  CONSERVATIVE 

C  OR  NONCONSERVATIVE  MATERIAL  DISCHARGED  INTO  A  RIVER  FROM 


C     A  DIFFUSER  OUTFALL 

C 

C  GORE  &  STORRIE     1986 

C 


IMPLICIT  REAL*8    (A-H,0-Z) 

COMMON/C/  C(25,101),CUI(25.101) 

COMMON/VAR/  X(25) ,XX(25) ,RKS(25) .QY(lOl) .RF(25) ,BS(25) ,HS(25) 

*.US(25).BETA(25),B(25).H(25).U(25).BSUM(25)   TOT  25   .V0L(25) 

CHARACTER*80  TITLE  ' 

C 

C  INPUT  DATA 
C 

CALL  CLRSCN 

CALL  SETUP 

0PEN(6,FILE='CAD0UT.DAT',STATUS='NEW') 

0PEN(5,FILE='PPINCAD.DAT',STATUS='0LD') 

0PEN(7,FILE='PLTPRED',STATUS='NEW') 

WRITE(6,2) 

2  FORMAT(/'  ENTER  TITLE  OF  STUDY') 
READ(5,3)    TITLE 
WRITE(6,3)TITLE 

ICAL=2 

WRITE(7,*)ICAL 

WRITE(7,3)TITLE 

3  FORMAT(A) 
35         WRITE(6.4) 

4  FORMATC  ENTER  QRS. BPWR.HPWR.UPWR.THETA, TEMPS, RBK' ) 
READ(5,*)   QRS,BPWR,HPWR,UPWR,THETA, TEMPS, R6K 

WRITE (6, 3000)QRS.BPWR,HPWR.UPWR,THETA, TEMPS, RBK 

3000  FORMAT(2X.F10.3.2X,F6.3.2X,F6.3,2X,F5.3.2X.F7.4,2X,F7.4,2X,F5.2) 

19    FORMAT('    ENTER  DESIGN  CASE:  QRUP,QEFL,CEFL,CBKG.TMP' ) 
READ(5,*)QRUP,0EFL,CEFL,CBKG,TMP 
WRITE(6,3001)QRUP,QEFL,CEFL.CBKG,TMP 

3001  f^0RMAT(2X,F10.3.2X,F6.2.2X.F8.2,2X.F5.2,2X,F6.2) 
33         WRITE(6.5)  .  .      .  ; 

5  FORMATC        ENTER  NTR,NYZ.N0UT,QCP1,QCP2') 
READ(5,*)    NTR,NYZ,N0UT,QCP1,QCP2 
WRITE(7,*)NTR 

WRITE(7,*)NYZ 

WRITE (6, 3003) NTR,NYZ,N0UT.QCP1.QCP2 
3003     F0RMAT(2X,I3,2X.I3,2X,I3.2X,F7.2,2X.F7.2) 
WRITE(6.6) 

6  FORMATC        ENTER  NTR  VALUES  OF   X,BS.HS,US') 
READ(5.*)    (X(I).BS(I),HS(I).US(I).I=1,NTR) 
DO   3004    1=1, NTR 

WRITE(7, 2345)1, X(I) 
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2345     F0RMAT(2X,I2,5X,F8.2) 

WRITE(6.3005)X(I).BS(I).HS(I).US(I) 

3005  F0RMAT(2X,F8.2,2X,F8.2.2X.F5.2,2X.F6.3) 
3004     CONTINUE 

47    WRITE(*.55)NTR 

55  FORMAT(//'    ENTER  ',12.'  VALUES  OF  DECAY  COEFFICIENT'/) 
DO  3077  1=1, NTR 

WRITE(*, 7717)1 
7717  FORMATC     TRANSECT  ',12,'  :  '\) 

READ(*,*)RKS(I) 
3077  CONTINUE 

CALL  CLRSCN 

DO  3006  1=1, NTR 

WRITE(6,3007)RKS(I) 

WRITE(7,*)I,RKS(I) 
3007  F0RMAT(3X,F10.8) 

3006  CONTINUE 
45    CONTINUE 

CALL  CLRSCN 

WRITE(*,8)NTR 
8  FORMAT(//'        ENTER   ',12,'   VALUES  OF  BETA'/) 

DO  409   1=1, NTR 

WRITE(*, 7778)1 
7778  FORMATC     TRANSECT  ',12,'  :  '\) 

READ(*,*)BETA(I) 

WRITE(6, 410)1, BETA(I) 

WRITE(7, 410)1, BETA(I) 
410  F0RMAT(2X,I3,2X,F9.7) 
409   CONTINUE 

WRITE(6,52) 
52    FORMATC    UN-IONIZED  AMMONIA:   ENTER  1  FOR  YES,   0  FOR  NO 
AND  EN 

*TER  PH  VALUE  ON  THE  SAME  LINE') 

READ(5,*)  AMONIA,PH 

WRITE(6,3010)AM0NIA,PH 
3010  F0RMAT(2X.F3.2,3X,F5.2) 
C  CALCULATE  FLOW  &  TEMP'R  PARAMETERS 

DELTA=0.0001 

QT=QRUP-K)EFL 

QRTO=QT/QRS 

CTDP=THETA**(TMP-TEMPS) 

RBT=RBK*CTDP 

DELQ=QT/NYZ 

NQ=NYZ+1 

KCP=QCP2/DELQ+1.5 

WRITE(6,24)TITLE 
24  F0RMAT(//2X,A) 

WRITE (6,56)   QRUP,CBKG,TMP,QEFL,CEFL,QCP1,QCP2 

56  F0RMAT(2X, 'UPSTREAM  FLOW= ' ,F8.3, '  BACKGROUND  CONC.= ' ,F8.3, '  DES 
*IGN  TEMP=',F5.1/2X, 'EFFLUENT  FLOW=' ,F8.3, '  EFFLUENT  C0NC.=',F8.2 
V2X,'DIFFUSER  OUTFALL  LOCATED  BETWEEN ',F7. 2,'      AND      ',F7.2/) 

CA=CEFL*QEFL/QT 
DO   12   1=1, NTR 
BSUM(I)=0. 
12         T0T(I)=0. 

DO  14   1=1, NTR 
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C       CALCULATE  B.H.U  FOR  FLOW=QT.    FROM  LEOPOLD-HADDOCK  EQNS. 
C 

B(I)=BS(I)*QRTO**BPWR 
H(I)=HS(I)*QRTO**HPWR 
U(I)=US(I)*QRTO**UPWR 
C 

C   CALCULATE  WEIGHTED  MEAN  VALUES  BW.HW.UW  FROM  OUTFALL  TO  TRANSECTd) 
C  ' 

IF(I.GE.2)  GO  TO  60 
XX(1)=X(1) 
BW=B(1) 
HW=H(1) 
UW=U(1) 

BSUM(1)=B(1)*XX(1) 
V0L(1)=B(1)*H(1)*XX(1) 
T0T(1)=XX(1)/U(1) 
GO  TO  62 
60    11=1-1 

XX(I)=X(I)-X(I1) 

BSUM(I)=BSUM(I1)+0.5*XX(I)*(B(I1)+B(I)) 
V0L(I)=V0L(I1)+0.25*XX(I)*(B(I1)+B(I))*(H(I1)+H(I)) 
T0T(I)=T0T(I1)+XX(I)/U(I) 
BW=BSUM(I)/X(I) 
HW=VOL(I)/(X(I)*BW) 
UW=QT/(BW*HW) 
C 

C       CALCULATE   PRODUCT  FUNCTION  FOR  DECAY.   RFfl). 
C 
62         RKT=CTDP*RKS(I) 

CBKX=CBKG*DEXP( -RBT*TOT( I ) ) 
A3=(RKT*XX(I))/U(I) 
R=DEXP(-A3) 
IF(I.GE.2)   GO  TO  64 
RF(1)=R 
GO  TO  66 
64  RF(I)=RF(I1)*R 

66         CONTINUE 

PHI=BETA(I)*X(I)/BW 
PHD2=2.0*DSQRT(PHI) 
CRPX=0.5*CA*RF(I)*QT/(QCP2-QCP1) 
WRITE(6,40)    I,BETA(I).RKS(I) 
40  F0RMAT(/5X. 'TRANSECT:    ' . I2,2X, 'BETA=' ,F9.6.2X, ' RKS= ' .F9.6./7X. ' X' , 

*9X,'BW',9X,'HW',8X,'UW') 
WRITE(6,23)    X(I), BW.HW.UW 
23  F0RMAT(2X.4(F9.3.1X)/) 

WRITE(6,42) 
42  F0RMAT(4X.  ■)(•  .4X.  'QY'  ,5X.  •C(X.QY)  '  .SX.'CUI '  .9X.  'C/CA'  .6X.  "QY/QT'  .2 

♦X.'DIL   FAC'//) 
DO   16  K=1.NQ 
QY(K)=(K-1)*DELQ 
IF    (QY(K).GT.QT)   QY(K)=QT 
DP1=(QCP1-QY(K))/QT 
DP2=(0CP2-QY(K))/QT 
DP3  =  (QCPl-fOY(K))/QT 
DP4=(QCP2-K)Y(K))/QT 
SUM1=0.0 
SUM2=0.0 
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DO  30  NN=1,9 

N=NN~1 

T1=DERF((DP2+2.0*N)/PHD2) 

T2=DERF((DP1+2.0*N)/PHD2) 

T3=DERF ( (0P4+2.0*N) /PHD2) 

T4=DERF((DP3+2.0*N)/PHD2) 
30    SUM1=SUM1+T1-T2+T3-T4 

DO  32  N=l,8 

T5=DERF((DP2-2.0*N)/PHD2) 

T6=DERF((DP1-2.0*N)/PHD2) 

T7=DERF((DP4-2.0*N)/PHD2) 

T8=DERF((DP3-2.0*N)/PHD2) 
32    SUM2=SUM2+T5-T6+T7- T8 

C(I.K)=CRPX*(SUM1+SUM2)+CBKX 
C 

C  CALCULATE  UN-IONIZED  AMMONIA  CONCENTRATIONS*  OPTIONAL  * 
C 

IF(AMONIA.LE.O)  GO  TO  15 

PKA=0. 090 18+2729. 92/(TMP+273. 2) 

PF=PKA-PH 

PCTU=l./(l.-t-10.**PF) 

CUI(I,K)=C(I,K)*PCTU 

15  CONTINUE 
GOTO  16 
NQQ=K 

GO  TO  18 

16  CONTINUE 
NQQ=NQ 

18    CONTINUE 

C 

C   PRINT  OUTPUT 

C 

DO  14  K=1,NQ 

CNET=C(I,K)-CBKX 

RC=CNET/CA 

IF(CNET.LE. 0.000009)  CNET=-CEFL 

DIL=CEFL/CNET 

RQ=QY(K)/QT 

IF(AMONIA.LE.O)  CUI(I,K)=0.0 

WRITE(6,25)  K,QY(K),C(I,K),CUI(I,K),RC,RQ,DIL 

WRITE(7,25)K,QY(K),C(I,K).CUI(I,K),RC,RQ,DIL 
25    F0RMAT(3X,I2,F8.2,2X,4(F9.4,2X,F9.1)) 
14    CONTINUE 

WRITE(5,10101)CHAR(26) 

WRITE(6.10101)CHAR(26) 

WRITE(7,10101)CHAR(26) 
10101  F0RMAT(1X,A.\) 

CLOSE (5) 

CLOSE (6) 

CL0SE(7) 

CALL  CLRSCN 

WRITE(*,1324) 
1324  F0RMAT(/////'    THE  OUTPUT  FILE  FOR  MIXCADIF  IS'//'    "  CADOUT.D 
*AT  "') 

STOP 

END 
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C 
C 

FUNCTION  OERF(X) 

REAL*8  X,XX,A1,A2,A3,A4,A5J,P,AA,DERF,SX 

Al=. 254829592 

A2=-. 284496736 

A3=l. 421413741 

A4=-l. 453152027 

A5=l. 061405429 

P=. 3275911 

XX=DABS(X) 

IF(XX.LT.10.)G0  TO  800 

DERF=1. 

GO  TO  801 

800  T=1./(1.+P*XX) 
AA=T*(A1+T*(A2+T*(A3+T*(A4+T*A5)))) 
DERF=1.-AA*DEXP(-XX*XX) 

801  IF(X.LT.O)DERF=-DERF 
RETURN 

END 
C 

c 
c 

SUBROUTINE  SETUP 

DIMENSION  RKS(20),BS(20),ZAV(20).VAV(20),X(20),DY(15) 

CHARACTER*80  TITLE 

OPEN( 7, FI LE=' PINCAL.DAT', ST ATUS=' OLD') 

0PEN(8,  F I  LE  =  ' PINCAD.DAT', ST  ATUS=' NEW) 

READ(7.1)TITLE 

WRITE(8,1)TITLE 

1  FORMAT(A) 
READ(7.919)QRS,BP,HP,UP,THETA 
WRITE(8.919)QRS,BP,HP,UP.THETA 

919  F0RMAT(2X,F9.2,2X,F5.3,2X,F5.3.2X,F5.3.2X,F6.3) 
READ(7,920)QRUP,QEFL,CEFL,CBKG.TMP 
WRITE(8.920)QRUP,QEFL,CEFL,CBKG,TMP 

920  F0RMAT(2X,F9.2,2X,F9.2,2X,F10.3,2X,F10.3,2X,F6.1) 
READ(7,*)NTR,NYZ,QCP 

C 

c 
c 

DO  2  1  =  1, NTR 
READ(7,*)X(I),BS(I),ZAV(I),VAV(I) 

2  CONTINUE 
CALL  CLRSCN 
WRITE(*.4008) 

4008  FORMATC     HYDRODYNAMIC  DATA  ENTRY  AREA'/'     ***************** 
*************^*^***,) 

WRITE(*,3)BP,HP,UP 

3  FORMAT(////'   THE  EXPONENTS  FOR  THE  LEOPOLD-HADDOCK  EQNS  ARE'/ 
♦'  WIDTH  EXP:  ',F5.3/'  DEPTH  EXP:  ',F5.3/'  VEL.  EXP:  ',F5.3) 

WRITE(*,4) 

4  FORMAT(/'   DO  YOU  WISH  TO  CHANGE  THEM?   YES=1   N0=0 ' ) 
READ(*,*)ICH 

11    IF(ICH.EQ.1)THEN 
WRITE(*.5) 

5  F0RKAT(/'   WIDTH  EXP=   '\) 
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READ(*,*)BP 
WRITE(*,6) 

6  FORMAT(/'      DEPTH   EXP=      '\) 
READ(*,*)HP 

WRITE(*.7) 

7  FORMAT(/'      VEL.    EXP=        '\) 
READ(*.*)UP 

ENDIF 

TOT=BPt-HPHJP 

IF(T0T.NE.1.)THEN 
WRITE(*.10) 
10  FORMAT ('  THE   EXPONENTS  MUST  SUM  TO  1.0'/'      RE-ENTER  THEM') 

GO  TO   11 

ENDIF 

CALL  CLRSCN 

WRITE(*,4000) 
4000     FORMATC  ENTRY  AREA  FOR  DESIGN  PARAMETERS'/'  ********* 

icic  ^^  Altitik  ic^  irir  irit  'ir'ir  iric  1c  ^t  Ik  it  icic  icli  '  \ 

WRITE(*,8)QRUP.QEFL 

8  FORMAT(////'  RIVER  FLOW  RATE  ABOVE  OUTFALL=' .F10.2/'  OUTFALL 
♦FLOW  RATE=  ',F10.2//'  DO  YOU  WISH  TO  CHANGE  EITHER?  YES=1 
*     N0=0      '\) 

READ(*,*)IC 
IF(IC.EQ.1)THEN 
WRITE(*,9) 

9  FORMATC        RIVER  FLOW  RATE=      '\) 
READ(*,*)QRUP 

WRITE(*,12) 
12  FORMATC        OUTFALL  FLOW  RATE=      '\) 

READ(*,*)QEFL 
QRS=QRUPt-QEFL 
ENDIF 

CALL  CLRSCN 
C 
C 

WRITE(*,1110) 
1110     FORMAT(////'  ENTER  DIFFUSER  OUTFALL  LOCATION'/'  IN  METERS  FR 

*0M  BANK') 
WRITE(*,1011) 

1011  F0RMAT(/'        INPUT  DISTANCE   FORM  REF.   BANK  OF  PROXIMAL  END   :    '\) 
READ(*.*)Y0UT1 

WRITE(*.1012) 

1012  FORMAT(/'        INPUT  DISTANCE  FROM  REF.   BANK  OF  DISTAL  END   :    '\) 
READ(*,*)Y0UT2 

WRITE(*,1013) 

1013  FORMATC/'        INPUT  NUMBER  OF   PORTS  ON  DIFFUSER   :    '\) 
READ(*.*)NOUT 

0PEN(9,FILE= 'SCALE. DAT', STATUS= 'OLD') 
READ(9.*)(DY(I),I=1,11) 
CLOSE (9) 
IG=0 
9097     IG=IG+1 

IF(Y0UT1.LE.DY(IG))THEN 

A=DY(IG)-Y0UT1 

DELY=DY(IG)-DY(I6-1) 

A1=FL0AT(IG-1)/10. 

QCP1=(A1-.1*(A/DELY)) 
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QCP1=QCP1*QRS 
ELSE 

GO  TO  9097 
ENDIF 
IG=0 
9098     IG=IG-t-l 

IF(Y0UT2.LE.DY(IG))THEN 

A=DY(IG)-Y0UT2 

DELY=DY(IG)-DY(IG-1) 

A1=FL0AT(IG-1)/10. 

QCP2=(A1-.1*(A/DELY)) 

QCP2=QCP2*QRS 
ELSE 

GO  TO  9098 
ENDIF 
C 

WRITE(8,661)NTR,NYZ,N0UT.0CP1,QCP2 

661  F0RMAT(2X,I2,2X,I2,2X,I4,2X,F8.2,2X,F8.2) 
DO  662    1=1, NTR 
WRITE(8,*)X(I),BS(I),ZAV(I),VAV(I) 

662  CONTINUE 
C 

CALL  CLRSCN 
WRITE(*,4002) 

4002  FORKATC     CONCENTRATION  VALUE  ENTRY  AREA'/'    ***************^ 

WRITE(*,28)CEFL 
28    FORMAT(/'    THE  EFFLUENT  CONCENTRATION  IS  :'.FI0.2,/'    DO  WISH 
*T0  CHANGE  IT?   YES=1  N0=0   '\) 
READ(*,*)IC 
IF(IC.EQ.1)THEN 
WRITE(*,30) 

30  FORMAT(/'    THE  NEW  EFFLUENT  CONC.  =   '\) 
READ(*.*)CEFL 

ENDIF 
WRITE(*,31)CBKG 

31  FORMAT(/'     THE  BACKGROUND  CONCENTRATION  IS  :'.F10.2,/'   00  YOU 
♦WISH  TO  CHANGE  IT?  YES=1   NO  =0   '\) 

READ(*,*)IC 
IF(IC.EQ.1)THEN 
WRITE(*,32) 

32  FORMAT(/'    ENTER  NEW  BACKGROUND  CONC:   '\) 
READ(*.*)CBKG 

ENDIF 

CALL  CLRSCN 

WRITE(*,4003) 

4003  FORMATC  DECAY   RATE   ENTRY  AREA'/'  *************^****i 

C 

WRITE(*.20) 

20  FORMAT(/'    ENTER  A  DECAY  RATE  FOR  THE  RIVER  BACKGROUND:  '\) 
READ(*,*)RBK 

WRITE(*.21) 

21  FORMAT(/'    AT  WHAT  TEMPERATURE  IS  THIS  RATE  KNOWN?  IN  C  :  '\) 
READ(*.*)TEMP 

WRITE(*.22) 

22  FORMAT(/'    WHAT  IS  THE  RIVER  TEMPERATURE?  IN  C  :   '\) 
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READ(*,*)"mP 
C 

WRITE(*,50) 

50  FORMAT(/'  DO  YOU  WISH  TO  CONSIDER  AMMONIA?     YES=1     N0=0      '\) 
READ(*,*)IAM 

IF(IAM.EQ.1)THEN 
WRITE(*,51) 

51  FORMAT(/'        ENTER  PH      '\) 
READ(*,*)PH 

ELSE 

PH=7.0 
ENDIF 
C 

0PEN(7,FILE='PPINCAD.DAT',STATUS='NEW') 

WRITE(7,1)TITLE 

WRITE(7,111)QRS,BP,HP,UP,THETA,TEMP,RBK 

111  F0RMAT(2X,F9.2,3(2X,F5.3),2X,F5.3,2X,F5.2,2X,F10.8) 
WRITE(7,112)QRUP,QEFL,CEFL,CBKG,TMP 

112  F0RMAT(2X,F9.2.2X,F8.2,2X,F8.2,2X,F8.2,2X.F5.2) 
WRITE(7,113)NTR,NYZ,N0UT,QCP1,QCP2,Y0UT1,Y0UT2 

113  F0RMAT(2X,I2,2X,I3,2X,I3,2X,F8.2,2X,F8.2.2X,F8.2,2X,F8.2) 
DO  4033  1=1, NTR 
WRITE(7,*)X(I).BS(I),ZAV(I),VAV(I) 

4033  CONTINUE 

WRITE(7,*)IAM,PH 

WRITE(7,4O07)CHAR(26) 
4007  F0RMAT(1X,A) 

CL0SE(7) 

RETURN 

END 
C 

c 

c 

SUBROUTINE  CLRSCN 
WRITE(*,101)CHAR(27),'[2J' 
101   F0RMAT(1X,A,A\) 
RETURN 
END 
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c 

C  PROGRAM  CONMIX. FOR  * 

C  THIS  PROGRAM  PLOTS  THE  OUTPUTS  FROM  MODELS  MIXCALBN  * 

C  AND  MIXCADIF.  IT  PRODUCES  THREE  TYPES  OF  PLOTS.  * 

C  1.   LATERAL  CONCENTRATION  PLOTS  AT  EACH  TRANSECT  * 

C  2.   A  TWO  DIMENSIONAL  CONTOUR  PLOT  < 

C  3.   A  THREE  DIMENSIONAL  MESH  PLOT  * 

C  * 

C  WRITTEN  BY  R.  JARVIS  * 

C  < 

C  GORE  &  STORRIE  1986  * 

COMMON/AA/  JK 

COMMON/B/  CC(50,50) 

COMMON/L/  XG(50).C(10, 50), TITLE 

COMMON/PLA/  XH(150) ,DLX(50) ,ZLEV(10) 

COMMON/PASS/  XDX(10),BS(10),HS(10),US(10),QRS,QRUP,QEFL,CEFL,CBKG 
*,QCP,P0LLU,ICAL,N0UT,QCP1,QCP2.Y0UT1,Y0UT2,Y0UT 

INTEGER*2   LWGT( 10) ,KNXT(2000) ,LDIG(10) , ICOL,IROW,NX,NY,NARC 
l.NDIV.NLEV.NPTS.NRNG 

CHARACTER*80  TITLE 

CHARACTER*20   FILIN 

CHARACTER*30  POLLU 
C 

NLEV=10 

NLEVL=NLEV 

IC0L=50 

IR0W=50 

c 

WRITE(*,6777) 

6777  F0RMAT(//'  PLOT  THE   FOLLOWING  TYPE   OF   PLOT        1=YES       0=N0'/) 
WRITE(*,6778) 

6778  FORMAT(/'  CROSS  STREAM  CONCENTRATION  PROFILES:      '\) 
READ(*,*)ILAT 

WRITE(*,6779) 

6779  FORMAT(/'       2  DIMENSIONAL  CONTOUR  PLOT:   '\) 
READ(*,*)I2D 

WRITE(*,6780) 

6780  FORMAT(/'       3  DIMENSIONAL  MESH  PLOT:   '\) 
READ(*,*)I3D 

C 

C     READ  PLCALPC 

C 

0PEN(5.FILE='PLTPRED',STATUS='0LD') 

READ(5,*)ICAL 

READ(5,112)TITLE 
112        FORMAT(A) 

READ(5.*)NTR 

READ(5.*)NYZ 

DO  606   1=1, NTR 

READ(5,*)II,XG(I) 
606   CONTINUE 

NTR2=2*NTR 

DO  7779   11=1, NTR2 

READ(5.*) ID. DUMMY 
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7779 

CONTINUE 

XHIGH=XG(NTR) 

DO  900  J=1,NTR 

DO  905  1=1.11 

READ(5.*)II,F,C(J.I) 

905 

CONTINUE 

900 

CONTINUE 

C 
C 

CL0SE(5) 

READ  PINCAL.DAT  FOR 

C 

IF( ICAL.EQ.1)0PEN(5,FILE=' PPINPRE.DAT', ST ATUS=' OLD') 

IF(ICAL.EQ.2)0PEN(5,FILE='PPINCAD.DAT',STATUS='0LD') 

READ(5.112)DTITLE 

READ(5,*)QRS,BEX,DEX,VEX,TEMC 

READ(5,*)QRUP,QEFL,CEFL,CBKG 

IF(ICAL.EQ.1)READ(5,*)NNTR,NNYZ,QCP.Y0UT 

IF(ICAL.EQ.2)READ(5,*)NNTR,NNYZ.N0UT,QCP1,QCP2.Y0UT1,Y0UT2 

DO  7336  I=1.NNTR 

READ(5,*)XDX(I).BS(I),HS(I),US(I) 
7336  CONTINUE 

CLOSE (5) 
C 

WRITE(*,1011) 
1011   F0RMAT(//'    ENTER  POLLUTANT  NAME  (30  LEHERS  MAX)  :  '\) 

READ(*,112)P0LLU 
C 

DLX(1)=XG(1) 

DO  502  J=2,NTR 

DLX(J)=XG(J)-XG(J-1) 
502   CONTINUE 

NTRP1=NTR+1 

SCALE=XHIGH/50. 

DO  3  1=1,50 

XH(I)=FLOAT(I)*SCALE 
3     CONTINUE 

1  =  1 

DO  501  J=l,50 

IF(XH(J).GT.XG(I))I=I+I 

IF(I.EQ.NTRPl)  GO  TO  501 

A=XG(I)-XH(J) 

B=DLX(I)-A 

DO  500  K=l,ll 

IF(I.EQ.1)THEN 

CC(J.K)=0. 

ELSE 

CC(J,K)=(B*C(I,K)+A*C(I-1.K))/DLX(I) 

ENDIF 

500  CONTINUE 

501  CONTINUE 
C 

WRITE(*,30) 
30         FORMAT(//'    PLOT  LOCATION   1=SCREEN  2=PL0nER  3=PRINTER   :    '\) 
READ(*,*)JK 

IF(JK.EQ.1)CALL  PL0TS(0,0,99) 
IF(JK.EQ.1)CALL  FACT0R(.5) 
IF(JK.EQ.2)CALL  PLOTS(0.9600,80) 
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IF(JK.EQ.3)CALL   PL0TS(0,0.11) 

IF(JK.NE.2.AND.JK.NE.1.AND.JK.NE.3)G0  TO  35 

CALL  SIMPLX 
C 

IF(ILAT.EQ.1)THEN 

CALL  TRANS(NTR) 

CALL  PL0T(0.,0.,-999) 
ELSE 

CONTINUE 
ENDIF 
IF(I3D.EQ.1)THEN 

CALL  THREED(NTR) 

CALL  PL0T(0.,0..-999) 
ELSE 

CONTINUE 
ENDIF 

IF(I2D.EQ.1)CALL  TWOD(NTR,XHIGH) 
CALL  PL0T(0..0.,999) 
35    STOP 

END 
C 

C  SUBROUTINE  TWOD  PLOTS  THE  TWO  DIMENSIONAL  CONTOUR  PLOT 

SUBROUTINE   TWOD(NTR.XHIGH) 

COMMON/AA/  JK 

COMMON/A/   X(3000),Y(3000) 

COMMON/D/  Z(50.50),ZZ(30O0) 

COMMON/B/  CC(50,50) 

COMMON/L/  XG(50),C(10, 50). TITLE 

COMMON/KKL/  ZPIJ(2000) .DLX(50) .ZLEV(IO) 

^Tr^Jn^'^^^^  XDX(10).BS(10).HS(10).US(10),QRS.QRUP.QEFL  CEFL  CBKG 
*,QCP.P0LLU.ICAL.N0UT.QCP1.QCP2.Y0UT1.Y0UT2,Y0UT        '^'''•'^'^^'^'^'^^ 

CHARACTER*30  SNAME.POLLU 

CHARACTER*80  TITLE 


NLEV=10 

NLEVL=NLEV 

IC0L=50 

IR0W=50 

NXL=50 

CMIN=99999.0 

CMAX=0. 

DO  99   1=1, NXL 

DO  99  J=l.ll 

IF(CC(I.J).LT.CMIN)CMIN=CC(I.J) 
IF(CC(I.J).GT.CMAX)CMAX=CC(I.J 
99         CONTINUE 

IF(CMAX.LT.1)LG0=3 
IF(CMAX.GT.1)LG0=2 
IF(CMAX.GT.10)LG0=1 
DO  98   1=1.10 
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LDIG(I)=LGO 
98    CONTINUE 

LWGT(1)=101 
LWGT(2)=201 
LWGT(3)=301 
LWGT(4)=401 
LWGT(5)=501 
LWGT(6)=601 
LWGT(7)=701 
LWGT(8)=801 
LWGT(9)=101 

cr=cmax-cmin 

zlev(1)=cr/10.+cmin 

zlev(2)=2.*cr/10.+cmin 

zlev(3)=3.*cr/io.-k:min 

zlev(4)=4.*cr/10.+cmin 

ZLEV(5)=5.*CR/10.+CMIN 

ZLEV(6)=6.*CR/10.+CMIN 

ZLEV(7)=7.*CR/10.+CMIN 

ZLEV(8)=8.*CR/10.+CMIN 

ZLEV(9)=9.*CR/10.+CMIN 

ZLEV(10)=CMAX 

DO  801  L=l,2000 

ZPIJ(L)=0. 

KNXT(L)=0 

X(L)=0. 

Y(L.)=0. 

ZZ(L)=0. 
801   CONTINUE 

CALL  PLOT(0.5,0.7,-3) 

CALL  SIMPLX 

IF(JK.EQ.1)CALL  FACT0R(.8) 

K=0 

NX=50 

NXL=NX 

NY=11 

NPTS=NX*11 

DX=(XHIGH-XL0W)/FL0AT(NX-1) 

DY=(YHIGH-YL0W)/FL0AT(NY-1) 

DO  10  1  =  1, NXL 

DO  11  J=l,ll 

K=K+1 

X(K)=FL0AT(I-1)*DX 

Y(K)=FL0AT(J-1)*0Y 

ZZ(K)=CC(I,J) 

Z(I,J)=0. 
11    CONTINUE 
10    CONTINUE 

CALL  ZGRID(Z.IROW.ICOL,NY,NX,XLOW,YLOW,XHIGH,YHIGH,X,Y,ZZ, 

1NPTS,CAY,NRNG,ZPIJ,KNXT) 

XFACT=(XHIGH-XLOW)/(XHPLOT-XLPLOT) 

YFACT=(YHIGH-YLOW)/(YHPLOT-YLPLOT) 

CALL  OFFSET(XLOW,XFACT.YLOW,YFACT) 

CALL  PL0T(XL0W,YL0W,13) 

CALL  PL0T(XHIGH.YL0W,12) 

CALL  PL0T(XHIGH,YHIGH,12) 

CALL  PL0T(XL0W,YHIGH,12) 
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CALL  PL0T(XL0W,YL0W.12) 

TICKY=5. 

TICI0(=5. 

NYP=NY-t-l 

NXP=NX+1 

DELX=(XHIGH-XLOW)/FLOAT(NX) 

DELY=(YHIGH-YLOW)/FLOAT(NY) 

XG0=XG(NTR)/9. 

CALL  STAXIS(.07..1,.01..03.2) 

CALL  AXIS (0..0.,' METERS  DOWNSTREAM  FROM  SOURCE ',-29.9.. 0  0  XGO^ 
CALL  PL0T(XL0W.YL0W.13)  ^3.^.  ,u.  .u.  ,At.Uj 

CALL  STDASH(.05,.05) 
DO  709  1=1,9 
YY=FL0AT(I)*.3 
CALL  PL0T(-.1.YY,3) 
CALL  PL0TD(9.+.1,YY,2) 
709   CONTINUE 

CALL  PL0T(XHIGH.YL0W,13) 
DO  711  1=1.11 
YYY=FL0AT(I-1)/10. 
YY=YYY*3.-.ll 

CALL  NUMBER(-.18,YY,.07,YYY.0.,1) 
711   CONTINUE  ^ 

DO  6067  1=1, NTR 
CALL  PL0T(XG(I).YL0W,13) 
TICKT=13. 

CALL  PL0T(XG(I),YHIGH+TICKT,12) 
CALL  WHERE(XL0C,YL0C,FUD) 
TNUM=FL0AT(I) 

CALL  NUMBER(XLOC-.03.YLOC+.05,.075,TNUM,0.,-l) 
6067  CONTINUE  •   »  ^; 

CALL  PL0T(-0. 5,-0.5,-3) 

CALL  SYMBOL(.25,1.5,.10,'FLOW  FRACTION' ,90. .13) 
CALL  PL0T(0.5,0.5,-3)  ' 

CALL  ZCSEG(Z.IROW.ICOL.NY,NX.XLPLOT,YLPLOT,XHPLOT  YHPLOT 
l.ZLEV.LDIG,LWGT.NLEV.HGT.NDIV,NARC) 
35    CONTINUE 
C 
C 

CALL  COLOR(0,IERR) 
IF(ICAL.EQ.1)THEN 
QPL0T=(QCP/QRS)*3. 
CALL  SYMBOL(O..QPLOT, .1,1,0..-!) 
ELSE 

QPL0T1=(QCP1/QRS)*3. 

QPL0T2=(QCP2/QRS)*3. 

CALL  SYMB0L(0.,QPL0Tl,.l,1.0.,-n 

CALL  SYMB0L(0..QPL0T2, .1,1,0..-!) 


ENDIF 


CALL  SYMB0L(2.75.3.6,.09.'TRANSECT  NUMBER ' ,0. .15) 
CALL  SYMBOL(1.0.6.0..!5.TITLE.0..80) 
CALL  SYMB0L(1.0.5.7,.12.'POLLUTANT:'.0..10) 
CALL  SYMB0L(5.0.5.7,.12.P0LLU.0..3O) 
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CALL  SYMB0L(1.0,5.5,. 

CALL  NUMBER(5.0,5.5,. 

CALL  SYMB0L(6.25,5.5. 

CALL  SYMB0L(1.0,5.3,. 

CALL  NUMBER(5.0,5.3,. 

CALL  SYMB0L(6.25,5.3, 

CALL  SYMB0L(1. 0.5.1,. 

CALL  NUMBER(5. 0,5.1,. 

CALL  SYMB0L(1.0,4.9,. 

CALL  NUMBER(5.0,4.9,. 

IF(ICAL.EQ.1)THEN 
CALL  SYMBOL (1.0. 4. 7 
CALL  NUMBER(5.0,4.7 
CALL  SYMB0L(6.0,4.7 
CALL  SYMB0L(1.0,4.5 
CALL  SYMB0L(5.20,4. 

ENDIF 

IF(ICAL.EQ.2)THEN 
CALL  SYMB0L(1.0,4.7 
CALL  NUMBER(5.0,4.7 
CALL  NUMBER(6.5,4.7 
CALL  SYMB0L(6.0,4.7 
CALL  SYMB0L(1.0.4.5 
P0UT=FL0AT(N0UT) 
CALL  NUMBER(5.0,4.5 
CALL  SYMB0L(1.0,4.3 
CALL  SYMB0L(5.0,4.3 

ENDIF 

RETURN 

END 


12, 'RIVER  FLOW  RATE  ABOVE  OUTFALL  :  ',0.,32) 

12,QRUP,0.,2) 

.12. 'CMS', 0., 3) 

12, 'EFFLUENT  FLOW  RATE  :  ',0.,21) 

12,QEFL,0.,2) 

.12   'CMS'   0.   3) 

i2,'EFFLUENT'c0NCENTRATI0N   :    ',0..25) 

12,CEFL,0..2) 

12. 'BACKGROUND  CONCENTRATION  :  ',0..27) 

12.CBK6,0.,2) 

,.12, 'PIPE  OUTFALL  LOCATION  :  ',0..24) 

..12.Y0UT,0.,2) 

,.12.'   METERS  FROM  BANK ',0., 18) 

,.12, 'OUTFALL  LOCATION  MARKED  AS  :  ',0.,29) 

5, .09, 1,0. ,-1) 


..12,'DIFFUSER  OUTFALL  LOCATION  :  ',0.,28) 

,.12,Y0UT1,0.,2) 

,.12,Y0UT2,0.,2) 

,.12, 'TO  METERS  FROM  BANK ',0. ,26) 

,.12, 'NUMBER  OF  DIFFUSER  PORTS  :  '.0..27) 


,.12,P0UT.0.,0) 

,.12, 'DIFFUSER  ENDS  MARKED  AS 

,.12,1.0..-1) 


,0..26) 


SUBROUTINE   TRANS  PLOTS  CROSS  STREAM  CONCENTRATION  PROFILES 

SUBROUTINE   TRANS(NTR) 
COMMON/AA/   JK 

COMMON/L/  XG(50).C(10. 50) .TITLE 
COMMON/JJ/  CSAV(50).CSC(4).D(50) 

COMMON/PASS/  XDX(10),BS(10).HS(10).US(10),QRS,QRUP,QEFL,CEFL.CBKG 
^QCP,P0LLU,ICAL,NOUT.QCPl.QCP2,Y0UTl,Y0UT2,Y0UT 
CHARACTER*40  SNAME 
CHARACTER*30  POLLU 
CHARACTER*80  TITLE 
NXL=50 

CALL  C0L0R(0,IERR) 

CALL  PL0T(1.3,2.4.-3) 

IF(JK.EQ.2)CALL  FACT0R(.75) 

IF(JK.EQ.1)CALL  FACT0R(.6) 

IF(JK.EQ.3)CALL  FACT0R(.75) 

CALL  SYMB0L(1. 0,6.0,. 15, TITLE. 0., 80) 

CALL  SYMB0L(1. 0,5. 7. .12. 'POLLUTANT: ',0., 10) 

CALL  SYMB0L(5.0,5.7..12,P0LLU,0..30) 

CALL  SYM80L(1.0,5.5,. 12, 'RIVER  FLOW  RATE  ABOVE  OUTFALL   :    '.0.,32) 

CALL  NUMBER(5.0.5.5,.12.QRUP,0.,2) 

CALL  SYMB0L(6.25,5.5,.12.'CMS'.0..3) 

CALL  SYMB0L(1. 0.5. 3, .12, 'EFFLUENT  FLOW  RATE   :    ',0.,21) 

CALL  NUMBER(5.0.5.3..12.QEFL,0.,2) 
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CALL  SYMB0L(6.25,5.3,.12.'CMS'.0.,3) 

CALL  SYMeOL(l. 0.5.1, .12. 'EFFLUENT  CONCENTRATION  :  '.0..25) 

CALL  NUMBER(5.0.5.1,.12,CEFL,0.,2) 

CALL  SYM80L(1.0,4.9,. 12, 'BACKGROUND  CONCENTRATION  :  '.0.  27) 

CALL  NUMBER(5.0,4.9,.12,CBKG,0.,2) 

IF(ICAL.EQ.1)THEN 

CALL  SYMB0L(1.0,4.7,.12,'PIPE  OUTFALL  LOCATION  :  ',0..24) 

CALL  NUM6ER(5.0,4.7,.12,Y0UT,0.,2) 

CALL  SYMB0L(6.0.4.7,.12,'   METERS  FROM  BANK', 0., 18) 

ENDIF  ' 

IF(ICAL.EQ.2)THEN 

CALL  SYM60L(1.0.4.7..12,'DIFFUSER  OUTFALL  LOCATION  :  '  0  28) 
CALL  NUMBER(5.0,4.7,.12,Y0UT1,0.,2)  " 

CALL  NUMBER(6.5,4.7,.12,Y0UT2,0..2) 

CALL  SYMB0L(6.0,4.7,.12.'T0  METERS  FROM  BANK'. 0.. 31) 

CALL  SYMBOL (1.0. 4. 5.. 12, 'NUMBER  OF  DIFFUSER  PORTS  :  '.0  27) 
POUT=FLOAT(NOUT)  *^"^' ) 

CALL  NUMBER(5.0.4.5,.12,POUT.O..O) 
ENDIF 

CALL  PL0T(-1. ,-2.4,-3) 
C 

IF(JK.EQ.2)CALL  FACT0R(.3) 

IF(JK.EQ.1)CALL  FACT0R(.2) 

IF(JK.EQ.3)CALL  FACT0R(.3) 

SCA=10. 

CMAXX=0. 

DO  600  1=1, NTR 

DO  601  J=l.ll 

IF(C(I.J).GT.CMAXX)CMAXX=C(I.J) 
601   CONTINUE 
600   CONTINUE 

CSC(l)=-l.*CMAXX/3. 

CSC(2)=CMAXX 

CALL  PL0T(0..0..-3) 

DO  700  1=1, NTR 

DIST=XG(I) 

IF(I.6T.1)CALL  PL0T(-X0R.-Y0R,-3) 

IF(I.GE.1.AND.I.LE.4)Y0R=12. 

IF(I.GE.5.AND.I.LE.8)Y0R=3. 

IF(I.EQ.1.0R.I.EQ.5)X0R=1. 

IF(I.EQ.2.0R.I.EQ.6)X0R=9. 

IF(I.EQ.3.0R.I.EQ.7)XOR=17. 

IF(I.EQ.4.0R.I.EQ.8)X0R=25. 

CALL  PL0T(X0R.Y0R.-3) 

DO  500  J=l,ll 

CSAV(J)=C(I,J) 

D(J)=FL0AT(J-1)*SCA 
500   CONTINUE 

CALL  SCALE(CSC,5.,2,1) 

CSAV(12)=CSC(3) 

CSAV(13)=CSC(4) 

CALL  SCALE(D,5..11.1) 

CALL  STAXIS(.2..3..0C5,.1,1) 

CALL  AXIS(0.,0., 'PERCENT  FLOW' ,-12.5. .0. ,0(12) .D( 13) ) 

CALL  AXIS(0.,0. .'CONCENTRATION'. 13. 5., 90. ,CSAV(12),CSAV(13)) 
C      CALL  C0L0R(4.IERR)  ^      " 

CALL  LINE(0.CSAV. 11. 1.0.0) 
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C      CALL  COLOR(0,IERR) 

CALL  SYMB0L(1.5,.25,.2,'TRANS.  #',0..8) 

TNUM=FLOAT(I) 

CALL  NUMBER(3.5..25..2,TNUM.0.,-1) 
700   CONTINUE 

CALL  PL0T(-X0R,-Y0R,-3) 

RETURN 

END 
C 

C     SUBROUTINE  THREED  PLOTS  THREE  DIMENSIONAL  MESH  PLOTS 
C 

SUBROUTINE  THREED(NTR) 

COHMON/AA/  JK 

COMMON/B/  CC(50,50) 

COMMON/L/  XG(50),C(10, 50), TITLE 

COMMON/Q/  CCP(50,50) 

COMMON/PASS/  XDX(10),BS(10),HS(10).US(10),QRS,QRUP,QEFL,CEFL,CBKG 
*,QCP,P0LLU,ICAL,N0UT,QCP1,QCP2.Y0UT1,Y0UT2,Y0UT 

REAL*8  CT,RQ(21).TMTEMP,TM(6,51),MASK(1000),VERTEX(16) 

CHARACTER*25  XTITLE,YTITLE,ZTITLE 

CHARACTER*30  SNAME,POLLU 

CHARACTER*80  TITLE 

INTEGER*2  NX 

DIMENSION  NTM(6) 
C 

CMAX=0. 

NX=50 

NXL=NX 

DO  4  1=1, NXL 

DO  3  J=l,ll 

IF(C(I,J).GT.CMAX)CMAX=C(I,J) 

3  CONTINUE 

4  CONTINUE 
IF(CC(20,1).GT.0.)THEN 

DO  553  1=1, NXL 

DO  554  J=l,ll 

JJ=12-J 

CCP(I,JJ)=CC(I.J) 
554     CONTINUE 
553     CONTINUE 

ISI  =  1 
ELSE 

DO  5540  1=1, NXL 

DO  5530  J=l,ll 

CCP(I,J)=CC(I,J) 
5530    CONTINUE 
5540    CONTINUE 

ISI=0 
ENDIF 
NXSIZE=50 
NYSIZE=50 
AZIMUT=330. 
ELEVAT=30. 
XLPLOT=0. 
YLPLOT=0. 
XHPL0T=10. 
YHPL0T=5. 
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IED6E=0 

IDIR=3 

IPR0J=0 

IFRAME=1 

ZL0W=1.0E35 

ICUT=0 

ITRIM=0 

CALL  SIMPLX 

IF(JK.EQ.2)CALL  FACT0R(.75) 

IF(JK.EQ.1)CALL  FACT0R(.6) 

IF(JK.EQ.3)CALL  FACT0R(.75) 

CALL  PL0T(3.0,3.3.-3) 

CALL  C0L0R(0,IERR) 


CALL  SYMBOL 
CALL  SYMBOL 
CALL  SYMBOL 
CALL  SYMBOL 
CALL  NUMBER 
CALL  SYMBOL 
CALL  SYMBOL 
CALL  NUMBER 
CALL  SYMBOL 
CALL  SYMBOL 
CALL  NUMBER 
CALL  SYMBOL 
CALL  NUMBER 
1F(ICAL.EQ 


1.0,6.0,.15,TITLE.0.,80) 

1.0,5.7,.12,'POLLUTANT: ',0.,10) 

5.0,5.7,.12.P0LLU,0.,30) 

1.0,5.5,. 12. 'RIVER  FLOW  RATE  ABOVE  OUTFALL  :  ',0.,32) 

5.0,5.5.. 12, QRUP.O. ,2) 

6.25.5.5,.12,'CMS',0.,3) 

1.0, 5. 3, .12, 'EFFLUENT  FLOW  RATE    :    ',0..21) 

5.0,5.3,. 12, QEFL,0. .2) 

6. 25, 5. 3, .12, 'CMS', 0., 3) 

1.0,5.1,.12, 'EFFLUENT  CONCENTRATION  :  ',0.,25) 

12,CEFL,0..2) 

12, 'BACKGROUND  CONCENTRATION  :  ',0.,27) 

12,CBKG,0.,2) 


5.0,5.1,, 
1.0,4.9,, 
5.0,4.9,, 
)THEN 


CALL  SYMB0L(1.0,4.7, 
CALL  NUMBER(5.0,4.7, 
CALL  SYMB0L(6.0,4.7, 

ENDIF 

IF(ICAL.EQ.2)THEN 
CALL  SYMB0L(1.0,4.7, 
CALL  NUMBER(5.0,4.7, 
CALL  NUMBER(6.5,4.7, 


12, 'PIPE  OUTFALL  LOCATION  :  ' ,0. ,24) 

12.Y0UT.0.,2) 

12,'   METERS  FROM  BANK'. 0., 18) 


12,'DIFFUSER  OUTFALL  LOCATION  :  ',0.,28) 
12,Y0UT1.0.,2) 
12,Y0UT2.0..2) 

CALL  SYMB0L(6.0,4.7,.12,'T0       METERS  FROM  BANK', 0., 26) 
CALL  SYMB0L(1.0,4.5,.12,'NUMBER  OF  DIFFUSER  PORTS  :  ',0.,27) 
POUT=FLOAT(NOUT) 

CALL  NUMBER(5.0,4.5,.12,P0UT,0..0) 
ENDIF 
CALL  PL0T(-2.,-2.5,-3) 

NYZ=11 

CALL  MESH(CCP,NXSIZE,NYSIZE,NX,NYZ,AZIMUT,ELEVAT,XLPLOT,YLPLOT. 

1  XHPLOT,YHPLOT,IEDGE,IDIR,IPROJ,IFRAME,ZLOW.ICUT.ITRIM, 

2  MASK, VERTEX) 

SET  UP  THE  AXIS  PARAMETERS  AND  PLOT  THE  AXES 

RNTP=FLOAT(NX) 

RNY=11. 

XMIN=0. 

XMAX=1000. 

CALL  P3D2D(1.,1.,0.,X0RG,Y0RG) 

CALL   P3D2D(RNTP.1.,0..XX,YX) 

CALL   P3D2D(1.,RNY,0.,XY,YY) 

CALL  P3D2D(1.,RNY.CMAX.XZY,YZY) 


NOV.  14.  1986   16:00       

CONMIX  __._ 


XrLTA;rX^X-XMIN)/XLEN 
YDELTA=1./VLEN 

^f,^[r.rs(A«CE  DO-N  STREAM. 

mni' FRACTIONAL  DISCHARGE' 

NVT=20 
ZTITLE='C0NCENTRAT10N' 

NZT=13 


Jfa  STAXIS(a25..p,.0625..0625  1)  „^ 

CALL  AXISCXORG  YORG.XTITLE   -NXT  XL^^^^^^ 

^C^L  fxl^sSrv;r.;rTrTU;:ZT!ilEK.ZAHGLE.O.  .ZDELTA) 

RETURN 
END 
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MISSISSIPPI   RIVER 
AhMONIA 

TRANSECT     1 
PARAMETER   1: 
QRIVER=  2374.020 
QEFL     =        76.000 
UPSTREAM  FLUX= 


.00000      15.00000 


.00000 


.00000 


.00 


200.0  METERS   FROM  OUTFALL 
AMMONIA 

BACKGROUND  CONC.=  .000 

EFFLUENT  CONCN.=       15.000 
,00       EFFLUENT  FLUX=   1140.00       TOTAL  FLUX   =    1140 


VELOCITIES  SIMULATED  FROM  RESISTANCE   EQN. 

VEL  CONC  SUMA       SUMQ       SUMF       Y/B       QY/QT     C/CAVG 


.00 

.00 

.00 

8 

.00 

.00     .00 

.00 

.000 

.000 

16 

.660  1 

16.67 

1.67 

.55 

6 

.00 

13.92    3.83 

26.80 

.030 

.002 

12 

.495  1 

44.44 

8.89 

1.69 

1 

.00 

160.55  167.80 

600.71 

.081 

.071 

2 

.082 

83.33 

7.78 

1.54 

.50 

484.69  691.13 

993.21 

.151 

.291 

1 

.041 

132.22 

6.11 

1.31 

.00 

824.23  1175.69 

1114.35 

.240 

.495 

.000 

150.00 

5.56 

1.23 

.00 

927.98  1307.62 

1114.35 

.272 

.551 

.000 

177.78 

4.44 

1.06 

.00 

1066.88  1466.71 

1114.35 

.323 

.618 

.000 

222.22 

3.00 

.81 

.00 

1232.20  1621.59 

1114.35 

.403 

.683 

.000 

333.33 

2.22 

.67 

.00 

1522.19  1836.22 

1114.35 

.605 

.773 

.000 

444.44 

4.44 

1.06 

.00 

1892.19  2155.33 

1114.35 

.806 

.908 

.000 

551.11 

1.11 

.42 

.00 

2188.20  2374.02 

1114.35 

1.000  1 

.000 

.000 

TRANSECT  1 

VARIANCE  FROM  DIFFERENT 

METHODS: 

PARAMETER 

VCMAX 

VCN       VUF 

VCQ 

VPQ 

162.04 


1409.69 


3182.71     163469.30 


3087.82 


NONDIMENSIONAL  VARIANCE 
PARAMETER   VCN/BB    VCN/HH 


1 
AMMONIA 
1 


X/B=  .36  X/H=  50.4 

VUF/BB  VUF/HH  VCQ/QQ 


.0046  89.419 

.00000     15.00000 


.0105         201.883 
.00000         .00000 


.0290 
.00000 


.  00000 


TRANSECT  2 
PARAMETER  1: 
QRIVER=  2374.020 
QEFL  =   76.000 
UPSTREAM  FLUX= 


2000.0  METERS  FROM  OUTFALL 
AMMONIA 

BACKGROUND  CONC.= 

EFFLUENT  CONCN.= 
.00   EFFLUENT  FLUX= 


.000 

15.000 

1140.00 


TOTAL  FLUX 


1140 


VELOCITIES  SIMULATED  FROM  RESISTANCE  EQN. 


VEL   CONC    SUMA   SUMQ   SUMF   Y/B   QY/QT  C/CAVG  C/ 


.00  .00 

.00 

3.00 

.00 

.00 

.00 

.000 

.000 

6.247 

6. 

16.67  2.22 

.70 

2.70 

18.50 

6.51 

18.57 

.036 

.003 

5.623 

5. 

100.00  3.56 

.97 

2.00 

259.33 

207.63 

491.18 

.216 

.087 

4.165 

4. 

168.89  5.00 

1.21 

1.00 

554.18 

528.89 

973.08 

.365 

.223 

2.082 

2. 

188.89  5.38 

1.27 

.50 

657.98 

657.97 

1069.89 

.408 

.277 

1.041 

1. 

205.56  5.56 

1.30 

.00 

749.16 

775.44 

1099.26 

.444 

.327 

.000 

228.89  5.56 

1.30 

.00 

878.88 

944.39 

1099.26 

.494 

.398 

.000 

300.00  7.78 

1.63 

.00 

1353.18 

1640.13 

1099.26 

.647 

.691 

.000 

341.11  5.00 

1.21 

.00 

1615.87 

2013.72 

1099.26 

.736 

.848 

.000 
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441   11  2  28         72         .00  1979.87  2364.94     1099.26     .952     .996       .000 
463!33     .00       .00         .00  2005.20  2374.02     1099.26  1.000   1.000       .000 


AVG.  CONC.  JUST  BELOW  OUTFALL.  CAVG=  .480 
AVG.  CONC.  AT  THE  TRANSECT,  CATRN  =  .463 
MEAN  DEPTH=  4.328     MEAN  VELOCITY=  1.184 


TOTAL  FLUX  AT  TRANSE 
SHAPE-VELOCITY  FACTOR= 


TRANSECT  2 
PARAMETER    VCMAX 


:  VARIANCE  FROM  DIFFERENT  METHODS: 

VCN       VUF       VCQ       VPQ 


2366.60    8691.42   13245.78  109118.10   21367.13 


NONDIMENSIONAL  VARIANCE     X/B=    4.32     X/H= 
PARAMETER    VCN/BB     VCN/HH     VUF/BB     VUF/HH 


1 
AMMONIA 
1 


.0405    464.039 
.00000  15.00000 


.0617         707.199 
.00000         .00000 


462.1 

VCQ/QQ 

.0194 
.00000         .0000 


TRANSECT      3 
PARAMETER   1: 
QRIVER=   2374.020 
QEFL     =       76.000 
UPSTREAM  FLUX= 


4500.0  METERS  FROM  OUTFALL 
AMMONIA 
BACKGROUND  CONC.=  .000 

EFFLUENT  C0NCN.=       15.000 
.00       EFFLUENT  FLUX=      1140.00 


TOTAL  FLUX  = 


VELOCITIES  SIMULATED  FROM  RESISTANCE  EQN. 


Z   VEL  CONC    SUMA   SUMQ   SUMF   Y/B   QY/QT  C/CAVG  C/C 


.00 

.00 

.00 

2 

.90 

.00 

.00 

.00 

.000 

.000 

6 

.039 

9.9 

61.11 

1.67 

1.24 

2 

.40 

51.03 

31.55 

83.60 

.090 

.013 

4 

.998 

8.2 

111.11 

2.22 

1.50 

2 

.00 

148.28 

164.43 

375.94 

.164 

.069 

4 

.165 

6.8 

127.78 

2.22 

1.50 

1 

.70 

185.28 

219.80 

478.38 

.189 

.093 

3 

.540 

5.8 

142.22 

2.22 

1.50 

1 

.10 

217.34 

267.77 

545.53 

.210 

.113 

2 

.291 

3.7 

205.56 

1.33 

1.06 

.50 

329.77 

411.56 

660.56 

.303 

.173 

1 

.041 

1.7 

283.33 

1.67 

1.24 

.00 

446.42 

545.60 

694.07 

.418 

.230 

.000 

.0 

311.11 

2.67 

1.69 

.00 

506.71 

633.90 

694.07 

.459 

.267 

.000 

.0 

381.11 

2.67 

1.69 

.00 

693.61 

950.37 

694.07 

.562 

.400 

.000 

.0 

433.33 

3.89 

2.18 

.00 

864.89 

1281.99 

694.07 

.639 

.540 

.000 

.0 

477.78 

3.33 

1.96 

.00 

1025.35 

1614.33 

694.07 

.705 

.680 

.000 

.0 

527,78 

3.67 

2.10 

.00 

1200.35 

1969.48 

694.07 

.779 

.830 

.000 

.0 

550.00 

2.28 

1.52 

.00 

1266.46 

2089.09 

694.07 

.811 

.880 

.000 

.0 

577.78 

1.89 

1.34 

.00 

1324.38 

2172.11 

694.07 

.852 

.915 

.000 

.0 

662.22 

1.67 

1.24 

.00 

1474.68 

2365.99 

694.07 

.977 

.997 

.000 

.0 

677.78 

.00 

.00 

.00 

1487.67 

2374.02 

694.07 

1.000 

1.000 

.000 

.0 

AVG.  CONC.  JUST  BELOW  OUTFALL,  CAVG=    .480 
AVG.  CONC.  AT  THE  TRANSECT,  CATRN  =  .292 
MEAN  DEPTH=  2.195    MEAN  VELOCITY=  1.596 


TOTAL  FLUX  AT  TRANSECT= 
SHAPE -VELOCITY  FACTOR= 


TRANSECT  3       :  VARIANCE  FROM  DIFFERENT  METHODS: 
PARAMETER    VCMAX      VCN       VUF       VCQ       VPQ 

1     2924.33   11368.83   15000.90   45450.23    9116.05 
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NONDIMENSIONAL  VARIANCE     X/B=    6.64 
PARAMETER    VCN/BB     VCN/HH     VUF/BB 


1 
AMMONIA 
1 


.0247   2359.812 
.00000  15.00000 


.0327 
.00000 


X/H  = 

VUF/HH 

3113.715 
.00000 


2050.2 

VCQ/QQ 

.0081 
.00000    .0000 


TRANSECT  4 
PARAMETER 
QRIVER=  2374.020 
QEFL  =   76.000 
UPSTREAM  FLUX= 


10150.0  METERS  FROM  OUTFALL 
1:  AMMONIA 

BACKGROUND  CONC.=  .000 
EFFLUENT  C0NCN.=  15.000 
.00   EFFLUENT  FLUX=  1140.00 


TOTAL  FLUX  =   1 


VELOCITIES  SIMULATED  FROM  RESISTANCE  EQN, 
Z   VEL   CONC   SUMA   SUMQ   SUMF   Y/B   QY/QT  C/CAVG  C/CA 


.00  .00 

.00 

2.30 

.00 

.00 

.00 

.000 

.000 

4.790 

6.5 

26.67  2.89 

.66 

2.00 

38.54 

12.69 

27.29 

.075 

.005 

4.165 

5.6 

62.22  3.56 

.76 

1.60 

153.19 

93.88 

173.43 

.175 

.040 

3.332 

4.5 

94.44  5.56 

1.02 

1.10 

300.11 

224.55 

349.83 

.266 

.095 

2.291 

3.1 

127.78  7.22 

1.22 

.90 

513.15 

462.92 

588.21 

.359 

.195 

1.874 

2.5 

150.00  8.11 

1.32 

.20 

683.47 

678.52 

706.78 

.422 

.286 

.416 

.5 

222.22  7.78 

1.28 

.10 

1257.26 

1422.77 

818.42 

.625 

.599 

.208 

.2 

255.56  8.22 

1.33 

.00 

1523.98 

1770.31 

835.80 

.719 

.746 

.000 

.0 

283.33  8.44 

1.35 

.00 

1755.30 

2080.03 

835.80 

.797 

.876 

.000 

.0 

344.44  1.89 

.50 

.00 

2070.93 

2371.42 

835.80 

.969 

.999 

.000 

.0 

355.56  .00 

.00 

.00 

2081.44 

2374.02 

835.80  1 

.000  1 

.000 

.000 

.0 

AVG.  CONC.  JUST  BELOW  OUTFALL.  CAVG=    .480 
AVG.  CONC.  AT  THE  TRANSECT,  CATRN  =   .352 
MEAN  DEPTH=  5.854     MEAN  VELOCITY=  1.141 


TOTAL  FLUX  AT  TRANSECT=  83 
SHAPE-VELOCITY  FACTOR= 


TRANSECT  4       :  VARIANCE  FROM  DIFFERENT  METHODS: 
PARAMETER    VCMAX       VCN       VUF        VCQ        VPQ 

1      1494.55    6824.84   13716.91  286194.10   21015.44 


NONDIMENSIONAL  VARIANCE     X/B=   28.55 
PARAMETER    VCN/BB     VCN/HH     VUF/BB 


1 

AMMONIA 
1 


.0540         199.154 
.00000     15.00000 


X/H  = 
VUF/HH 


1733.9 
VCQ/QQ 


.1085         400.270  .0508 

.00000         .00000         .00000         .00000 


TRANSECT  5         17450.0  METERS  FROM  OUTFALL 

PARAMETER  1:  AMMONIA 
QRIVER=  2374.020     BACKGROUND  CONC.=    .000 
QEFL  =   76.000     EFFLUENT  CONCN.=   15.000 
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UPSTREAM  FLUX=     .00  EFFLUENT  FLUX=  1140.00 
VELOCITIES  SIMULATED  FROM  RESISTANCE  EQN, 


TOTAL  FLUX 


1140. 


Y   Z   VEL   CONC    SUMA  SUMQ   SUMF   Y/B   QY/QT  C/CAVG  C/CAT 


.00 

.00 

.00 

1.80 

.00 

.00 

.00 

.000 

.000 

3.748 

5.16 

77.78 

3.56 

.93 

1.30 

138.45 

64.49 

99.97 

.104 

.027 

2.707 

3.72 

144.44 

3.40 

.90 

1.00 

370.43 

277.34 

344.74 

.193 

.117 

2.082 

2.86 

233.33 

3.35 

.89 

.80 

670.43 

547.03 

587.46 

.311 

.230 

1.666 

2.29 

250.00 

3.33 

.89 

.70 

726.11 

596.74 

624.74 

.333 

.251 

1.458 

2.00 

303.33 

3.11 

.85 

.60 

897.83 

746.30 

721.96 

.404 

.314 

1.249 

1.72 

338.89 

2.22 

.68 

.30 

992.60 

818.80 

754.58 

.452 

.345 

.625 

.86 

383.33 

2.79 

.79 

.20 

1103.92 

900.64 

775.04 

.511 

.379 

.416 

.57 

483.33 

3.67 

.95 

.10 

1426.92 

1182.00 

817.25 

.644 

.498 

.208 

.28 

550.00 

3.33 

.89 

.00 

1660.26 

1396.89 

827.99 

.733 

.588 

.000 

.00 

611.11 

5.00  1 

.17 

.00 

1914.79 

1659.14 

827.99 

.815 

.699 

.000 

.00 

727.78 

5.00  1 

.17 

.00 

2498.14 

2341.53 

827.99 

.970 

.986 

.000 

.00 

750.00 

.00 

.00 

.00 

2553.69 

2374.02 

827.99 

1.000  1 

.000 

.000 

.00 

AVG. 

CONC. 

JUST 

BELOW 

OUTFALL, 

.  CAVG= 

.480 

AVG. 

CONC. 

AT  THE  TRANSECT,  CATRN  = 

.349   TOTAL  FLUX  AT 

TRANSECT= 

MEAN 

DEPTH= 

3.405 

MEAN  VELOCITY= 

.930 

SHAPE- 

VELOCITY  FACT0R= 

TRANSECT 

5 

:  VARIANCE  FROM 

DIFFERENT  METHODS: 

PARAMETER 

VCMAX 

VCh 

1 

VUF 

VCQ 

VPQ 

1 

6718.22 

38319.33   46227.59  245477.00 

33674.34 

NONDIMENSIONAL  VARIANCE     X/B=   23.27     X/H=   5124.9 
PARAMETER    VCN/BB     VCN/HH     VUF/BB     VUF/HH     VCQ/QQ 
1        .0681   3305.256     .0822   3987.387      .0436 
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ENTER  TITLE  OF  STUDY 
MISSISSIPPI  RIVER 

ENTER  QRS.BPWR.HPWR.UPWR.THETA, TEMPS, RBK 


2374.02  .05         .50 

ENTER  DESIGN  CASE: 
2298.02         76.000  15.00 

ENTER  NTR,NYZ,QCP 

5    10   .0 
ENTER  5  VALUES  OF  X,BS,HS,US 


45       1.030       20.00 
QRUP,QEFL,CEFL,CBKG,TMP 


,00 


.00   22.90 


1         200.00 

551.11 

3.97 

1 

0800 

2       2000.00 

463.33 

4.33 

1 

1800 

3       4500.00 

677.78 

2.19 

1 

6000 

4     10150.00 

355.56 

5.85 

1 

1400 

5     17450.00 

750.00 

3.40 

9300 

ENTER     5  VALUES 

OF  DECAY 

1    .0000231 

2   .0000231 

3   .0000231 

4   .0000231 

5   .0000231 

ENTER  VALUES  OF 

BETA 

TRANSECT        1    : 

.0040000 

TRANSECT        2 

.0017000 

TRANSECT       3 

.0007800 

TRANSECT        4 

.0010000 

TRANSECT       5 

.0007000 

ARE  YOU  CONSIDERING  UN-IONIZED  AMMONIA:  YES=1   N0=0 
ENTER  THE  PH  OF  UN-IONIZED  AMMONIA 
THE  PH  OF  UN-IONIZED  AMMONIA  IS   8.3 


TRANSECT: 

1 

BETA= 

.004000     RKS=      . 

000023 

X 

BM 

HW 

UW 

200.000       551. 

110 

3.970 

1.080 

QY           C(X, 

QY) 

CUI 

C/CA 

QY/QT 

OIL 

FAC 

1 

.00 

7.0777 

.6304         14 

.7392 

.0000 

2.12 

2 

237.40 

1.2646 

.1126           2 

.6334 

.1000 

11.86 

3 

474.80 

.0072 

.0006 

.0150 

.2000 

2079.68 

4 

712.21 

.0000 

.0000 

.0000 

.3000 

-1.00 

5 

949.61 

.0000 

.0000 

.0000 

.4000 

-1.00 

6 

1187.01 

.0000 

.0000 

.0000 

.5000 

-1.00 

7 

1424.41 

.0000 

.0000 

.0000 

.6000 

-1.00 

8 

1661.81 

.0000 

.0000 

.0000 

.7000 

-1.00 

9 

1899.22 

.0000 

.0000 

.0000 

.8000 

-1.00 

10 

2136.62 

.0000 

.0000 

.0000 

.9000 

-1.00 

11 

2374.02 

.0000 

.0000 

.0000 

1.0000 

-1.00 

TRANSECT: 

2 

BETA= 

.001700     RKS=      . 

000023 

X 

BW 

HU 

UW 

2000.000       511. 

609 

4.131 

1.123 

QY           C(X 

QY) 

CUI 

C/CA 

QY/QT 

DIL 

FAC 

1 

.00 

3.1833 

.2835           6 

.6291 

.0000 

4.71 

2 

237,40 

2.1853 

.1946           4 

.5507 

.1000 

6.86 
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3 

474.80 

.7069 

.0630 

1.4722 

.2000 

21.22 

4 

712.21 

.1078 

.0096 

.2244 

.3000 

139.18 

5 

949.61 

.0077 

.0007 

.0161 

.4000 

1937.32 

6 

1187.01 

.0003 

.0000 

.0005 

.5000 

57222.88 

7 

1424.41 

.0000 

.0000 

.0000 

.6000 

-1.00 

8 

1661.81 

.0000 

.0000 

.0000 

.7000 

-1.00 

9 

1899.22 

.0000 

.0000 

.0000 

.8000 

-1.00 

10 

2136.62 

.0000 

.0000 

.0000 

.9000 

-1.00 

11 

2374.02 

.0000 

.0000 

.0000 

1.0000 

-1.00 

TRANSECT: 

3 

BETA=     , 

.000780     RKS= 

.000023 

X 

BW 

HW 

uw 

4500.000       544.: 

357          ; 

3.624 

1.204 

I 

QY           C(X, 

QY) 

CUI 

C/CA 

QY/QT 

OIL 

FAC 

1 

.00 

3.1071 

.2767 

6.4705 

.0000 

4.83 

2 

237.40 

2.1085 

.1878 

4.3909 

.1000 

7.11 

3 

474.80 

.6589 

.0587 

1.3721 

.2000 

22.77 

4 

712.21 

.0948 

.0084 

.1975 

.3000 

158.19 

5 

949.61 

.0063 

.0006 

.0131 

.4000 

2387.14 

6 

1187.01 

.0002 

.0000 

.0004 

.5000 

78222.83 

7 

1424.41 

.0000 

.0000 

.0000 

.6000 

-1.00 

8 

1661.81 

.0000 

.0000 

.0000 

.7000 

-1.00 

9 

1899.22 

.0000 

.0000 

.0000 

.8000 

-1.00 

10 

2136.62 

.0000 

.0000 

.0000 

.9000 

-1.00 

11 

2374.02 

.0000 

.0000 

.0000 

1.0000 

-1.00 

TRANSECT: 

4 

BETA  = 

.001000     RKS= 

.000023 

X 

BW 

HW 

uw 

10150.000       528. < 

)45 

3.839 

1.169 

QY           C(X, 

QY) 

CUI 

C/CA 

QY/QT 

OIL 

FAC 

1 

.00 

1 . 5899 

.1416 

3.3109 

.0000 

9.43 

2 

237.40 

1.3957 

.1243 

2.9065 

.1000 

10.75 

3 

474.80 

.9442 

.0841 

1.9662 

.2000 

15.89 

4 

712.21 

.4922 

.0438 

1.0250 

.3000 

30.48 

5 

949.61 

.1977 

.0176 

.4118 

.4000 

75.86 

6 

1187.01 

.0612 

.0055 

.1275 

.5000 

245.04 

7 

1424.41 

.0146 

.0013 

.0304 

.6000 

1027.14 

8 

1661.81 

.0027 

.0002 

.0056 

.7000 

5586.98 

9 

1899.22 

.0004 

.0000 

.0008 

.8000 

39434.46 

10 

2136.62 

.0000 

.0000 

.0001 

.9000 

359252.00 

11 

2374.02 

.0000 

.0000 

.0000 

1.0000 

-1.00 

TRANSECT: 

5 

BETA= 

.000700     RKS= 

.000023 

X 

BW 

HW 

UW 

17450.000       538.916 

4.176 

1.055 

QY           C(X, 

QY) 

CUI 

C/CA 

QY/QT 

OIL 

FAC 

1 

.00 

1.2007 

.1069 

2.5003 

.0000 

12.49 

2 

237.40 

1.0753 

.0958 

2.2392 

.1000 

13.95 

3 

474.80 

.7723 

.0688 

1.6084 

.2000 

19.42 

4 

712.21 

.4449 

.0396 

.9266 

.3000 

33.71 

5 

949.61 

.2056 

.0183 

.4281 

.4000 

72.96 

6 

1187.01 

.0762 

.0068 

.1587 

.5000 

196.89 

7 

1424.41 

.0226 

.0020 

.0472 

.6000 

662.43 

a 

1661.81 

.0054 

.0005 

.0112 

.7000 

2778.83 
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9  1899.22  .0010  .0001  .0021  .8000         14532  04 

10  2136.62  .0002  .0000  .0003  .9000         93644  91 

11  2374.02  .0000  .0000  .0001  1.0000       385323.80 
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ENTER  TITLE  OF  STUDY 
MISSISSIPPI   RIVER 

ENTER  QRS.BPWR.HPWR.UPWR.THETA, TEMPS, RBK 
3100.000         .050         .500       .450         .0300     20.0000 
ENTER  DESIGN  CASE:   QRUP.QEFL.CErL.CBKG.TMP 
3000.000     100.00  20.00         .00       22.90 

ENTER  NTR,NYZ,N0UT,QCP1,QCP2 


.00 


10 


10     1117.73     2188.74 


ENTER  NTR  VALUES  OF  X.BS.HS.US 


200.00         551.11 

3.97 

1 

.080 

2000.00         463.33 

4.33 

1 

.180 

4500.00         677.78 

2.19 

1 

.600 

10150.00         355.56 

5.85 

1 

.140 

17450.00         750.00 

3.40 

.930 

.00002310 

.00002310 

.00002310 

.00002310 

.00002310 

1        .0040000 

2        .0017000 

3        .0007800 

4        .0010000 

5        .0007000 

UN-IONIZED  AMMONIA: 

ENTER 

1 

FOR 

YES, 

AND  ENTER 

PH  VALUE   ON  THE   SAME   LINE 

1 

8.3 

0  FOR  NO 


MISSISSIPPI  RIVER 

UPSTREAM  FL0W=3000.000   BACKGROUND  CONC.=    .000  DESIGN  TEMP=  22.9 

EFFLUENT  FLOW=  100.000   EFFLUENT  CONC.=   20.00 

DIFFUSER  OUTFALL  LOCATED  BETWEEN1117. 73  AND  2188.74 


TRANSECT:   1  BETA= 
X        BW 
200.000   551.110 


.004000  RKS=  .000023 

HW       UW 
3.970     1.080 


QY 


C(X,QY) 


CUI 


C/CA 


QY/QT  OIL  FAC 


I           .00 

.0000 

.0 

.0000 

.0     -1.0000 

2     310.00 

.0000 

.0 

.0000 

.1     -1.0000 

3     620.00 

.0027 

.0 

.0042 

.27426.6480 

4     930.00 

.2437 

.0 

.3778 

.3     82.0535 

5   1240.00 

1.4340 

.1 

2.2227 

.4     13.9469 

6  1550.00 

1.8583 

.2 

2.8803 

.5     10.7628 

7   1860.00 

1.8216 

.2 

2.8235 

.6     10.9794 

8  2170.00 

1.0171 

.1 

1.5765 

.7     19.6637 

9  2480.00 

.0758 

.0 

.1175 

.8  263.7737 

10  2790.00 

.0003 

.0 

.0005 

9**  *-*♦♦♦** 

11  3100.00 

.0000 

.0 

.0000 

i!o    -1.0000 

TRANSECT: 

2     BETA= 

.001700 

RKS=     .000023 

X 

BW 

HW 

UW 
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2000.000   511.609     4.131 
K    QY     C(X,QY)     CUI 


1     .00 

.0033 

2  310.00 

.0223 

3  620.00 

.1529 

4  930.00 

.5593 

5  1240.00 

1.1763 

6  1550.00 

1.5869 

7  1860.00 

1.4981 

8  2170.00 

.9697 

9  2480.00 

.3875 

10  2790.00 

.0870 

11  3100.00 

.0201 

1 

.467 

C/CA 

QY/QT  OIL  FAC 

0 

.0051 

.06073.3723 

0 

.0346 

.1  896.9067 

0 

.2369 

.2  130.8418 

0 

.8668 

.3  35.7622 

1 

1.8232 

.4  17.0029 

1 

2.4598 

.5  12.6029 

1 

2.3221 

.6  13.3499 

1 

1.5031 

.7  20.6245 

0 

.6006 

.8  51.6177 

0 

.1348 

.9  229.9990 

0 

.0312 

1.0  993.4742 

TRANSECT:  3  BETA=  .000780  RKS=  .000023 
X        BW        HW       UM 
4500.000   544.357     3.624     1.572 


K    QY 

C(X,QY)     CUI 

C/CA 

QY/QT  OIL  FAC 

1     .00 

.0028 

.0 

.0043 

.07148.1117 

2  310.00 

.0204 

.0 

.0316 

.1  981.9149 

3  620.00 

.1470 

.0 

.2278 

.2  136.0911 

4  930.00 

.5541 

.0 

.8589 

.3  36.0915 

5  1240.00 

1.1808 

.1 

1.8302 

.4  16.9382 

6  1550.00 

1.5975 

.1 

2.4761 

.5  12.5198 

7  1860.00 

1.5076 

.1 

2.3367 

.6  13.2665 

8  2170.00 

.9707 

.1 

1.5046 

.7  20.6032 

9  2480.00 

.3809 

.0 

.5904 

.8  52.5084 

10  2790.00 

.0823 

.0 

.1276 

.9  242.9518 

11  3100.00 

.0180 

.0 

.0279 

1.01111.1721 

TRANSECT:  4  BETA=  .001000  RKS=  .000023 
X        BU        HW       UW 
10150.000   528.945     3.839     1.527 


K    QY 

C(X,QY) 

CUI 

C/CA 

QY/QT 

OIL  FAC 

1     .00 

.1221 

.0 

.1892 

.0 

163.8056 

2  310.00 

.1869 

0 

.2898 

.1 

106.9853 

3  620.00 

.3799 

0 

.5888 

.2 

52.6453 

4  930.00 

.6720 

1.0417 

.3 

29.7598 

5  1240.00 

.9724 

1.5072 

.4 

20.5681 

6  1550.00 

1.1490 

1.7809 

.5 

17.4070 

7  1860.00 

1.1115 

1.7229 

.6 

17.9929 

8  2170.00 

.8813 

1.3660 

.7 

22.6935 

9  2480.00 

.5773 

.8948 

.8 

34.6441 

10  2790.00 

.3365 

0 

.5216 

.9 

59.4311 

11  3100.00 

.2472 

0 

.3832 

1.0 

80.9018 

TRANSECT:  5  BETA= 
X        BW 


.000700  RKS=  .000023 
HW       UW 
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17450.000       538.916  4.176 

K         QY  C(X,QY)  CUI 


1           .00 

.1670 

2     310.00 

.2306 

3     620.00 

.4125 

4     930.00 

.6736 

5   1240.00 

.9306 

6  1550.00 

1.0780 

7   1860.00 

1.0475 

8  2170.00 

.8562 

9  2480.00 

.5978 

10  2790.00 

.3873 

11  3100.00 

.3076 

1.377 

C/CA 

QY/QT 

OIL  FAC 

.0 

.2589 

.0 

119.7269 

.0 

.3574 

.1 

86.7255 

.0 

.6393 

.2 

48.4879 

1.0442 

.3 

29.6890 

1.4425 

.4 

21.4904 

1.6709 

.5 

18.5530 

1.6237 

.6 

19.0926 

1.3272 

.7 

23.3582 

.9266 

.8 

33.4566 

!o 

.6003 

.9 

51.6429 

.0 

.4768 

1.0 

65.0209 
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ENTER  TITLE  OF  STUDY 
MISSISSIPPI  RIVER  -  TEST  DATA 

ENTER  POLLUTANT  NAME 
CHLORINE 

ENTER  QRS.BPWR.HPWR.UPWR, TEMPS, NTR 


2374.0 


1 


,50 


.45  22.9 


ENTER  NTR  VALUES  OF   X,BS,HS,US 


200.0 

551.1 

3.97 

1.080 

2000.0 

463.3 

4.33 

1.180 

4500.0 

677.8 

2.19 

1.600 

10150.0 

355.6 

5.85 

1.140 

17450.0 

750.0 

3.40 

.930 

ENTER  NTR  VALUES  OF  BETA 
.00400 
.00170 
.00078 
.00100 
.00070 
ENTER  MQ  &  QRUP  VALUES 
QRUP(        1)=      1000.00 
ORUP(       2)=     3000.00 
ENTER  MT  &  TEMP  VALUES 

TMP(      1)=      22.0 
ENTER  MF  &  QELF  VALUES 

QEFL(    1)=   75.00 
UN-IONIZED  AMMONIA:    ENTER   1   FOR  YES 
0  FOR  NO 
AMONIA=0. 
ENTER  QCP,CEFL,CBKG,CS,THETA,RBK,XWCP 

.00  20.00  .00  .02 

ENTER  NTR  VALUES  OF  RKS 
RKS(  1)=  .000023 
RKS(  2)=  .000023 
RK:S(  3)=  .000023 
RKS(  4)=  .000023 
RKS(   5)=   .000023 


1.10 


.00  20000.0 


PREDICTIONS  OF  RUNS  FOR  MANAGEMENT  OPTIONS 
MISSISSIPPI  RIVER  -  TEST  DATA 


*  *  RUN  NO.: 

1 

QEFL=     75.000 

QRUP= 

1000.000 

TEMPR=22 

CEFL=       20.00 

CS= 

.02 

,0  PH=  7.0 


EY 


.0 


.1 


.2 


.3 


.4 


200.0         1.602       20.144 

3.848 

.027 

.000 

.000 

2000.0            .625         8.991 

6.263 

2.117 

.347 

.028 

4500.0            .569         8.707 

5.998 

1.961 

.304 

.022 

10150.0            .273         4.345 

3.834 

2.633 

1.408 

.586 

17450.0            .329         3.154 

2.837 

2.064 

1.215 

.578 

XS   (WITH  CE)=         -999.0 

MIXING   ZONE   LEN6TH=1029815 

.5 

CONC= 

.00 

DIST.   TO  D/S  WPCP=  20000.0 

SHORE  CONC. 

AT  D/S  WPCP= 

2. 
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AVG.  CONC.  AT  D/S  WPCP= 


,79 


CRITICAL  POINT  METHOD  RESULTS 


QY/QT 
.10 
.20 
.30 
.40 

*  *  RUN  NO.: 
QEFL=     75.000 
CEFL=       20.00 

X  EY 


XL 

1348.5 

7364.9 

13483.1 

20000.1 


CL 

6.504 

2.711 

1.461 

.856 


2 
QRUP=  3000.000 
CS=  .02 

.0 


CEA 
.06 
.15 
.27 
.47 


XSCEA 
6536.9 
19898.0 
-999.0 
-999.0 


TEMPR=22.0     PH=  7.0 


.1 


.2 


.3 


.4 


200.0         2.710         7.245         1.266 

.007 

.000 

.000 

2000.0         1.058         3.290         2.248 

.717 

.107 

.007 

4500.0           .963         3.243         2.190 

.674 

.095 

.006 

.0150.0            .461          1.713         1.501 

1.010 

.522 

.207 

.7450.0           .556         1.359         1.216 

.869 

.497 

.227 

XS   (WITH  CE)=         -999.0 

MIXING  ZONE   LENGTH=1085378.6 

CONC= 

.00 

DIST.   TO  D/S  WPCP=  20000.0           SHORE  CONC. 

AT  D/S  WPCP= 

1.; 

AVG.   CONC.   AT  D/S  WPCP=          .34 

CRITICAL  POINT  METHOD  RESULtS 

QY/QT             XL                 CL 

CEA 

XSCEA 

.10          1453.5           2.303 

.17 

7399.9 

.20         8428.7           1.015 

.39 

-999.0 

.30       15934.9              .579 

.69 

-999.0 

.40       30085.1              .290 

1.38 

-999.0 

20 


SUMMARY  OF  RUNS  FOR  MANAGEMENT  OPTIONS 
MISSISSIPPI  RIVER  -  TEST  DATA 


RUN  QEFL  QRUP  TEMP  PH  CAWP  XSCE  QY/QT  CEA  XSCEA  CBKG  CSIJC  CBIOT 


1.  75.000  1000.0  22.0  7.0  .787  -999.0  .10  .06  6536.9  .00  .0  .0 

.20  .15  19898.0 

.30  .27  -999.0 

.40  .47  -999.0 

2.  75.000  3000.0     22.0     7.0     .341  -999.0     .10  .17  7399.9   .00  .0     .0 

.20  .39  -999.0 

.30  .69  -999.0 

.40  1.38  -999.0 


